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ACOUSTIC PARAMETERS 
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DENSITY OF FROZEN MATERIALS 
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SECTION 2 
AIR BLAST ‘ 


Traditionally, the air blast parameter which has 
attracted the most interest is the maximum static overpressure. 
For the typical static o -erpressure vs time profise as measured 
by a pressure sensor, the maximum precsure occurs at the shock 
front, or almost coincident with the arrival of the wave at the 
sensor location. If one is concerned about damage or injury 
from air blast, one must, in addition to maximum static over- 
pressure, be interested in the static overpressure impulse, the 
Maximum dynamic overpressure, the dynamic overoressure impulse, 
and the time of arrival of the air blast shock front as a function 


of distance from the explosion. 


2.1 aretic Environmental Differences 
id The basic parameters of interest in determining the 
ree field air blast values are the pressure and the sound 


velocity, which depends on the temperature and wind velocity. 

As shown in Section 1, the standard pressure for the Arctic is 
essentially the same as the midlatitude vaiue. The temperature, _ 
however, is markedly colder during the winter months. The January 
standard 75° sea level temperature is given as -24°C and inspec- 
tion of Figure 1-4 shows that the mean tempereture is below this 
value for much of the Arctic. The extreme that can reasonably be 
expected is about -57°C. The effect of these decreased temper- 
atures will be noted in Section 2.2. 


oy Temperature inversions are more probable, stronger, 
and more extensive in Arctic than in temperate climates. This 
can enhance the propagation of low overpressure values to long 
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= Wind also affects the transmission of the low over- 


pressure shock wave and causes an enhancement to low overpressure 
damage in the downwind direction. 


| A major environmental difference in the Arctic is the 
high probability of snow or ice cover and frozen ground. Air blast 
over snow can be strongly affected as will be discussed in Section 
2.3. The attenuation of the shock in snow can affect th: covpling 
of the blast energy to the ground or structures. The or.cence of 
the ice layer over the sea can influence the air bDiuc~ €. -.vec 


from underwater bursts. Surface effects wili be discusteu ir 
Section 2.4. 


2.2 | og Free Air Blast Prediction 


Free air blast predictions for nonstandard atmospheric 
conditions are generated as described in EM-l1 (DNA, 1978) from the 
standard 1 kt curves by using Sachs scaling relationships. The 


effects of Arctic meteorological phenomena on predictions will be 
discussed. 


2.2.1 OD sachs Scaling Techniques 


Kod Two basic assumptions are inherent in the Sachs rela- 
tions. First, it is assumed that the air blast wave propagates 
in a homogeneous atmosphere with the ambient conditions at the 


altitude of the observation point. Second, the total energy 
available for air blast is independent of altitude; that is, the 
energy partition is unchanged. 


| The maximum static, maximum dynamic and total pressures 
are related by the expressions 
02 
Po? ea. Pas (2.1) 


where the ranges are given by 
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® is the appropriate maximum pressure, 

?, is the ambient atmospheric pressure, 
Bis the distarce from the explos' on, and 
Wis the yield of nuclear explosio.. 


The subscripts 1 and 2 refer to conditions for the "reference" 
explosion {wsually considered as 1~kt yield at standard sea 
level conditions) and the “problem” explosion, respectively. 


The time of arrival of shock front and the positive 
phase duration are given by 
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where fo is the speed of sound in ambient atmosphere and the 


ranges are related by Equation 2.2. 


The tota. positive phase overpressure impulse and the 
dynamic pressure im.’ se are given by the expression 


Wo\a, 3 2/3 [Cc 
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ty = [— a} I,, 
7 2 Wr Pou Coo 1 
where the variables are as previously defined and the ranges 


are related by Equation 2.2. 


In our application, the subscript 1 refers to the mid- 
latitude standard va.‘«s and subscript 2 refers to the Arctic 
values of interest. The yield will be taken as 1 kt So we are 
interested in the changes that will occur when the 1 kt mid- 
latitude standard curves are scaled to ] kt Arctic conditions. 
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The pressures in the Acctic at sea level are virtually identical -° 
to those found in the midlatitudes. The variations from the 
standard values caused by meteorological perturbations is of the 
same order as for temperate climates. Thus, the pressure ratio 
Po2/Poi is essentially unity, and no differences are expected in 
the pressure radius curves in the Arctic. 


Note that the time and the impulse scaling relations 
also involve the ratio of the sound speed which is related to 
the temperature by the expression 


¢ (2.5) 


where the temperatures must be degrees Kelvin. For the mean 
January Arctic temperature at sea level this ratio is 1.075, 
implying a 7.5% increase in th2 time and impulse values in the 
Arctic. Por the extreme temperature case (~60°C}) the increase 
is 15.5%. 


In Figure 2-1 the change in the shock front arrival 
time is noted for the extreme case. In Figure 2-2 the change in 
the impuise values for the extreme case is shown. Even these 
changes for the extreme case are of marginal interest since a 
15% increase in the impulse would not i1 general cause any 
practical systems effects, and it woul¢ occur with only a small 
probability. The mean 7.5% increase wi ch can be expected in 
the coldest months is within the basic ncertainties in the 
impulse predictions and the resulting damage effects. 
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aay Inspection of Table 1~1 shows a small deviation of 

the Arctac pressures from the midlatitude values as a function 
of altitude. From equations 2.1 and 2.2 the ccaltitude ranges 
° to various overpressure values were calculated as a function : 
of burst altitude for the Arctic and midlatitude pressure - ; 


s ae altitsde profiles. For all overpressure values considered ws 

h , between 1 and 1000 psi there were insignicicant differences (<5%) a 
; in the Arctic and midlatitude coaltitude ranges. . iy 
rf 

ay The conclusion is that no significant differences will . 

’ be found in the free air blast values under Arctic conditions. oy 

Sachs scaling can be used to provide the free air values if Ki 

precise time and impulse values are required. " 


The reliability of Sachs scaling under Arctic condi- 
tions may be questionable, The Sachs relations can be derived ri 
rigorously from theoretical considerations. However, the 1 kt ' 
free air curve is based on 4 combination of theory, calculations A 
and experimental data. For the low overpressure values there has : 
always been some uncertainty. Scaling this curve to conditions 


far removed from the experimental data on which it is based must 4 

; be treated cautiously. c 
ae 

} There is some evidence that Sachs scaling at depressed i, 
‘ temperatures is valid. The technique has been used to correlate : 
i dats in all of the high explosive (HE) tests that have been N 
‘ performed over snow and ice. In the Distant Plain events to be at 
i described in the next section, Sachs scaling was used to correlate 4 
summer and winter results and no inconsistencies were found. “4 

|g Modified Sachs scaling tetween altitudes using the ‘Jy 
4 - . atmospheric parameters at the target location has been used to aa 
‘ correlate and predict blast values in inhomogeneous air with a high : 

. degree of success. Comparisons of computer code calculations in : 

4 ; non-unifuicm air (Wells, 1971) with Sachs scaled blast parameters J 
et 
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indicate that the technique can be used reliably for these cases. 
This would imply that modified Sachs scaling can be used for pre- 
Gicting the blast environment if a temperature inversion is present 
if the pressures are nigh enough to ignore refraction effects. 


2.2.2 OP z- Effect of Temperature Inversion 


Cay A temperature inversion causes a sound speed gradient 

$ xist at low altitudes resulting in refractive effects and 
can, therefore, amplify the overpressure at the ground from a 
burst occurring below an inversion. Converseiy, surface over- 
presssres are reduced if the detonaticn is above the inversion. 
These refractive effects are important only for very low over- 
pressagres (<l psij. The effects are serious enough in consider- 
ation of safety from HE tests, to restrict shots when inversions 
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exist to inhibit long range damage to windows etc. This may be 
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of interest militarily since in the very severe arctic winter 


mn ¢ 
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losing building integrity due to window breakage is much more 


ae 


important than in temperate climates. 


Gay The lapse rates of Arctic inversions are more severe 
than is typical of temperate areas, as described in Section i. 
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It is therefore likely that inversions will exert a more 
significant influence on blast phenomena in the Arctic than 


elsewhere. The increase2 incidence of inversions in Arctic 


2Pe 


areas will increase the vrotcbility of seeing these effects. 


ee 


ae Although corrections for inversions are small, the 
enhancement of low static overpressure at long ranges may somewhat 


OPP RE 


increase the possibility of damage to blast-sensitive targets for 
bursts below the inversion. Later this year a report (Reed, 1980} 
of an extensive experimental study will be puviished detaiiing 

the effect of inversions and wind velocity on air blast. This 
study will supersede anything available at this time. Quanti- 
tative predictions should be delayed until the report is available. 
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In addition to the temperature, the wind velocity 
= causes a change in the relative sound speed and, therefore, on 
. the blast parameters at very low overpressures. No direct effects 
‘ would be expected at higher pressures. The effects of wind will 
be omsidered in a report to be published during 1980 (Reed, 1980). 


a i 2.265 PP = tftect of Wind r 
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‘ 
The dry snow of cold regions is easily lifted by ‘ 
a turbelent winds to create a dense cloud that obscures vision and J 
can become integrated with an air blast wave. Any wind of 
velocity over 15 miles per hour causes blowing snow if the temp- 
erature is well oelow the freezing point. As examples, periods 
during which Slowing snow has reduced visibility to less than 

cs . 1000 yards extend from 75 hours in one area to as long as 266 
consecutive hours in another area. In Sub-Arctic forests, such 


ar a oe a oe eee 


as grow in eastern Siberia, surface winds are impeded by the 
trees and blowing snow is 1l2ss prevalent. 
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qa Tne reduced visibility would have the most direct effect 
on the amount of thermal radiation from a nuclear weapon buict 
reaching the ground. This, in turn, would have an indirect effect ia 
upon the air blast phenomer.a; that is, the possibility of the 
formation of a vrecursor under these conditions would be very 
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remote. 


; | og h nore Significant aspect of the presence of dry snow 

is the fact that a olast wave could carry many snow particles as 

¥ ; it propagates along the surface of the ground (or ice/water sur- ‘ 
faces). This might lead to enhanced damage, which will be 


discussec in Section 2.5. z 


2.2.4 es .:.. Effects of Precipitation, Fog and Clouds “7 
‘ The effects of atmospheric moisture on blast propagatior se 


' are not well known; however, theoretical studies agree qu lita- 
tively with the small amount of experimental data. As a strong 


’ 
blast wave propagates through air containing water droplets it 
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vaporizes some or all of the water. Vaporization of the water 
absorbs energy that otherwise would be available for the blast 
wave to propagate through the air. As a result, the blast wave 
is attenuated more rapidly in air that contains water droplets 


than ém air that does not. 


a The effect of water droplets on peak overpressure may 
be calculated in terms of effective yield. This procedure is 
used to obtain lower calculated overpressures at some distance 
from the burst. Rain or fog has a negligible effect on the 
amount of available energy clese to the nuclear source. The 
energy density within the fireball is orders of magnitude higher 
than the energy required to vaporize whatever water may be 
present, and the amount by which the suspended liquid increases 
effective air density, even under tne extreme conditions within . 
Clouds producing severe thunderstorms, is not likely to exceed 


2 percent. 


Fg Figure 2-3 shows the effective yield for three yields 
and two conditions of moisture content. The water densities 
used in the calculations correspone roughly to precipitation 
rates of 0.1] (light rain) and 0.5 (heavy rain) inches per hour. 


—_ The curves shown in Figure 2~3 are based on the 


assumption of uniform water content between the source and the 
target. In an actual rainstorm, this assumption is artificial, 
but without such an ascumption the analysis of rain's effect 
would be unduly complex. Typically, water content is several 
times as high within a rain cloud as it is below the cloud. 
Actual water distribution patterns are complex, different for 
different rainstorms, and generally unpredictabie,. 
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As stated in EM-1, rain or fog effects should be 
eva d only when the optimization of blast against soft 
targets is important, and then only if the rain or fog extends 
througtcut a volume that includes both the target and the burst. 
HOB curves for thermally near-ideal surface conditions should be 
used with Pigure 2-3 since thermal energy is attenuated by rain 
or fog and precursor effects would not be expected above a wet 
surface. 


| | The effects of atmospheric moisture on other blast 
parameters, such as time Of arrival, positive-phase duration, and 
dynamic pressure are not well known; however, theoretical consider- 
ations indicate that arrival times will remain essentially un- 
changed, positive-phase durations will be slightly reduced, and 
dynamic pressures will be slightly increased. Calculations of 
these other air blast parameters should be made in the normal 
manner, without applying any effective yield factors. Enhanced 
effects on dynamic overpressures are discussed in Section 2.3.6. 
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one can derive some conclusions related to the applications 


Referring to Figure 2-3, and recalling that (¥ 
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to Arctic environments: 
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(1) For light rain or fog, the 125 KT and 1 MT curves 


indicate effective yields of 90% or ebove for <he 
1/3 _ 
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peak overpressures of interest. Since (0.90) 


. ee od 


0.97, it is evident that light rain cr fog is not 
going to cause a significant perturbation to the 
ordinary air blast effects. 


For heavy rain and for peak overpressures in the 
5-20 psi range, effective yields can be in the 
70-80% range for the larger yields. Since 
(0.7)273 = 0.89, it is unlikely that, even for 
this extreme case, the deviations in blast effects 
from normal would be considered significant. 
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No test Gata from nuclear bursts in snow are available 
to the U.S. A possibl2 estimation of the general effect of snow 
can be made by an extension of the reasoning of the preceding 
paxagraphs if we assume that the amount of water in heavy and 
light snows is similar to the amount of water in heavy and light 
rains. The snow particles would have to be first melted and 
then heated to evaporation with the resultant transfer of more 
of the blast energy. This could result in an increase in 
attenmation over that noted in Figure 2-3 since an additional 
energy of about 100 calories per gram of water would be required 
to melt the snow and evaporate it. The interaction may involve 
breaksp of the snow flakes and water droplets for more efficient 
energy transfer. The forve required to shatter the crystaline 
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structure is probably larger, but the effect of this on the 
energy transfer is unknown. There is, however, no positive 
evidence that this reduction snould be greatly different than 
that accasioned by temperate forms of precipitation at militarily 
significant ranges. It should be emphasized here that no valid 
numerical evaluation of this aspect of Arctic environment can be 
made without further experimentation. 


| od Since low dense clouds are very prevalent over the 
polar ice during the summer, the effect might be worth studying 


—_ 


in more detaij. A recent review and analytical consideration of 
this effect (Friedberg, 1976) points out that the attenuation 


in fegs and clouds is more severe because of the smaller water 
drops and more efficient transfer of energy to the water and 
subsequently larger attenuation of blast energy. No work in 
this area was referenced after the 1950s in the above report. 
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2.3 | eae Blast Over Frozen Surfaces 
2.3.1 ww Reflection Characteristics of Snow Layers 


(U)} When a shock front enters a layer of snow it is ateenuas 
ted straagly. Drag forces on the snow crystals dissipate energy 
containe® in the wind behind the shock front. The energy trans- 
mitted tm the snow crystals is then consumed in compacting the 
snow layer. 


— Reflection occurs at the top surface of a deep snow 
layer just as it does at a ground surface. Momentum is conserved 
in the imteraction. A blast wave striking the earth transmits 
only a Small fraction of its energy as ground shock; consequently, 
the earth's surface apprexinates an ideal reflector. A blast 
wave Striking a snow surface is analogous to a ball bouncing 
from a heavy rug. The reflecting surface has a cushioning effect 
that makes it a poorer reflector. 


In the case of a thin layer ot snow, the cushioning 
effect ceases when the pressuce wave ..1etrates the snow layer, 
reflects from the ground surface, and prunagates back to the 
snow surface. At this time, the snow layer iv supported by an 
internal pressure as high as the pressure proa.ced by the blast 
wave reflecting from the surface; the reflecting qualities of 
the snow leyer then approach the near-ideal reflecting qualities 
of the underlying surface. 


Neither theoretical nor experimental date are available 
on the effects of thin snow layers on a blast wave, hcwever, 2 
rough calculation is enlightening. If a shock front in snow moves 
with a speed comparable to that of sound in air, a layer of snow 
one foot thick, struck by a normally incident blast wave, will 
absorb energy from the blast wave for about 2 milliseconds and 
il] have the properties of a near~ideal reflecting surface after 
that time. This 2-millisecond interval is appreciably long only 
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when conpared with relatively short duration blest waves. For 
example, it might alter a 750 psi blast wave from a 1 kt source. 
The overpressure pulses of this blast wave have an effective 
triangalar duration of about 20 millisecends. At lower over- 
pressures, the pulse becomes broader, and ..¢ snow layer wouid 
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have less effect. For a given overpressure, larger yields than is 
1 kt also produce broader pulses. It should also be noted that, ve 
for a ®00-1b HE detonation, the triangular duration of about ¥ 
20 msec occurs at a maximum overpressure of only 20 psi. For HE et 
detonations of smaller charges, these durétions would correspond @ 
to even lower peak overpressures. This discussion indicates = 


the following: 
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= ° If a blast wave with a very short-duration pressure a 

I pulse strikes a thin layer of snow, the snow may e 
“ alter the leading edge of the pressure pulse enough ee 

ome 
j to reduce peak reflected overpressure. The short aa, 
pressure vs time pulse corresponds to high over- ayy 
H pressures from relatively low yield nuclear detonations ny 
7 and/or virtual*’yv all overpressures from small-charge S$. 
] HE detonations. at ng 
ite 
i é ae +.) 
° For a situation where interest is in lower over- ma 
] pressures and yields greater than 1 kt a thin snow atte: 
cover affects such a small portion of the overpressure ] 
F | pulse that peak reflected overpressure is essentially Aa 
the same as for a near-ideal surface. a 
hy A 
; . ental. 
{ , Measurements of the properties of Snow under dynamic ate 
loading have been made (Napadensky, 1964} which indicate that Le 
: : ; ‘ a rar 
| relatively small amounts of energy will be absorbed by a snow = 
layer because the snow is compacted to densities equivalent to Ay 
ry ice by pressures in the 20 - 40 bar region. As one might expect, or 
aa 
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. : 
= a very large variation in snow properties was found for dif- is j 
i 


ferent types of snow in different stages of compaction. The ‘1 
experimemts were not taken to large pressure values so the : Me 
%& integral of PdV cannot be obtained with any degree of accuracy. Hy, 


If the approximation 1/2 PAV is used, which will overestimate 


the integral, then the energy loss due to this mechanism could 4 
be significant in reducing the effective blast yield for a ! <4 
3 1 kt burst detonated over 1 m of snow, which is a reasonable ; + 
upper limit for Arctic winter conditions. | va 
The ©.S. has never detonated a nuclear weapon in the mt 
atmosphere in an Arctic environment. Therefore, all predictions i 
related to the effects of an Arctic environment upon air blast | KN 
parameters from nuclear explosions must be deduced either from . 
theoretical calculations or from the results of experiments y , 
using HE sources. For many years, we have been interested in oy 
predicting the behavior of air blast phenomena from nuclear Ri 
bursts in temperate environments; during that time, these ‘ . 
Geductive methods have proven effective, except for cases where Me 
thermal/air biast interactions are important, e.g. precursor ony 
wave formation and propegatiorn. Experience and advancing devel- 0 
opments in instrumentation techniques have revealed the utility 4 
of and the limitations on the data obtained from the small and aH 
A) large charge HE tests. 
i: 
ay Other than the inability of the HE charge tests to : : 
properly simulate the nuclear bursts' thermal/air blast inter- : ' > 
action, the most important “sin of omission” in HE tests is ; 
insufficient band width of the instrumentation system used. : a 
Sometimes this is referred to as “inadequate frequency response : a 
of the transducer circuitry". In effect, this limited frequency : nf 
response has a Similar effect on the pressure vs time measure- ~ 
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2 ment aS was Gescribed above for the snow layer case. That is, 
the gagz electronic system would "chop off" the true peak over- 
z : pressure and the recorded result would be in error. The magnitude 
; ; of the error depends on the bandwidth of the circuit and the 
, 2 magnitufe of the peak pressure. 


; Fg Recent HE experimental programs have emphasized these 

band wifith aspects; in particular the TRW 3-lb 9404 experiments 

. (Carpenter and Brode, 1974) and the BRL (Dipole West) 1000~-1b TNT 

tests (Reisler et al, 1975 and 1976) employed instrumentation band 
widths which were compatible with the sizes of the explosive 
charge. Unfortunately, the same was not the case for many HE 

( experimeats performed during the 1950's and 1960's. 


aay To explain this concept further, Table 2-1 is presented. 
; The Table lists the instrumentation band widths normalized to 
Carpenter's experiment, which are required to be compatible with 
each Size of explosive charge used for an HE test. It is obvious 
that as the charge sizes increase the band width requirement 


| relaxes. 
1 
] 2.3.2 OD xx Blast Over Shallow Snow 
/ An interesting pair of HE events was conducted as part 


oF the DISTANT PLAIN test series (Reisler et al, 1967). These 

events weve 20 ton TNT surface bursts with the same conditions 

except that Event 3 was a Summer shot and Event 5 was a winter 

shot. The temperature for the summer shot was 110°F and for the 

winter shot was 33°F. The winter shot had a snow cover of ebout 
\ y 4 inches over soil frozen to a depth of about 9". 
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TABLE 2~1 . i 
*REQUIRED" BANDWIDTH VS. TNT CHARGE SIZE eS 
REQUIRED BANDWIDTH 
CHARGE SIZE (TNT) * (NORMALIZED TO CARPENTER'S}) 
1 lb. 800 kHz 
Carpenter 8 lb. 400 kHz 
ae 32 Ib. 252 kHz 
256 lb. 126 kHz 
-o Dipole West 1000 lb. 80 kHz 
Suffield, etc 20 tons 23.4 kHz 
100 tons 13.7 kHz 
£00 tons 8 kHz 
‘ 
\ 
*i.e., These are the bandwidths required so that the data 
system woule be equivalent to Carpenter's system used 
for the 8-lb. exneriments. 
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= Tne comparison of the overpressu.2s obtained on the 
two shots is shown in Figure 2-4. Note that the data prints } 
a@zree closely except at the high overpressure values where the 


the curve for the summer event. There is obviously some scatter 
in the data points, and there are only four gage positions in 
the high pressure region. The interesting fact is that the 
pressure-time records for these high pressure positions indicate 
a very narrow pulse of the order of shock traversal time through 


er 


f 
4 experimenters drew the pressure curve for the winter event below 
j the shaliaw snow while the time width for the stock at the lower 
pressures is significantly ionger than the snow shock transit 
] time. This may be only an interesting coincigence. Additional 


experiments or calculations could resolve the ques#ion. 


: The dynamic pressure and impulse measurements indicated 
fae good agreement between the two events. In this case chere was 
BO « no increase in dynamic pressure due to entrainment of snow by 

7 s the blast wave. 


4 
: : 2.3.3 ae: - Blast Over Deep Snow 


{ ° | Denver Research Institute (DRI) (Wisotski, 1966) and 
U.S. Army Waterways Experiment Scation (WES) (Ingrar, 1960 end 
* . 1962, Joachim, 1964 and 1967) performed HE tests which are 
: . most applicable to our Arctic énvironment situation. The LRI 

; tests employed l-ib and 8-lb size charges, while WES used 32~Ib 
A ; and 256-l1b charges, primarily. In both cases, the band width 
‘g | of the instrumentation was in the region or 0-20 kHz, too limited 
with respect to the size of the sources used. fortunately, most 
of these meaSurementsS were confined to the lower overpressures 
(less than 20 psi) .inere the limited hand width would have less 
effect on the accuracy of the measurements. However, because it 
is difficult to determine the magnitude of the errors due to the 
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limiced band width, one must be very cautious when attempts 

are made to compare data collected by one agency with similar 
data collected by another group using different instrumentation. 
Thus, the most valid conclusions come from the DRI bare ground 
vS snow~corered ground air blast data; more tentative conclusions 
are derived from the WES data taken in the Arctic compared with 
Gata taken over bare ground by BRL at Suffield, Canada. 


ay Of course, we must not forget that all of these con- 
clusions are based upon HE test data; therefore, the implied 


assumption is that the thermal radiation from the nuclear burst 
fireball affects the air blast parameters similarly in both 


hae 
EA 
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temperate and Arctic environments--an assumption which requires 


ove 
we 


Sy 


much more thorough investigation. 


ne vee 


a 


ae DRI performed a series of small-charge HE tests over 
bare ground and over snow-covered ground using the same gage 


tals 
aye 


arravs and electronic instrumentation on each test. These data 
comprise the most complete set of results available on the 
effects of a deep snow layer on air blast parameters. Although 
the charges used were only 1-lb and 8-lbs, since the same 
instrumentation was used for all tests, the lack of sufficient 
band width is probably not serious as far as the overall com- 


parisons are concerned. 


oP e 
5 rethtats 4 


ap The effect of snow and bare ground surfaces on Mach~ 
region peak static overpressures is summarized in Figures 2-5 

. through 2-8. Note that the plotted data are "as read" and they 

Vo correspond to an average ambient atmospheric pressure of 510 mm Hg 

| (9.86 psi). The results shown in Figure 2-5 are typical; the data " 

: indicate that the peak static overpressures for the snow-covered eg 

surface are depressed from those meaSured over bare ground. For 

| the Hc = 1/2 ft case, the two curves are very close to parallel, 
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so the decrease in pressure over snow is independent of pressure 
Magnitude. For He = 4 ft, there is some variation with pressure , 
level indicated; particularly at the higher overpressures where : 
the above-smow curve appears to “turn over" slightly. This 4 
latter behavior is noted also on Figures 2-7 and 2-8 at the 
higher overpressures. Also, the figures indicate that the 


K peak overpressures for the low burst height (He = 1/2 ft) are 
i depressed the most by the snow cover. ! 
ay The effect of snow-covered and bare ground on static ye 

p overpressure impulse is shown on Figures 2-9 thrcugh 2-12 for , 
: the various burst heights. In geieral, the comparisons indicate _ 
: that the snow layer tends to suppress the total impulse; however, . - ) 
Le the scatter in the data is quite severe, and it is difficult to 4 

detect a consistent amount of suppression cue to the different “4 


Surfaces. Looking at Figure 2-12, it is evident that the 7 
variations between the snow-covereé and bare ground values are ; 
reduced as the burst height is increased. : 1 


+ ae ee 8, 
0 st ae as er ty: 


|g The reflection coetficients from snow, bare ground 
’ and concrete are plotted vs scaled charge height in Figure 2-13. a 
( Qualitatively, the results are as expected; one would expect 
that the least amount of energy of the explosive would be trans- 
ferred to the concrete surface and that the most would be 
absorbed by the Snow. BecauSe we are comparing data (concrete) 
taken on another test, using instrumentation witn an unknown 


\ bandwidth, we must be cautious in using the values shown for 
q prediction purposes. 2 
” ; Guy The effect of snow-covered and bare ground on the path J 


of the Mach triple point is shown in Figure 2-14. In general, 
the triple poi.t rises faster over bare ground than over the 
r snow cover. Data from tests having burst heights higher than 


|: ae 


¢ : . esse eOald 2 edit, PSR ee 
A ‘ as ‘ ties i cee mda ha ete ee ae! ea eit tte eee 5 btine athe wi Letalianern ek «nfidlanes. 5 anata dale yAaaare Miettinen tii te! sdanstatticatcel 
4 
1 . ‘a . 
’ se ~ 
“ . 
/ \ Je — as -~ 
“a - 
‘ i . — 
/ oe Fe ig A ; 2 os 


. ¢ , Re wr vr . wr , - am gee AA 
Mier k hd tdi Sd) cee. ta ees 2 ny é ee tale ph Mt, AA rey neh a trata ore a. DD Ro 1.8 O on 5 rata x Re ted 2S arate ah re 


' TTT Tn ere an m NE T RA A CA UR P PUR 
TF RIYUWILT LF OU rw 


: 


' 


‘ ON IMPULSE VALUES 


IMPULSE AS A&A FUNCTION OF BOTH 
HORIZONTAL DISTANCE FROM GRCUND ZERO 
AND CHARGE HEIGHT (Hc) ABOVE GIVEN SURFACE 


BWERAGE BAROMETRIC PRESSURE ~ 510mm. of MERCURY 
ONE POUND SPHERICAL PENTOLITE CHARGES 


. : piss tpethy: Ea adil hoot Geaad satey team cae 228 
SB EES Ge ee’ as Pe cers erode opted sian comet esi oe 
grea T ess eee a iis at 
, SteEisrsts mbt tet | 
¢ apt, at > . oes - 


EFFECT OF SNOW AND BARE GROUND SURFACES | | 


eo wo 


Fk tte ener 


IMPUL SE (psi msec) 
wt 


& 
° P+} -o— —o- He ABOVE GROUND = 1/2 tt B35 
. Phi —-—@ He ABOVE SNOW = 4 ft Seneaae 
T3-O----O-He ABOVE GROUND = 4 ft =ecrss 
. SS Ce cei ee cllohoe yee Sa Ee 
Sadessenusdeniocesseccssss ml feessszenspeteel ttt +] 
‘ 2 aaavecdsecesssase ttririte tte trinttictst ott et 
3 4 5 6 7 6 9W 


HORIZONTAL DISTANCE FROM GROUNO ZERO,. 


; SP cure 2-9. (Wisotski and Snyer, 1966 


| a i ee ee ee eee ee ce A oe eT Oy 


IMPUL SE (psi msec) 


gapmoreneniamen rt ao 


EFFECT OF SNOW AND BARE GROUND SURFACES 
ON IMPULSE VALUES 
IMPULSE AS A FUNCTION OF BOTH 


HORIZONTAL DISTANCE FROM GROUND ZERO 
AND CHARGE HEIGHT (Hc) ABOVE GIVEN SURFACE 


AVERAGE BAROMETRIC PRESSURE - SIOmm of MERCURY 
ONF POUND SPHERICAL PENTOLITE CHARGES 


© 


~ 


tit Hodes 


CG. 


rentftes peese 


trees 
bed e -etberee 


7 1 WEE Saogsecs 
FHI -O——-o- He ABOVE GROUND: | nes =f EES 
FE ~--— He ABOVE SNOW = CO ftp Pere 
+t -| ---O-4ce ABOVE GROUNO : 6 ft ieee ly 
Ss pea 5 OOS SS cnet dacees esencenen So = EGESsG=© 
TSE HUE n isin dies iether em SSCS SE GE 
plaaccctedassace soddcasuscecedesssisrest bso cot et EH 
3 4 5 6 7 8 9 0 15 


HORIZONTAL DISTANCE FROM GROUND ‘a 
GB sve 2-10. (wisotski and snyer, 1966 


Seeicdiieh hese ed keds & Be 
202.47 NN Ye 


Pats 


i 


/ 


GRR e tte oss tee 


b. 


z 


t{MPULSE (psi msec) 


tS a rR aeRO = toe mes —— 
e. 
a “< 
7 a 
“ —_—— 
#08 


= N,N ee Aa er ASRS ASSESS RSS ACAD 


EFFECT OF SNOW AND B4Re GROUND SURFACES 
ON IMPULSE VALUES : 


{MPULSE AS A FUNCTION 93F BOTH 
HORIZONTAL DISTANCE FROM GROUND ZERO 
AND CHARGE HEIGHT (Hc) ABOVE GIVEN SURFACE 


AVERAGE BAROMETRIC PRESSURE - $10mm ut MERCURY 
ONE POUND SPHERICAL PENTOLITE CHARGES 


eb tpe leaves rere mae aT. 
Mm pee eqpeconnes somroeses.ccses coreseens serseeuens soneren i ress 
caboose: aes ies 


poneevess 


Bi 4 
eee 


Heoease! 


Hae 


i -@——@ He ABOVE SNOW = 
Ft.~O——O-— He ABOVE GROUND? 3 ftk&= 


i 


aharieaiy 
' 
i 
| 
HES 
[1] 
0 
age 
Ty 


. Retest =s 

a -~# Hc ABOVE SNOW * 3ft ——— SESEGanes: 
= = tee ip cra +4 
Sasgsisagsissseegian iia = an meeeeeee Deea Sa ree 
pl Soc ca coussscvscscghtsssctiitionfunmnninunment ta aeae sods 
3 4 be 6 7 686 9$ WO & 


HORIZONTAL DISTANCE FROM GROUND ‘a 


mye 2-11. (Wisotski and Snyer, 1966) 
2-29 


6; Pe Ww hi ‘ Heater otaratd ptr Wralh ohn ge Me Me ie te He te te ha te BL Me 
: fe ee Osta stag : fan CHAU CN BOE EY Se CORR R RI re TS 


2 nt ew alee waren cad ete. welll 


tL 
5 FR 


70 


eae = i a ae i 2 


EFFECT OF SNOW AND BARE GROUND SURFACES . 
ON IMPULSE VALUES 
IMPULSE AS A FUNCTION OF BOTH 


HORIZONTS. DISTANCE FROM GROUND ZERO 
AND CHARGE HEIGHT (Hc) ABOVE GIVEN SURFACE 


MMERAGE BAROMETRIC PRESSURE - 510mm. of MERCURY 
ONE POUND SPHERICAL PENTOLITE CHARGES 


bal SERSUizotgotesteae 
Her seret art 


r 
. 


t 


dobre oeet 


i 


trrperts 


IMPULSE (psi msec) 


oe 
BP SESs 5555 Ca sareescornens ro 
Sheets gett iesgecp ee 
4: aa sess: 2 
pas 
. . - 
i : rare eae t 


Heh 
E 


eS 

ae 
Seeesri 
Beesacce 
SiSsasss 
ABOVE SNOW SsScosees3 
ABOVE GROUND = 8 ft Sseesessss 
joasnsds ated seeeerieesteeni gags aoeSencsss 
gas cregeccsss sesss Hestiestinee tire raonne oan ee eee 
4 5 6 7 8 9W 15 

HORIZONTAL DISTANCE FROM GROUND - 


ae AaaeCtne seal amelie ee 


GP cure 2-12. (Wisotski and Snyer, 1966) 
2 


5s spuanennren mcs am ~~ etn e BS Sea Sas ese ae ek en ane PE 
_— + 
. i -\ a 
_ 
Ps ee ery - - 2 a 
ae fe a ie XN aa — A a “ 
~ ee . Ps —_——o _ — i~ 
ee Pe ae 7 BAe, on ae ~~ 
’ ie teem are 
eect oe —_——o =. - _ 7“ 


TT ats ata aimee PU eta Se heen Tat PECs Cee, we dacmete 0s hethine ant scl sembticniee tess alias mmatinne .imrenid neetnin tine ented ARES 


ane See 


: — ~~ TO TT Te a Ce a a 
sa walataF Bo tal tel ae ta! Stal Etat trl Bet TUT ae Dae tot RD WP Ue ate stat De Wi De aE Ae Masia a a aaah 


ee <<< 


e 


*. 


a 


yg tp 
ASS 


2 Mgr ane 


fw 


ad , 
bP ae 


oe 


cM, 


ww Pee. " oD hall 


es BE eae 


hap omg ~- 


AVERAGE BAROMETRIC PRESSURE -S1Omm. of MERCURY 


es 


Were 
Itb 


-O----O- ABOVE CONCRETE (ARL} 


-O- —O- A‘OVE GROUND 


—O-—o— ABOVE: SNOW 


(¥/"@Y)'s - 1N319194309 NOIL93793y 


EFFECT OF SNOW AND BARE GROUND SURFACES 
ON REFLECTION COEFFICIENT 


FIGURE 2-13. (Wisotski and Snyer, 1966) 


is 


SCALED HEIGHT OF CHARGE - Ac,{tt/w”?) 


3 ty7 ox, 


~ 


a PL 


Ie 


ie 


rt tot 5.9 8,8 O.09,8 28k: 


¥ 9,9 gph det gat gat” 


a 
Prien aed 


re 


6 


ee 


4) > OW3Z ONNOYNS 


WOd JONVLSIO IWLNOZIGOH 


QHNOYD 3A08V 
bp 00ed bhobe SEE 0s 205835650) i aoe 


Ss9YVHD ALIVOLN]Ad IvoiWddHuS Od JNO 


ABNIUZN 10 “UAU OIE - 


JUNSS3ud DwAINOUVG 3OVUERAT © 


32 


( 
N 


(Wisotski and Snyer, 1966) 
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3RNNS 3ANOBV LHOIZH 
EFFECT OF SNOW AND BARE GROUND SURFACES 


ON PATH OF THE TRIPLE-POINT 
GURE 2- 
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ty 

are mot definitive and may not follow this trend. Data ee 

on the ground range at which the triple point forms are oe 
incomplete; so, no comparison is possible for the snow-covered 

and bare ground tests. ana 

Syl 

Ga Finally, the DRI data are flotted on a height-of-burst 

(BOB) chart shown in Figure 2-15. The above-snow curves are oh 

supported with more data, and it is possible to be fairly con- . 

wy 


fident as to the form of these curves. The bare~-grouna data 
are less extensive, but it is again evident that the snow 
reduces the distance at which a particular peak overpressure 
is observed. The magnitude of the distance reduction appears 
to increase as the overpressure level decreases. There is 

7 strona evidence of the over snow contours "pulling in" for the 
Ad surface burst case (HOB = 0); this is consistent with the fact 


nen Co] ‘meal ont =. Cd Lec a f ae 


that a saerface detonation over snow loses a large portion of 
its explosive energy to the snow which is close to the explosion. 


7 The Greenland HE series involved a large number of 

* tests from about 1958 to the middle of the 1960s. A large number 
; of WES and Cold Revions Research and Engineering Laboratory 

] (CRREL) reports which were referred to previously were written to 


describe the results of the various tests. Included were tests 
over and under the deep snow on the Greenland ice cap, and over 
and under ice. Shock transmission through snow and ice were 
measured as well as a large number of cratering shots in snow 
and ice. A report never widely distributed summarizes these 


results (Smith, undated). 
The HOB related shots were primarily 32 and 256 i 

pound charges with scaled heights of burst to 12 £t/1bt/3, The we 

instrumentation band width was too narrow to adequately resolve > 
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the narrow pulses; so, as is the case with the DRI experiments, 
one must be very sareful in comparing the WES data with other 
data. I. this case, however, no comparison measurements were 
Made over ground with the same instrumentation; so the compar- 
isons are more uncertain than with the DRI measurements. 


2.3.4 | Overpressure Contours from HOB Tests 


For the military planner, the air blast height-of- 
burst (BHCB) charts are the most useful for prediction purposes. 
Since the Arctic environment data we have for eir blast is from 
HE tests, we shall emphasize the HE HOB charts; also, maximum 
overpressure is the principal parameter we shall consider. 


fd A series of high explosive (HE) blast tests was con- 
ducted jeintly by the U.S. Army Ballistic Research Laborat ries 
(BRL) and the Canadian Defense Research Establishment Suffield 
(DRES) daring the fall of 1969. These tests, held at the 
Watching Hill test range at DRES in Alberta, Canada, were known 
as the 1969 Height-of-Burst Series (Reisler et al, 1976 and 
1969). Later, during the summer of 1975, another series of HCB 
tests ‘as conducted by BRL aS a part of the three-year DIPOLE 
WEST series (Reisler, 1975). 


qa Some of the results from these HOS tests are plotted 
in Figures 2-16 through 2-18, showing the peak overpressure 
contours for various overpressure values (Reisler, to be pub~ 
lished). These data correspond to air blast wave propagation 
over bare ground under "near-ideal" corditions, which implies 
that there tre no significant thermal effects. 


Looking at these figures, the plotted data and the 
sOlid-line contours correspond to the BRL tests referred to 
above. Additional curves are shown to correspond to data 


collected by other agencies on their tests using various HE 
i 
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OVERPRESSURE HOB CURVES, 2 and 1 PS] 
(Reisler, to be published) 
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| : = sizes; tests were performed by Sandia Corporation (SC) 


(Vortman & Shreve, 1976), the Naval Ordnance Laboratory (NOL) 
(Hartman and Kalanski, 1952), and the Atomic Weapons Research 
Establishment (AWRE), United Kingdon (UK) (Worsfold, 1957 and 
1963). The BRL contours indicate that there is some data 
scatter around the actual contour lines; as is usually the case, 
the data scatter is more pronounced for the lower overpressure 
contours. It is also significant to note that HOB data from 
other agency tests do not always agree with the BRL curves. 

In fact, for overpressures of 10 pr:'i and lower, the deviations 
are significant. For comparison purposes, we shall use the BRL 
contours, but we should remember that an error band of + 10% is 
estimated for the data. 


The data plotted on Figures 2-16 through 2-18 are 
"as read", and although they are scaled to 1 1b TNT, they are 
not scaled to sea level conditions. The atmospheric pressure 
at the test site varied from about 13.38 to 13.87 osi. The 
pressure scaling factor (S,) for this test series varies from 


about 1.060 to 1.098. “This means that the correction to sea level 


conditions would be between 6% and 10% for the data shown. 


|g Data from both the WES and DRI HE studies have been 


combined in Pigures 2-19 and 2-20 to show how the data over 

Snow compare with the BRL bare-ground HE data. It should be 
noted that the smali-charge data have been Sachs<scaled to BRL 
average pressure P. = 13.63 psi. As was discussed in some 
detail in Section 2.3.3, such data comparisons can be misleading, 
if taken too literally. This is because the WES and DRI data 


were obtained by using instruments with inadequate freguency 
response. Therefore, it is likely that a portion of the 
obvious displacements of the over-snow overpressure contours 
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_ the ware ground contours is due to the limited bandwidths, 


ard it is difficult to determine what portion of each displace- 
ment is “aeal®. The conclusion is that there is an effect, 

shown qualitatively in Figures 2-19 and 2-20; however, to attempt 
to quantify that effect based on the data available, will prob- 
ably lead tt larger effects than actually exist. 


2.3.5 GP vicrs Scaling of Saow Depth Effects 


| acd The minimum snow depth on the various DRI HOB measure- 
n 


S Over snow was about 6"/1bt/3, If this snow depth is 


1/3 relation, then these EOS 


scaled to muclear yields by the W 
curves for a 1 kt would correspond to snow depths of at least 
sixty feet, which is much deeper than snow encountered in the 
Arctic exoept for the snow/ice depths found in the highly 


glaciated areas. 


The DISTANT PLAIN winter event snow depth of 4" is 
equivalent to a depth of about one foot when scaled for a kt. 
The typical snow depth can range up to 60 cm to 1 m near the 
end of the winter season over much of the Arctic region. Thus, 
we are left in a quandary. The HOB curves over deep snow show 
a marked drawing-in of the curves for surface bursts over deep 
snow with mo dependence on snow depth, while the surface burst 
over shallow snow showed no effect or at most a questionable 
effect at Sigh overpressures. 


Phere is no real reason to expect a priori that the 
standard w/3 scaling should be used when considering surface 
interaction effects due to the snow which is far from an ideal 
reflecting surface. For an ideal reflecting surface with no 
energy loss at. the surface or for near-ideal situations where 
only minor effects are expected then the wi/3 
justified. 


relation can be 
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Measarements of the response of snow to loading 

(Napadensky, 1964) indicate an elastic response at overpres- 
: sures below 10 to 30 atmospheres depending on the snow type, 
: then a crushing region where a large volume decrease occurs 
with small increases in pressure, then a region with relatively 
small volume decrease as the pressure increases to 160 atmos- 
pheres or so until the density of ice is approached. Thus, 


for pressures below the yield threshold no permanent deforma- 
tion of the surface would result. 


The snow surface does not act like a rigid boundary 
even in this elastic region. In Figure 2-21 (Ingram, 1962) 
the magnitude of the reflected shock measured over a snow 
surface is compared with the theoretical value over a rigid 
surface for normally incident shock waves. The values of the 
incident shock are considerably less than the yield strength of 
snow. Note that the measured shock pressure is about 70% of the 
theoretical value and the difference seems to be increasing at 
the higher overpressures. No data were given for non-normal 
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incident shock waves. These meaSurements were taken in Greenland rw 
with 100 fcot snow depths; so extrapolation to shallow snow cases ie, 
is uncertain. The DRI experiments involved snow depths as small iy 
j as 6"/1b2/3, The reduction of the pressure over snow as com- " 
pared to bare ground was about 11% averaged over ull ground j 
ranges and burst heights. The DRI bare ground values were less os 
than the rigid surface values as indicated by Figure 2-13, where ! "4 
the reflection coefficient for ground is less than for concrete. | . 
. No calculations are available to indicate the depth of snow nyt 
: required to induce these effects as a function of yield and t 
' specifically to indicate the magnitude of the effect expected ; , 
' for the nuclear case. in 
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For incident pressures above the yield limit, PV 
work is done by the crushing process and energy is removed 
from available blast energy. Porzel (1962) gives 
Q= 1/2(P-P iV -V) as an estimate of the energy absorbed 
by an ideal absorber which will overestimate the energy 
absorbed. If we use Py as 150 psi or about 10 atmospheres and 
(Vo-V) of 2 for compressing snow of density of about 
-3 g/cm then we get the follow.ng estimate of the energy 
absorbed by a snow layer. The energy loss as a function of 
range is given by the expression 


R R 
AE AE 
SE f Am dam 24D f sn r ar 
Ro Ro 
R 
= 1.04x102p f (P-Po)e dr (2.6) 
Ry 


the integral can be evaluated from the 1] kt standard pressure 
radius curve. If the fractional energy loss is considered and 
if yields other than 1 kt are allowed we have 


rywl/3 
AE __D -10 AP rar 
We cI7s * {1-04x107? f 73 73} - 42e7) 
Ry pyl/3 


7 and the ranges are 


where D represents the snow loading in g/cm 
in cm. The integral has been evaluated from RS corresponding 
to the charge radius, and the expression in the braces is shown 
in Figure 2~22. Beyond the range corresponding to 150 psi the 
ey we is zero; the value of the braces is essentially 

7x10 ~. 
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Reduction Factor 


1G 10 104 
Overpressure (psi) 


GR risvre 2-22. ENERGY REDUCTION FACTOR DUE TO SNOW LAYER 


Thus, one might expect if the above assumptions 
are correct that the reduction in yield for overpressures 
below about 150 is given by 


D -4 
= x 7x10 . 
wi73 


The ranges to these overpressure values might be given by 
scaling by the expression 


AE 
Ww 


R'(p) = (w- AE)?/3 RB, VA (P) (2.8) 
rather than R(P) = w/? Ri ql?) (2.9) 
so that R‘/R = (1 - 48)1/9, (2.10) 


Consider the HE charges over deep snow. The snow 
depth was at least 6" scated to 1 pound charge. Therefore 
D = 6.6 g/em” and AE/w ~.4. Therefore, R'/R~(1 - 24373 = .64 
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: > pie reduction in range of about 20%, which is of the order 


of the changes noted in the experiments. In practice one might 
.- expect that the ranges would be depressed more for stations 
closer tm the ground and less at the higher altitudes, whereas 
the above estimate is an average reduction assuming that the 
blast ware is developing symmetrically from the burst point. 

' Details of the interaction at the surface such as the effect 

a of angle of incidence of tne shock wave have been ignored. 


Recall that the experiment showed no effect of snow 
depth for snow depths considerably larger than the 6" scaled 
minimum. Making the same calculat‘ons for the 20 ton HE shot 
with a maximuin snow depth of 4" or about 3 g/cm? gives 

~s AE/W ~7.7x107> or essentially no reduction in yield and no 
reduction in the pressure-radius relations, confirming the 
experimental results. 


Note that the above relation does involve a wi43 


scaling of snow depth. Extrapolating to the nuclear 1 kt case 
and a snow depth of 1 m or a loading of 30 g/en* we obtain 
AE/W = .0@21 or a negligible effect. The effect would be even 
smaller for larger nuclear yields. The above general agree- 
ment may, of course, be fortuitious and a thorough theorctical 
investigation of the subject considering the air shock inter- 
action with the nonideal surface should be made. 


2.3.6 QB era Effects and Precursors 
al Observations on the low-altitude nuclear weapons 
e over bare ground show that at a thermal exposure level 
of 10-30 cal/em* a popcorning effect occurs where particles 
of the soil are forcibly ejected into the air. This apparently 


' occurs due to the very rapid heating and vaporization of the 
; water entrained in the sand (or other) crystals in the soil. 
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- ejected particles are heated and form a very efficient 


mechanism for heating the layer of air for a few feet above the 
surface. A similar effect occurs when rapid heating of orgenic 
materials takes place on the surface. The natural convective 
heat tramsfer will also be very high and will assist in heating 
the air layer. These types of effects are certainly strong 
enough tc lead to the formation of a precursor wave. 


The precursor is characterized by a highly turbulent 
flow behind the wave front. Dense dust clouds raised by this 
turbulence tend to follow the shock front az it propagates 
outward. 


a No empirical evidence is available to indicate the 

effect of the thermal and shock environment from a nuclear 
burst over snow. The following assumptions have been made in 
determining the effect of the thermal pulse on snow. First 
the energy is assumed to be deposited in the top centimeter of the 
snow layer. This thickness is arbitrary and the thermal energy 
is undoubtedly transmitted deeper than this in new light snow 
and to shallower depths for old packed snow. The actual depth 
is not critical; however, the point is that very high tempera- 
tures that would be obtained by assuming the energy to be 
deposited in a very thin surface layer are not realistic. 
Secondly, it is assumed that any melted snow is not heated 
above the melting point because of the very high conductivity 
of the slush that will result from surface melting. This means 
that the energy contained in the thermal pulse will result in 
melting the maximum depth of snow possible instead of raising 
the temperature of the melted snow. Of course, if the snow 
melts completely, the temperature of the surface may begin to 
increase above 0°C. 
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; ar: characteristics of snow cover a wide range. The 

° ivity can vary from .5 tc .9, depending upon the condi- 

a tion of the surface so that the absorptivity may vary from .5 
» ‘- to .1. Fresh snow, then, will reouire about .9 cal/em? deposited 
to reach the temperature of 0°C and another 8.1 cal/em* to melt 
each centimeter layer for a total of 9 cal/en* for each centimeter 
of snow depth. Since only about .1 of the energy is absorbed, an 
incident exposure of ahout 90 cal/em* will be required for each 
centimeter of depth. Assuming packed dirty surface conditions, 
: the required exposure is about the same since the density and the 
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. absorptivity can increase about a factor of 5 each. 


i The above estimates indicate that about 2700 cal/em* 


vs would be needed to ccmpletely melt one foot of snow. No mechanism 
is availabie to transfer the energy to the air. This is far above 

‘ the 30 calfem? of tnermal energy that typically will produce 
popcorning and other surface effects which serve to transfer 
energy to the air layer. The conclusion from this discussion 
is that under most arctic environments, conditions will not be 
favorable for the formation of a precursor ' "ast wave; that is, 
the thermal/air-blast interaction effects will be minimal. This 
contlusion may be substantiated by experimental measurements 
being performed presently in solar furnaces (Knasel, 1980). 


2.3.7 GRP: 2c 2vence of Snow and Water on Dynamic Pressures 


The air blast dynamic pressure is defined by the 
relation 1/2 § v’, where 5 is the density of the air behind the 
: shock front and V is the particle velocity of the air. Experi- 
t ment has shown that blast waves which are “loaded” with dust, 
es e@.g., precursor waves, can produce higher-than-expected damage 


. 
, 


to drag-sensitive targets. 
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cal the explanation is that the dust picked up by the 

last is accelerated to near shock-front velocities, and the 
increased average density of the air/dust combination result; 
in enhanced pressures. 


It is expected that the same would be true to some 
extent for the Arctic environment; however, in this case the 
blast wave would be loaded with ice crystals and/or water 
particles. The net effect would be similar to the dust case 
with density and dynamic pressures increased. In order to 
Getermine the magnitude of these increases under various 
conditions, thorough investigation is needed; some data are 
available from blast waves propagating over water. Other useful 
information could be obtained from --omputer code esults. 


2.4 We... Blast from Underwater Bursts 

= The air shock resulting from an underwater burst has 
een measured on a few underwater nuclear bursts and several 
series of small charge conventional explosives tests. 


2.4.1 P= arison of HE and Nuclear Tests 


| Al Chapter 7 of DASA 1200 gives analytical techniques for 
¢ uting the air shock expected from underwater bursts for sev- 


eral DOB, which take into account the available empirical evi- 
dence. Prediction curves are given to show the expected air 
shock for a 1 kt nuclear burst for a wide range of DOB. 


are A series of 5 ton HE tests were made (Pittman, 1970) 
to determine the air blast from underwater bursts and to cor- 
relate with the sparse nuclear data available. Very good 
correlation with the Baker and Umbrella nuclear data was obtained 
by using the water column or plume velocity as the scaling para- 
meter for shallow bursts. No correlation of the air blast 
: effects with cavitation closure was possible. - 
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omy NOL hasS a program to compute the airblast from under- 
water bursts by using two-dimensional hydrodynamic techniques, but x 


results have not been released for publication (Lorenz, 1980). 
No calculations including the effects of an ice cover have been a 


: made or are planned. vy 
i ~ 
: 224.2 Ber rece of Ice Cover ne 
No experiments have been done to determine the effect ) 
of an ice cover on the air blast from a nuclear weapon. Con- Nt 
sideration of the air blast production mechanisms described in Ny 
DASA 1200 lead one to expect, if anything, a decrease in the air i 
shock if an ice cover were present. It does not appear that an & 


increase in the air blast due to an ice cover could occur for 
equivalent DOB as compared with an underwater burst. 


Contribution to air blast arise from three different 


mechanisms, the relative importance of which depends upon the 4“! 


DOB. The initial air pulse results from the transmission of 
the water shock across the interface, another contribution arises oe 
from the spray dome, and the third from the plume. 


the water shock into the air 
1/4 


The direct transmission of 


is the dominant mechanism only for depths below about 700 W 


feet where the spray dome and plume effects are minimal. In \ 


this region the water pressures are iow enough that acoustic hn 
| theory can be used to provide an estimate of the coupling at the 
interface. DASA 1200 explains several techniques of varying 
complexity to descrihe the energy transfer across the interface 
and propagation into the air. The expected air shocks are very 


weak (< 1 psi). 


Replacement of a layer of water with ice at the sur- 
face would result in a decrease in the coupling efficiency 
because of the introduction of a second interface where mis- 
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matching and energy loss can occur. Using the values of the 
ice, water, and air acoustic characteristics given in Section 
1.2, we can estimate the size of the effect as follows: 


fg The overpressure in the air is given by the expression 


A : 2 F Cc. cos ce 
4p, Pa c. cos ¢,, + Ry Cy cos ?, 


(2.11) 


where a and w subscripts refer to air and water values of the 
Parameters, P is the censity, C is the sound speed and ¢ is the 
angle from the normal to the wave front. The angles are related 
by Snell's law: 
Sin %, 


sin 3) = (2.12) 


For simplicity consider normal incidence, then substitute values 


for parameters and we have AP /AP = 5.6 x 1074 wi.ich indicates 


the reason why such small air blast occurs with deep bursts. 


Ld If we have an ice layer between the water and air then 
e have 


iz 2x428 |, 2#2.95x10° 


2.95x10°  1.54x10°+2.95x10° 


= 3.8x1074. (2.13) 


fnerefore the effect of the ice layer is to reduce the air blast 
pressure by about 1/3. 


se] The spray dome results when the water shock pressure 

s strong enough when it reaches the surface that the resultant 
tension in the water from the combination of the reflected tensile 
wave and the incident compression wave exceeds the tensile strength 
of water. This results in cavitation and the separation of a 
layer of water from the surface with some imparted upward momentum. 


The spray dome then produces an air shock which can be predicted 
by the technigues noted in DASA 1200. 
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7 Introduction of an ice layer for an equivalent layer 
of water on the surface would obviously cause changes in spray 


dome development. The pressure pulse transferred to the ice and 
reflecting as a tensile pulse at the upper ice surface could 

lead to ejection of a layer of ice whenever the tensile strength 
of ice 1s exceeded. Since the tensile strength of ice is much 
larger than that of water, this will occur only for much larger 
values of water shock pressures than are needed for spray dome 
development. This probably will imply a smaller value of air 
shock than produced from the spray dome. If the censile strength 
of the ice is not exceeded, no air shock from this type of 
mechanism would be expected. 


The plume or water column is the dominant air blast 
mechanism when the DOB is less than about 75w/3 ft. The plume 
is treated as a supersonic hody moving through the air, and the 
air shock is computed as described in DASA 1200 by standard hydro- 
Gynamic considerations of the bow shock from a blunt body. At 
the deoths where this mechanism is important the water shock pres- 
sures are so large (104 psi) that a considerable thickness of ice 
would be shattered. If the entire thickness were shattered, the 
effect of the ice on plume development would probably be similer 
to an increased DOB equivalent to the ice thickness. If the ice 
layer were not completely shattered, then some of the energy of 
the plume would be expended in breaking up the ice layer and the 
air blast would be expected to be less. 


In the above considerations, the effect of the ice 
cover, if any, would reduce the magnitude of the air blast. It 
is not expected that more detailed calculations involving hydro- 
dynamic considerations would change these qualitative conclusions. 
Detailed calculations wouid be necessary to determine safe escape 


ranges for aircraft delivering for instance an ice penetrating ASV 


nuclear burst. 
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2.5 Gp =<: Couplin to the Surface from a Low Altitude Burst Be KY 
| ar The coupling of energy into the surface from a low a oy 
titude burst is obviously very intimately connected to the - oy 
cratering problem which is considered in Section 3 and aiso is = ii 
related to the air blast HOB curves which are considered in : .] 
Section 2.3. . s 
i 
2.5.1 ae Ground Coupling Effects i a 
ae fwo cases are of interest involving a snow-ice-ground - ‘! 
configuration. In the first the burst occurs above the snow oe i 
layer so that the shock must traverse the snow layer to reach : is 
the underlying ground or structure. In the other case @ burst ! 
occurs below the snow layer as might happen with en impact fuze br 
which is not actuated by the less dense snow lays’. In the first W 
case the snow layer will act as an attenuating medium and will on 
reduce the energy transferred to the underlying medium. In the Mi 
second case a tamping action might occur and an increase in a 
energy coupled into the underlying material may occur. « 
a 
Both the WES and DRI HE test series included shots in ~ 
snow with an attempt to measure shock wave parameters in the snow nN 
as well as the movement cf the snow (acceleration, velocity and + 
displacement). A common problem of these measurements was a Ry 
very large scatter in the data as evidenced in Figure 2~23 ow 
(Wisotski, 1966) and Figure 2~24 which shows the bounds for i 
the data points for shock measurements in ice and snow (Ingram, a 
1960). The long dashed lines in Figure 2-23 are the limit lines b 
for the snow data from Figure 2-24. The two sets of data are me 
seen to be in essential agreement and suffer from the same order on 
of uncertainty. The source of the data uncertainties include ae 
possible guenching of the charge by the snow surrounding the " 
charge and the difficulty of getting good coupling between snow & 
and the gages since snow is a mixture of air and suspended ice , 
, crystals. ; " 
mh 
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PEAK PRESSURE AS A FUNCTION OF BOTH 
HORIZONTAL DISTANCE FROM GROUND ZERO 
AND CHARGE DEPTH IN SNOW LAYER 


MEASUREMENTS TAKEN AT MiD-CEP" 4 OF SNOW LATER 


AVERAGE BAROMETRIC PRESSURE - SIOmm of MERCURY 


ONE POUND SPHERICAL PENTOLITE CHARGES. 
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FIGURE 2-24 
PRESSURE IN ICE AND SNOW FROM HE EXPLOSION 
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lr showiag that more attenuation of shock energy is taking 
has been suggested (Smith, 


Place. A deray of pressure as p38 


The slope of the snow curve is much steeper than that 


undated) as being a reasonable fit of the snow shock measurements. 
The DRI measwrements at pressures less than 1 psi show a marked 
reduction of the slope, but the curves drawn to represent the 


data are very subjective. 


ee Note that for pressures well under the yield limit there 
S appreciable attenuation of blast energy. Calculations of the 
attenuation of blast energy by precipitation referred to in Sec~ 
tion 2.2.4 considered energy transferred to water droplets and 


resulting in vaporization of the suspended water for overpres- 


sures as low as 13 psi. The cutoff pressure was essumed to be 


a function of water droplet size but independent of water concen- 


tration. However, the highest concentration considered was about 


5% by weight. Scaling of these results to a snow density of 


23 g/cm results in attenuations much larger than noted in Fig- 
ures 2-23 and 2-24. Friedberg considered evaporation of the 


water requiring about 700 cal/g. It is possible that the shocks 
in snow involve melting of the snow, which would require about 


80 cal/g and which might occur at lower overpressures since 


smaller temperature rises are involved. Then, however, one 


should ask why the shock in ice shows no indication of attenua- 


tion. If the energies involved are large enough to involve phase 


change effects then an attenuation in ice shocks would be expected. 


If we assume attenuation is due to energy lost in 


crushing the snow, then the effect can be estimated by using the 


same general procedure as in Section 2.3.4. The energy lost up 


to a range R is given by 


R R R 
AE = At am = 4mp Sf SE rar = 2.073x10°% f aPr’dr. (2.14) 
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Or a particular yield W (kt) we have if 
‘ 
“ 
= = 2.073x10 ° p f ap(r/w'/3)? a(rswi/?) (2.15) ne 
R 3 
ofall m1 
U 
Wi 
a 
a 
nt 
where AP is the overpressure minus the yield strength. If My 
the above expression is evaluated for a 1 pound charge, then the 
dash~-dot line in Figure 2-24 is obtained. The shock falls oy 
‘ 
progressively lower than the free air curve until the assumed Hy 
yield strength of 140 psi or about 10 atmospheres is reached “ 
st 


then parallels the free air curve. This is of course only a 
very crude estimate of the effect, but again it is interesting 
that it is in che range expected. 


P| In the nuclear case we do not have a burst in a large 
amourt of snow, but are interested in the attenuation of the 
blast wave crossing a depth of snow of order of a meter or less 
in thickness. The above calculation shows that the energy losses 
in a spherical case scales as wi/3, This would imply distances 
about 126 times larger for 1 kt than for the 1 pound HE charges 
and would indicate that the snow depths normally encountered in 
the Arctic would have essentially no effect on the coupling. 
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ey This is to be expected if the shock energy density is 
c red. A one thousand psi shock wave has an areul energy 
2 cal/cm?. The energy loss per gram of 


density of aout 5x10 
snow is 1.65a107° AP or 1.65 cal/g for 1000 psi. Since the 
snow loading is of order 30 g/em* 


insignificant. Of course we are assuming that no PV work is 


at most, the energy loss is 


done for presSures below the yield strength so the attenuation 

to low presszre blast waves would be zero. This does not agree 
with the experiments which do show attenuation as compared with 
the air curve. The stress-strain curves of Napadensky mey not 

be accurate at low overpressures and there may be no well defined 
yield point as he measured. Unconsolidated snow would be expected 
to have a very low yield strength. At the present time a quanti- 
tative measure of the protection of the snow layer is not possible 
but the effect is expected to be small for typical Arctic snow 
depths for nuclear yields. 


feat al The possible tamping acticn of snow if a burst is 


detonated below the snow Layer has been considered by Science, 
Systems & Software (Allen, et al, 1975). In Figure 2-25 

the results are shown for a snow depth of 6 g/em*. At 6 psec 
there is abozt a 15% enhancement of energy coupled to the ground 
and tne energy in the air is somewhat less for the snow case as 
would be expected. The calculations were not carried out to 
later times but the difference might well disappear by later 
times. However, note that the snow loading is considerably less 
than the 30 g/em* that can be present in the Arctic. A larger 
coupling efficiency might be found at lower yields. In practice, 
the snow layer above the burst would he perturbed which weuld 
tend to reduce the tamping effect. A sample calculation with 
deeper snow should be made to later times to determine the magni- 
tude of this effect even though a large effect is not expected. 
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rst bas been considered experimentally as well as theoretically. 


2.5.2 GP wacer Coupling Effects 


The coupling of energy from a low altitude or surface 


However, there is not a large amount of data in this area and 
certainly none thac considers the complications due to an ice 
layer. For low altitude bursts where the coupling of the air 
blast imto the water is of interest, one would expect the 
presence of the ice cover to decrease the shock transmitted 
into the water because there are two surfaces with inpedance 
missmateches instead of only one. 


Re ey For near surface bursts where there is interaction of 


the weapcn outputs with the surface, the situation is much nore 
complicated. There were several nuciear weapon tests involving 
very swall heights of burst over sea water in the Pacific. How- 
ever, the weapons were mounted on barges in the tests. The area 
covered by the barge was large enough to have a strong effect on 
the coupling to the sea water. Fer this reason, any underwater 
shock measurements in these tests would probably be different 
than for a burst directly over the water. 


Systems, Science and Software has performed a series of 


calculaticns to determine the early time coupling of energy from a 
1 MT burst to various surfaces. The coupling of energy to sea water 


was compared to that with NTS Tuff (Allen, et ai 1974). At very 
early times the energy in the sea water is about 50% higher than 
that in Tuff. The calculations did not conrinue to late times to 
consider the underwater shock formation and growth. The increase 
in coupling was due to the lower opacity of sea water as compared 


to soil. The presence of salts in sea water does affect the opac- 
ity. The salinity of sea ice is less than sea water but is highly 
variable depending on the ice history. Because of the vast energy 


available and the high temperatures that are reached, one would 
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expect the ice calculations to be very similar to the water 
calcultions. Radiation-hydrodynamic calculations of the 
subseqvent shock development would be necessary to determine 
the effectiveness of this method of coupling energy into the 
water pressure pulse as compared to an underwater burst. 


The presence of snow cover on the ice canopy could 
affect the coupling of energy to the ice then into the water as 4 
discussed in Section 2.5.1. The magnitude of the tamping action 
versus snow depth and yield is unknown. This effect could have 
implications in ASW. If a technique for locating Soviet sub- 
marines under the ice is developed, then the necessity of using 
an ice penetration weapon must be addressed. In this case the . 
coupling efficiency for the various ice surface configurations 
will be of great interest. 
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There is much uncertainty connected with the underwater . 


R shock from near surface nuclear bursts even if the ice cover is 
A not present. A large dependence on the details of the surface 
why configuration may exist even for the late time air and water 
“ shocks. In the vast there has been little incentive for work 
< in this area. Increased Soviet use of patrols and the current , 
‘ nonavailability of ASW techniques for the polar area may result 
% in an interest in these problems. In order to develop effective 
4 airborne tactical nuclear ASW techniques it will be necessary 
RX to consider the surface effects on underwater shock. j 
" te 2.6 GP ix prast Target Damage Sffects 
Kr et The two most important environmental effects on targets 
‘ or target response in the Arctic are the snow cover on targets and / 
iN the temperature of the materials used to build the targets. A . 
uJ possible effect is an. increase in dynamic pressure due to snow , 
‘ loading of the shock wave. ¢ 
" | 
: a 2.6.1 , gies Cover on Targets 
% It is a fact of life in the Arctic that target struc- 
y tures, even those built above ground, will be covered with a 
i layer of snow and/or ice. In fact, most structures designed 
4 for arctic use are built to take the most advantage of this 
‘* ic cover layer. Snow cover over suriace or buried structures 
A affords protection to the structures because it attenuates the 
a : air blast load transmitted to the structure. Air-blast-induced +: 
‘tg accelerations in a snow layer from detonations above the surface 
. attenuate rapidly with depth. Peak vertical downward accelera- 
: oe tions at 2 feet below the snow surface are 3 or 4 times greater : 
i \ than those at 5 feet. Much of the air blast energy is absorbed 
v in compac..ng the snow layer. 
f 
Por structures and equipment above ground, the most 
‘ al effective snow cover protection is afforded by a snow berm over 
XK the top of the structure. This berm eliminates any corners 
‘ ‘ fc 
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or vertical walls and presents a smooth aerodynamic surface 


to the air blast wave; this has the advantage of preventing 5 
* large reflected pressure loading of the structure. Also, as ft, 
= the snow berm becomes somewhat compacted, it can contribute r 
p re to the overall structural strength of the target. 
é : 2 7 Por underground (vr undersnow) structures and equipment, & 
‘ the snow cover, in addition to providing attenuation of the shock hi 
A 7 loading, again contributes to the structural strength. This 
‘4 e L strength contribution can be traced to the "bridging" effect of ‘ 
‘ the snow arch over the buried structure. In.a sense, this snow \ 
: , arch acts as an additional structural member when a load is 
: applied. . 
a . i oa The protection afforded by the snow cover is, of course, ' 
a functi:.. of the geometry of the snow cover in relation to the ri 
construction of the target in question and is dependent on the Ky 


properties (density, moisture/ice content, etc.) of the snow 
coyer vs depth. Therefore, it is not possible to present useful ~ 
generalized predictions of the effectiveness of snow cover pro- 

tection; each case must be considered individually. 7 


; , Ee During the Greenland HE test series, the resistance " 


¥ . of snow arches was considered (Smith). A summary of 
the results obtained is shown in Figure 2-27. As one would 


wa 


expect, the damage level depends primarily on the ratio of y 
arch span to the crown thickness. A strong word of caution. 

is necessary because these were 4E tests, and the width of the 
: pressure pulse is much less than would be experienced from 

( nuclear tests at the same overpressure levels. No calculations 
have been made predicting the magnitude of this effect. 
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2.6.2 Effect of Target Temperature 


Those structural materials which are used on targets 


i 
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essentially tne same in Arctic climates as they are elsewhere 
May react to cold in such a manner as co alter their vulner- 
ability to nuclear effects. Srecifically, metals, rubber, 
plastics, ceramics, and fabrics will undergo changes in their 
strength, elasticity, impact resistance, and other related 
characteristics. These changes will increase the susceptibility 
of the material to damage. Steel is an important material for 
military targets. The mechanical properties of steel vary with 
temperature in a non-uniform fashion. However, the most important 
effect of lcw temperatures on steel is to reduce its ductility. 
This property change can cause brittle fracture to occur in 
structures exposed to relatively small static loads. A signif- 


icant reduction in impact resistance will also accompany a loss 
in ductility. 


= During the series of HE tests conducted by WES on 
e Greenland ice cap, some military equipment was inadvertently 


exposed to air blast loading. However, during those tests, no 
measure of the loads and/or response of these targets was obtained. 
Thus, without any definitive data on the subject, we could only 
speculate on the quantitative effect of the reduced temperatures 
of the target materials with regard to damage criteria. This is 

a technical area which requires more thorough investigation. 


2.6.3 Enhanced Dynamic Pressure 

It has been speculated that air hlast dynamic pressures 
in the Arctic would be "loaded" with ice crystals and/or particles 
for many situations of military interest. These waves impinging 
upon drag-sensitive targets could impose enhanced forces, which 
would result in more severe damage than one would predict for 
the unloaded waves. 
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| aceral Virtually no pertinent data are available pertaining 

o this effect; before one could attempt to quantify the effect, 
7 a great deal of effort would be required to collect data and 

! perform computer code calculations to check the data consistency. 


2.7 Brvoncrusiors and_ Recommendations 


wf In the previous sections the current status of knowledge 

of air blast and surface effects predictions under arctic con- 
ditions were considered including the free air blast parameters, 
precursor effects near the surface, the effects of precipitation, 
clouds and inversion layers, changes in the height~of-burst curves 
over shallow and deep snow, and surface coupling considerations 
including the attenuation or possible tamping effect of a snow 
layer and the effect of the snow/ice canopy on shock transmission 
accross the air-water interface. The uncertainties in the various 
subjects can be corrected by a recommended research program. 


2.721 Qi corcs ions 


The cold temperatures in the Arctic cause a slight 


increase in the time of arrival, time duration and impulse 
expected from a free air burst. If one were considering the 
effect of attacking a specific impulse sensitive target in the 
coldest area of Siberia, the change might be worth including. 
The overpressure-radius and dynamic pressure-radius relations 
are unchanged since the atmospheric pressure and the variations 
in the Arctic are essentially the same as in temperate climates. 
In conclusion, the free air prediction values given in EM-1 

are adequate for Arctic free air values, and sealing to arctic 
pressure and temperate values is not necessary. 


The attenuation of blast weve energy by precipitation, 
fogs, and clouds is considered in EM-1 and has been treated in 
later studies. The amount of precipitation and the precipi-~ 
tation rates in the Arctic are in general less than in most 
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=. areas. The fact that the precipitation will likely 


be snow will not increase the magnitude of the effect. Light 
rain and fogs reduce the effective blast yield of at most 10% 
even for large yields at overpressures as low as 1 f51. For 
most studies the attenuation could be ignored since it is very 
@ifficult to have accurate knowledge of precipitation patterns 
and rates. 


fg If a burst occurs below a temperature inversion 


refractive effects can cause a focusing action and increase 

the extent of blast effects along the ground for very small 
overpressures (<1 psi). If the burst occurs above the 

inversion layer, the opposite effect is noted. The high 
probability of strong inversion layers in the Arctic would 

lead to an enhancement cf these effects noted in temperate 
climates. Since the effects are only noted at very small over- 
pressures, the military effects of inversions could be important 
only for low overpressure targets. Inversion effects should 
definitely be considered in determining fai.-safe ranges for HE 
testing. An experimental program concerning the effect of temper- 
ature inversions and wind on blast has recently been completed 
and a definitive report on this subject will be puplished during 
198@. Calculations have been proposed to determine whether 
temperature inversion effects can occur for overpressure values 
as high as 1 psi. 


fg In EM-1 the recommendation is made to treat frozen 
ground, snow, and ice as thermally ideal surfaces and, there- 
fore, not to expect any classical precursor effects as dis- 
cussed in EM~1. This is a resuit of the very large amounts 
of energy required to produce water vapor and heat up a layer 
of air near the ground and the fact that these surfaces will 
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: ustally have a large albedo implying absorption of a small . i 

fraction of the incident thermal energy. Experimental con- i 

firmation of this fact is expected in the very near future. .. | 
J 
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Frozen soil, ice and Snow samples will be exposed in the French 

solar furnace to determine their response to thermal loading and 
i their capability of transferring heat to the near surface air 
layer. 


q lg Bwen though snow is expected to be a thermally 
ideal surface, because of the strong attenuation in snow, 
Significant reduction of air blast cver snow has been measured 
for bursts over deep snow. Reductions in ranges of 25 - 40% to 
the overpressures less than 100 psi were noted in these experi- 
ments. Scaling the snow depths to nuclear yields would result 
in depths of about 50 et/kti/3, Depths this deep would only be 
§ found in Greenland or other highly glaciated areas. However, 
the validity of scaling the snow depths with yield in this 
Manner is very questionable. No theoretical work has been done 
in this area to determine exactly what interaction is occurring 
in the snow layer and to determine the proper scaling method. 


Comparisons have been made of the air blast from 20 ton 


: surface shots over bare ground and with a snow depth of 4", which 

: scales to 1.25 et /kt?/3, No change was noted in dynamic pres- 
sures on impulse. A possible reduction was noted in the over- 

fh pressure over 600 psi over the snow. This difference was not 


explainable and could be due to data uncertainties. 


a> Thus, we have two sets of data, one with scaled snow 
Gepths of 50 £t/kt?/3 
overpressure contours, and another with a scaled snow depth of 
1.25 £t/kt?/3 


where large differences were noted in the 
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where no change was noted for a surface burst. 
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The typical arctic snow depth will usually be less thin 3' which 


. is in the shallow scaled regime. However, there is no assurance 
s that this type of scaling is valid. 


call the presence of deep snow can affect the coupling of 

: ast energy in two ways.’ In the first, if a low altitude air 
burst is used to attack a hardened structure either above or 
below the ground covered by a layer of snow, then one would 

F expect a decrease in the blast energy coupled to the structure 

: and a decrease in the damage. Large numbers of measurements of 

: shock in snow from HE shots have been made. The scatter in the 

data range oyer an order of magnitude. The experimental uncer- 

tainties are large and involve an effect on the HE burning due 

to the snow and the difficulty of getting a good match between 

‘ the snow and the measuring instruments. There have been experi-~ 

ments performed to measure the basic shock properties of snow 

and the data again show a very wide scatter depending on the 

State of the snow. No calculations of the attentuation to be 

expected from snow layers have been located. Theoretical pre- 

dictions have been made that snow shock values are very similar 

to NTS Tuff. At the present time, no predictions on the attenua- 

tion properties can be made. 


—_ In the second case, if a weapon were detonated below a 
snow layer (for example by having an impact fuze that does not 
actuate in snow) the tamping action of the snow because of the 

= larger opacity as compared with air could result ‘an ae larger 
coupling of energy into the ground. Likewise, if a weapon were 
detonated on the surface of the snow, the decrease in opacity as 
compared to ground micht result in a larger transfer of energy 

into the srow and ultinately into the ground than for « ground 

; , . surface burst. A single calculation of the tamping effect showed 
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— a 15% increase in the (6 wsec) ground coupling for a 1 M7 


burst under & g/cm” of snow. Considering that a typical Arctic 
ease involves about 5 times this amount of snow, one might 
experience a significant increase in the grounded coupling. How- 
ever, until calculations are extended to later times, no pre- 
Gictions of the magnitude of the effect are possible. 


Batth of these snow effects could have implications for 
attacking tazgets in northern USSR. In the first case, the snow 
May tend to decouple the air blast energy from the target and 
lead to less damage and ettack effectiveness than expected. In 
the second case, the snow may enhance the damage and attack 
effectiveness. The second case may have implications in ASW 
also. Currertly, one desires a burst at a sufficient DOB so 
that little energy is dissipated above the water surface to 
maximize the submarine damage range. This would require an 

ice penetrating weapon in the Arctic. However, because of the 
snow and/or ice tamping effect, it may be possible to fuze the 
weapon to go off under the surface of the snow or ice and 
enhance the coupling of energy to the water so that an under 
water burst may not be required. 


2.7.2 GRP ze connensat ions 


Sicnificant uncertainties which may be important for 
systems in t... Aretic were found to exist in the following areas: 


o The effects of precipitaticn, fogs, clouds and 
temperature inversions on the air blast 

o HOB curves over snow for nuclear yields 

o Effects of snow cover in altering the coupling of 
energy to the surface versus HOB/DOB 

o Air blast from underwater bursts through an ice canopy. 
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‘ Note that alf of the above items involve shock propagaticn through 
: and interaction with lossy materials consisting of air mixed with 
: quantities of water in various states (vapor, liquid or solid). 

: No recent hy#rodynamic calcuJations were found considering these 
materials. fsolution of uncertainties in all of the above areas 
‘ could be obtained by a three part research program. 


Preliminary Paalytic. 4 aud 1~Dimensional Hydrodynamic Calculations 


During thie phase equation of state information should 


be collected an snow, ice and frozen ground materials. Analytical 
- calculations wsing th: developed theories of shock propagation 
through lossy materi .ls should be made to determine the attenua- 
tion of shocks throus:. these materials and to provide confirma- 
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tion of hydrodynamic runs. A series of led hydrodynamic calcu- 
lations shoulé be made addressing the attenvation of shocks in 
air-water mixtures and snow, the coupling of shocks from air 


<a 
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x to the ground arc structures through various snow depths, and 
to compare the response of frozen grounds with rocks. These 


ee em 


calculations aight prov: e resolution of some of the uncer- 
tainties in the above areas and would provide guidance in 
setting up multidimensional hydrodynamics runs. 


Specifically the effects of precipitation, fog and 
clouds in causing attenuation to the air shock could be deter- i 
mined and compared with current calculations and predictions. : 
The runs showing the coupling of air shock through various ‘ 


ae 


' snow depths into the ground and structural materials will 
show the degree of protection provided by snow cover. These : 
calculations should be done for various incident shock strengths : 
to show the effect on both very hard targets such as silos ; 


(1000 - 2000 psi) as well as softer structures (<100 psi). The 
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ated shock angle should be varied to see what effect this 


has on the coupling and air interaction process. This will 


provide guidance in the effect of the snow layer on the air 
blast and effects expected in HOB studies. 


ad The results of the phase 1 calculations would be used 
Oo determine what 2-d hydrodynamic runs are necessary to resolve 
the remaining uncertainties. 


a A series of runs may be necessary to produce HOB curves 
over snow. The effect of the snow depth on the air blast may 


be yield dependent ard it may be necessary to generate curves 
for more than one snow depth. 


|g Calculations of the shock transmitted to hard targets 
covered by a snow layer from a low altitude air burst may be 
necessary depending upon the results of the l-d coupling and 
attenuation calculations. 


The tamping effect of a snow layer should be deter- 


mined by repeating the calculations of s? 


for depths of snow 

\ representative of Arctic conditions. If a coupling significantly 
greater than 15%$ is noted at the early times, then the calcula- 
tions should be carried to later times to determine the increase 


in the ground shock. 


The tamping effect for smaller yields representative 
of ASW weapons should be determined in a snow-ice-water geometry 
to determine if ice penetrating weapons would be required in 
attacking submarines beneath the ice. 


PWigven:ige . 


' 2-74 
aay | 
ie 
h 
. aes eiitreretban. eommmmms erm on + a, 5.3 saat 7 
me Aileen De chanare” mm mn pants, mene amrcterbiatis eitatimmens 2 satsectitel ensabmend deamreinai ehh ncadankalaie 5: hake saben ia EN eS , 


| nd Depending or the results of the phase 1 calculations 
~d calculations of the air blast from underwater bursts with 


an ice canopy may be warranted. If ice under strong shocks 
loses its integrity then these calculations will not be neces~ 
sary, and the ice cau be treated as an increased equivalent 
water layer for air blast predictions. 


Experimental Program 

| During phases 1 and 2 of the program requirements 
experiments to define the basic physical properties of 

snow, ice, and frozen ground can be determined. Information 

in this area exists, but in the ten years since these experi- 

ments were performed better techniques have been developed, 


Depending on the results of the gomputer calculations 
a series of HE tests in an Arctic environment may be warranted. 
The subjects of interest would be effect of depth of snow on 
air blast measurements at the surface and above the surface, 
effect of inversions, correlation if any between yield and 
depth of snow, effect of ice canopy on water and air shock from 
underwater bursts and coupling of airshock through the snow to 
the ground for a burst above the ground. Such a series should 
include static overpressure and dynamic pressure versus time 
measurements and have an instrumentation system with an 
adequate band width to resolve the narrow pulse widths. 
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SECTION 3 
CRATERING PHENOMENA 


The mechanisms producing a crater for near surface, 
surface or subsurface bursts are closely allied to the air 
blast and surface effects considered in Section 2. There have 
been no nuclear tests by the U.S. in cold climates; so the U.S. 
has no experimental data base for nuclear cratering phenomena 
in the Arctic. As discussed in EM-1 the data base for nuclear 
craters consists entirely of large yield bursts in the Pacific 
and small yield bursts in Nevada. Thus, cratering from nuclear 
bursts is a very uncertain subject at best. Adding the complex- 
ity of Aretic conditions increases the uncertainty. 


3.1 GR arctic Environmental Differences 


The difference of importance in cratering is the large 
Obability of occurrence of snow, ice and frozen ground in 


the Arctic. In heavily glaciated areas the snow/ice thickness 
will be deep enough that the entire crater forms in these 
materials. Most of the area, however, wiil have only 1 m or less 
of snow or ice over frozen ground; so a layered geometry must be 
considered in the cratering predictions. For large yields the 
scaled depths of snow or ice are negligible, and as will be shown 
later, the crater in the ground will be little affected oy the 


snow layer. 


The ice canopy may influence the underwater crater 
development. The existence of underwater permafrost may be 
important. Adequate experimental data ire not available in 
these cases. 
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3.2 | Craterinag Mechanisms in Arctic Media 


In the Arctic environment, it is obvious that the 
medium in which the explosion crater is created can take many 
different forms. Some of the forms are bare ground (frozen 
and/or underlain with permafrost), frozen ground covered with 
snow and/or ice, thick ice layers over water, and shallow bodies we 
of water. In the latter case, the crater could form in the solid 
medium under the water. 


| Pigure 3-1 is a schematic illustration of a crater <a 
1 
ormed by a surface burst, showing descriptive nomenclature. Ss 


Since crater size varies primarily with charge yield, depth 


of burst (DOB), and the cratered medium, it is desirable that 
tests be conducted with as many different charge geometries 

and in as many different media as possible. This also involves 
the development of suitable scaling relations by which results 
of small-scale tests can be used to predict the results to be 
obtained with much larger yields. Thus far, attempts to 
correlate theory with empirically developed exponents, or 
scaling laws, have met with only limited success. 


Figure 3-2 shows some ideal crater cross sections from 
nuclear bursts, illustrating the effect of HOB and DOB on crater 
volumes. If a nuclear or HE burst is sufficiently high above 
the ground surface, only a Shallow compressional crater is formed 
and no ejecta produced. As the height of burst decreases, the 
crater volume increases and an increasing fraction of the crater 
is due to excavation and ejection of material from the crater 
recion. 
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Fg Tae large difference in energy density (ratio of 
explosive yield to explosive mass) between high explosive énd 
nuclear devices can cause substantial differences in cratering 
efficiency (the ratio of volume of crater to explosive yield) 
and in the relative importance of various cratering mechanisms 


between the two types of sources. 


The material properties of the medium in the crater 
region influence the crater volume primarily through their 
compressibility and shear strength under dynamic loading con- 
ditions. Water content plays a large role in determining shear 
sirength, especially in soils. The largest crater volumes are 
found in wet soils and the smallest crater volumes are found 
in rock. Jointina is an important factor in determining the 


1 ey ae 


crater size in rocx geologies, Geologic layering is a rough 
indicator of material properties with depth and must be con- 


anes’ 0 


sidered when predicting crater volumes. Frozen ground or 
ground interspersed with ice lenses and/or permafrost will 
behave like rock as far as crater formation is concerned. 
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‘ rg The failure process in snow differs from that in 

3 glacial ice, frozen ground, rock, and certain types of soil. 
, Characteristic features of this failure (referred to here as 
ys “viscous-damping failure”) are: 1) damping of the disturbance 
during the rise to peak pressure, and 2) substantial recovery 
of stored potential energy during unloading. Due to the 
unique physical properties of snow, craters formed by explo- 
sions in snow will be unusual in appearance and size 

compared with craters formed in other media. Snow is a 
composite material that consists of a relatively incompres- 
sible crystalline solid (ice) and a compressible gas (air). 
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—_ is found in the interconnecting voids in the ice matrix 
and comprises up to 79% of the volume of snow near the surface 
of the ice cap. Other properties of snow of importance in 
cratering are low melting and vaporizing temperatures. 


Immediately after detonation, as the hot gas bubble 
begins to form a cavity by vaporization, the surrounding snow 
is compacted radially, and the air in the voids is compressed. 
Cavity walls are fractured and an ice skin is formed by fusion. 
During this loading of the snow, a significant amount of the 
explosive energy is expended in compacting and deforming the 
snow withogt destroying cohesion. Some snow is dissociated 
and thrown cut as ejecta. 


Much of the energy used to compress the air during 
loading is recovered during unloading (after the pressure 
wave has passed), which resuits in fracturing and deforming 
the snow. Tne primary cavity then exhibits a reversal in the 
direction of displacement (implosion) as the snow attempts to 
regain its original location. This part of the mechanism is 
referred to as pseudo-elastic rebound. Simultaneously, the 
compacted snow zone and the ice skin are fractured. 


fg A sensitively balanced transition condition appears 
to exist zt critical depth. The balance determines under what 
conditions fractures during the rise of pressure and the outward 
expansion of the gas bubble predominate over fractures formed as 
a result of implosion. Implosion is closely followed by a vortex 
within the snow and scouring action as the gas bubble emerges 
from the rising column definec by the vortex, This scouring 
largely determines the final shape of the apparent crater. 
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At a charge depth less than that at which maximum 
scouring occurs, more of the energy of the explosion is expended 
in the atmosphere and less is available to the snow. An apparent 
crater is formed in the cir blast range and the secondary zone 
of the fragmentation range. Refer to Figure 3~l1. The volume 
of the apparent crates per pound of explosive charge is max- 
imum at the transition limit between the two ranges, where 
scouring is a maximum. Dimensions of the apparent crater are 
neither predictable with accuracy by conventional cube root 
sealing nor usable as a basis for predicting undersnow damage 
because 1) it is difficult to determine the proportion of the 
explosive emergy partitioned to loading tne snow, and 2) the 
apparent crater in snow occurs subsequently to loading and is 
the result of the scouring action of the vented gas bubble. 


3.3 


High Explosive Cratering Experiments 


Tere have been many HE cratering experiments performed 
in the Arctic or sub-Arctic. The surface materials include snow, 
ice and frozen ground of various types. Many of the experiments 
have been designed to cetermine the optimum depth of burst of 
various HE types and charge sizes for producing the largest 
crater for mining and excavating. In Figure 3-3 (Bauer et al, 
1973) representative cratering efficiencies are given as a 
function of depth of burst for several arctic materials. For 
Purposes of nuclear cratering emphasis on shallow or surface 
bursts would be of more interest. 


The large diffezences in the lower and upper limits 
for frozen materials noted in Figure 3-3 are typical in crater- 
ing experiments due to variations in local geology and material 
properties. The cratering efficiency cf HE charges increases 
with increasing water ccntent. As discussed in EM-1 this has 
also been noted in unfrozen ground materials. 
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fg In Figure 3-5 the scaled radius of snow and ice craters 


In conjunction with the WES and CRREL air blast experi-~ 
ments ix Greenland, extensive measurements (Livingston, 1970 and 
1968) were made of craters in deep snow and ice. The depths 
were such that the entire crater was in the snow or ice. Shallow 
and evem above ground burst heights were used in addition to 
depths of burst extending below the optimum depth of burst. 


In Figure 3-4 the efficiencies of HE cratering in ice 
and snow as a function of the scaled DOB are compared. The 
very wide bounds shown in the figure result From several dif- 
ferent types of explosives and several charge sizes. Thus, 
the Scatter is due to material properties and differences in 
efficiercies of explosives as well as possibly to inapplica- 
bility of ~ube-root scaling of the charge weight. Note that at 
optimum depth of burst the snow crater will be about 3 times as 
large as a crater in ice. The optimum depth of burst in ice is 
somewhat deeper. For @ surface burst the snow crater is about 
twi e 4s large as an ice crater. No data were provided for 
near surface air bursts over ice. 


are compared as a function of scaled depth of burst. The radius 
of the snow crater is much larger than ice especially at the 
deeper depths. For a surface burst the radius for snow is about 
50% larger than for ice. In Figure 3-6 the scaled crater depths 
in snow and ice are compared. Tne differences for deeper depth 
of burst are not as large as for the radii but for a surface 
burst the crater depth in snow is about twice that in ice. 


Surface bursts are very important militarily; so the 
Greenizrd surface burst experiments have been analyzed as a 
function of charge weight (Conway and Meyer, 1970). The apparent 
crater depths and radii are summarized in Figure 3-7; also 
included on these figures are data from Sager (1960 and 1961). 
Figure 3-8 shows the apparent crater volume as a function of 
charge weight. Other cratering data from surface events in 

6now are virtually non-existent. 
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Figure 3-9 shows the variation of apparent crater 
radius with charge weight for surface bursts in snow as compared 


1 with craters from surface TNT events in clay, sand, basalt and 
4 shale. Similarly, Figure 3~10 presents a comparison of apparent 
. crater depths versus charge weights. These figures show that 


eraters in snow tend to be larger than craters in other media 
for the same charge yield. 


. a This increased size appears to be due to the greater 
amount of material vaporized and compacted during the explosion. 

i Although no ejecta measurements were made, examination of the 
crater lig profiles indicates that the contribution (to volume) 
of the ejection mechanism in snow craters is correspondingly 
less than in craters in other media. Craters in snow have a 
characteristic wide shallow appearance. The magnitude of the 
pseudo-elastic rebound in snow is greater directly under the 
charge than in the material pushed laterally outward because 
of the greater lateral confinement of the material under the —~_ 
cnarge. 


3.3.1 OP scarin Considerations 


Equations for scaling crater dimensions in snow within 
a range of yields of 06.5 to 5,000 Jbs, as determined by the use 
of the method of least squares, are presented in Figures 3-7 
and 3-8. These equations show a significant departure from the 
common cube-root scaling. For the apparent crater radius, a 
Slightly smaller scaling exponent of 0.26 is indicated. The 
scaling exponent for apparent crater depth, 0.15, is considerahly 
i smaller than that normally applied to craters in soii. These 
unusual scaling exponents are probably best explained by the 
mechanism ef pseudo-elastic rebound in snow. A correspondingly 
‘ [ low scaling exponent of 0.75 is evident (Figure 3-8) for the 
apparent crater volume. It should be noted that these empirical 
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‘ 
scaling components are based on a limited amount of data and ae 
should be con-idered as approximations. Tne use of these ‘ 
exponents tc scale HE data to nuclear expiosions would be zy 
questionable because of the magnitudes of the NE yields and the 
differences in thermal enersy re2lease, which appears to influence iS, 
crater formation in snow significantly. oy 

Even though the scaling rules for the snow radius and si 
depth are very uncertain when extended to nuclear yields, it is rn 
instructive to compare these results with the wet soil EM-1 : 
pr2dictioms. Using the relations in Figures 3-9 and 3-10 a 4 
crater fran a 1 kT surface burst over snow would have a radius x 
of 127 ft and a depth of 19 ft. Using the relations given in ‘ 
the revised cratering section :£ EM-l, a crater from al kT 4 
burst in wet soil would nave a radius of U1 ft and a depth Ay 
of 42 ft. Thus, a wide shallow crater is predicted by the EF ht 
snow data when scaled to nucl:ir yields. i" 

j 

Snow efficiency varies from about 7x10" to 3x104 £t? ra 
per ton, being equivalent to wet sand or muck. Ice varies from oy 
4x10> to éx10? £t? per ton. Frozen soils range from 3x10 to y 
6x103 £t3 per ton. The highest efficiencies are found in frozen os 
silts which are equivalent to wet soft rock, and the lower ie 
efficiencies are for frozen aggregates which are equivalent to Mt 
hard rocks. As is normal for cratering measurements, a wide At, 
range of values is noted. It is suggested that for want of a a 
better method these efficiencies be used in conjunction with Cy 
the prediction methods in EM-1. . 
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- df The effect of the increased thermal yield of nuclear 
ursts aS compared with HE is unknown. For typical soil materials 

° small yield muclear devices (<1 kt) are assumed to be about 1/7 

; as efficient as HE and large yield devices (21 kt) are assumed 

ee to be 1/20 as efficient as HE charges. Thus, a 1 kt nuclear 

° burst would produce a crater volume of about 3x10° to 1.5x10° ft? 


e per kt. A kt of energy is capable of melting about 1.5x10° et} 


: of snow (p = .3 g/cm", melting energy about 80 cal/q) and vapor~ 

} izing about 1.7x10° et? of snow (vaporizing energy about 700 cal/g). 
If a fraction of this energy were available to increase the crater- 
: ing efficiency for nuclear bursts, then the efficiencies obtained 

' i using the current prediction methods may be too low by as much as 


i . a factor of two. 


Pi 3.3.2 Gp .: ered Geometry Considerations 


i see Consider the 19 fi. crater depth found for a 1 kT burst 
in snow. This is much deeper than typical snow depths except in 
highly glaciated areas. Thus, the cratering data considered 
above must be modified to include the effects of the shallow 
snow. HE experiments considering a layered geometry of dry 
soil over wet soil have been represented by the expression 


Pe eles wa AL os etme IF ek el ea Ree 2 lee SS er OA he eth 8 ee . 


(VV, /(Vy-Vp) = 1 - exp(-5.4 a/vl/3) 


where : 

a = depth to base material (water table or cemented : 

layer) ; 

Vv = apparent crater volume in the layered geology : 

Vy = apparent crater volume in the surface material j 

: when d = © ' 

Vy _ = apparent crater volume in the base material when : 
d= 0. 


In Figure 3-11 the data for a sand over a cemented soil layer is 


oh Rey Nae eR ae, 


shown. The upper curve will represent a case with a definite 
boundary between the surface and base layer such as would occur 
for snow over frozen ground. The dashed line is the equation 


given above. 
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iets Im using the technique given above an iterative 
i he volume appears on both sides of 


process is mecessary since t 

the equation (or on both axes in Figure 3-11). The volumes Vy 
and Vy are obtained by using tie HE efficiencies along with tie 
proper nuclear yield and efficiency ratio factor. When this is 
done for the typical Arctic case of 1 m of Snow over frozen 
ground then the following observations can be made. For 1 kt 
Vu = 1x10 £3 and Vy coe be as low as 2x104 et?. The value 
of V would be about 2.5x10° ft. The scaled snow depth is only 
about .05. In this case the volume of the crater in the base 
material would not be altered appreciably by the thin snow 
layer. For larger nuclear yields this depth of snow wouid be 


insignificant. 


ae The calculations for the lzyered geometry are very 


uncert.ain, and we are applying the results of layered geometry 
eases .ar outside the original configuration. 


3.4 Underwacer Cratering 
Urderwater cratering is discussed as a major topic in 


Chaptéc 2 of DNA EM=-1 and in Chapter 8 of the Underwater Hand- 


book. iLowever, existing manuals do not discuss any effects that 


a 


may be caused by conditions peculiar to cold-weather regions, 


nor do there appear to have been any experimental investigations 
into this matter. The discussion that follows must therefore be 
regarded as con;ectural. The factors that might cause variations 


under Arstic ccoaditions rrom what is predicted under temperate 
conditions are differences in bottom composition, if any, and 


the presence of ice. 


|g Reference to Table 2-12 in Problem 2-36 of DNA EM-1, 


wherein soil correction factors are given for various bottom 
materials, reveals that the range of bottom materials covered 
encompasses the range of materials expected to be found in any 
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of the world's oceans, including the Arctic. It is known, how- 
ever, that subsea permafrost exists in several of the seas of 
the Arctic, the Laptev, Kara, E. Siberian, Chukchi and Beaufort 
Seas for example (Lewellen, 1973, 1974, 1977). The extent of the 
substa permafrost is not known, and there i3 a limited amount of- 
infozmaticn as to its characteristics (Chamberlain et al., 1°78). 
Of particular importance is whether the permafrost in a given 

Ke area is bonded, with ice in the interstices, or unbonded, 
saturated with brines that have a depressed freezing point. 
The soil correction factors may well be different for the two 
States, and different rreom that applicable to the basic material 
forming the permafrost, 


Gi Underwater cratering from an underwater detonation 


occurs when the fiist expanding bubble interacts with the bottom. 
The presence of an ice cover would affect the early time history 
of this bubble only if t bubble interacts with the ice layer 
as well as with the bottom. .nd then only to tne extent that the 
energy required to vaporize ice differs from that required to 
vaporize seawater. Depending upon whether the ice is old ice 
of low salinity, or more recently-formed ice of higher salinity, 
the latent heat of fusion may vary from less than 40 cal/g to 


i! 


% the 89 cal/g of ice of zero salinity at -1°C {Neuman and Pierson, 
1966). The energy used to vaporize ice is thus some 6-12% more 
than would be used to vaporize an equivalen: zmount of seawater 
at 0°C. 


wal In analogous fashion, for a surifice burst or low air 

urst over water to create an underwater crater, the expanding 
fireball must vaporize the water layer beneath and interact with 

% the bottom. Again, the preseice of ice cover may require 6-12% 
more energy for vaporization than if no ice were present. Except 


j in extremely shalicw water, the volume of ice ton be vaporized, 
even in the case of a relatively thick solid ice pack, would be 
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a Small percentage of the water to be vaporized. The decrease 
in exergy available for cratering would therefore be expected 


to be extremely small in most cases. 


It is concluded that, except where the bottom is subsea 
permafrost, the methods of predicting underwater crater dimensions 
i gives in Problem 2~36 of DNA EM-1 are valid under Arctic condi-~ : 

tions, regardless of the amount of ice present. The uncertainties ; 
°- in crater dimensions given are of the order of plus 150-160% to e 
minus 50-60%. These are large enough to encompass any additional 
uncertainty due to ice cover. In very shallow water, if the 
crater lip height were such that it extended above the water or 
ice {unwashed crater), cr to just below the water surface, the 
scouring action of broken ice caused by its wave-induced motion 
would be expected to hasten the erosion of the crater lip. If 
the bottom is composed of subsea permafrost, the proper soil 
correction factor to usSe is not currently known. 


3.5 GR conciusions _and_ Recommendations 
The following conclusions are drawn from the rather 


“8 Cant information applicable to arctic cratering and recommenda- 
tions of research necessary to reduce the uncertainties are given. 


3.5.1 Conclusions 


Equations for scaling crater dimensions in snow within 
a range of yields of 0.5 to 5000 lbs TNT show a significant 
departure from the customary cube-rcest scaling. For apparent 
crater radius, the scaling exponent 0.26 is indicated, which is 
close to fourth-root scaling. The scaling exponent for apparent 
erater depth is 0.15, which is considerably smaller than that 
usually applied to craters in soil. A correspondingly low scal- 
ing exponent of 0.75 is derived for apparent crater volume. 


oo 


It should be noted that these empirical scaling expo- 
nents are based upon a limited amount of data and should be con- 
sidered as approximations, The use of these exponents to scale 
up to nuclear explosion craters would be questionahle because 
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‘ magnitudes of the nuclear burst yields and the pr. fou } 
differences in thermal energy release, which appears tc infiuence 

h erater formation in snow significantly. When comparisons of } 
i ; apparent crater radii for TNT tests in snow are made with 

a i craters from TNT tests in clay, sand, basalt, and shale, the i 
é \ comparisons show that craters in snow tend to be larger than } 
\ craters in other media for the same charge size. This increased t 
« size appears to be due to the greater amount of material 


. vaporized and compacted during the explosion and to a scouring 
+" _ action. ¥ 


= Although no ejecta measurements were made in the WES 


pe ee 


4 
tests, examination of the crater lip profiles indicate that the 
* contribution of the ejecticn mechanism (to volume) in snow 
« craters is correspondingly less than in craters in other media. , 
4 Craters in snow have a characteristic wide, shallow appearance. : 
. The magnitude of the pseudo-elastic rebound in snow is larger ' 
§ directly under the charge than in the snow pushed laterally ; 
oe outward, because of the greater lateral confinement of the , 
z snow directly under the charge. ; 
& : 
‘" Craters in frozen soil or permafrost have similar ! 
‘ appearance and dimensions to craters formed in hard rock. It . 
is speculated that the mechanisms for the formation of these 
a craters are similar to those in other soils. 
. Since all of the available cratering data obtained ‘ 
a under Arctic environmental conditions have been collected from 
’ HE charge tests, the main questions which remain when one uses 
MY these data to predict craters from nuclear weapon bursts are: 
e 
2 1. Does tne enhanced thermal radiation associated wicn a 
- the nuclear burst have a profound effect on the i) 
zs partition of the total energy going into the snow ' 
« and/or ice cap? ' 
£ ‘ 
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a Z. Further, does this enhanced thermal radiation 


\ i significantly alter the mechanisms of crater . 
formation which are associated with the crater 
| formed in snow by the relatively small HF charg2s? 


ai In areas where subsea permafrost is not present, ; 
the methods of EM-] are adequate for predicting crater 


dimensions, using the soil factor appropriate to the bottom 
composition. Cther arctic environmental factors are not 
expected to have a significant effect on underwater cratering. 


4 
« 
4 
: ! 
3.5.2 a Recommendations v 
Ss 


Cay To establish whether or not it is valid to scale : 
rx, 7 HE-charge crater data to the nuclear burst cases, computer My 
5 ‘ code calculations should be performed for the different charge f 
i - output characteristics and sizes. The results of the HE-charge 
‘ . d calculations can be checked against test data to verify the 
accuracy of the code(s). 


om oe 


Should a large-Charge Hi test (100-con TNT or more) be 


should be obtained. Early-time photography of the crater formation 

* should be obtained also. Additional crater dimension data from 
smali-charge HE explosions should be collected on a "test of 
Opportunity" basis, but a test series performed specifically to 
obtain crater data is not recomnended. 


Cot etal 


j 
t implemented in the Awcti®y, the Grat@r and ejecta measurements 
i 


The recommended series of calculations, experiments 
and field tests described ir. Section 2.7.2 should be planned 
such that the above questions will be answered. 


i i a On is 


i A research program to narrow the uncertainties in 
; underwater cratering should take a dual approach - to support 
7 | the collection of data to delineate the areas of bonded and 


unbonded permafrost, and to determine the appropriate soil 
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— factor to use for each type. Whether seismic io 
methods can be used to delineate the two types of permafrost f Hi 
areas or whether a program of subsea coring is required should : 
be investigated by experts in the field. The determination of 
soil correction factor should be accomplished by analytic means ' “9 
if possible, but laboratory or field testing might be required. 4 
This question should be studicu by experts in the area before H 
embarking on a program of field testing. 
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SECTION 4 
THERMAL RADIATION 


7 


About a third of the yield of a nuclear weapon deto- 
nated at Tow altitudes in the atmosphere is emitted as infrared, 
visible amd ultraviolet radiation with a pulse ‘width depending 
on the yield and altitude typically lasting for a few seconds 
for a megaton weapon. This thermal energy can be transmitted 
to large ranges in the atmospher2 and is usually readily 
absorbed in a tiiin Surface layer on most target material causing 
large surface temperature increases which can cause damage to the 


J 
Wo 
i 
x 


targets. The types of damage usually of concern result from 
fires that start due to ignition of combustible materials and 
from burns to personnel but also can involve thermo-mechanical 
loading due to very large fluences incident on hardened facil- 
ities such as radars. 
4.1 GP acctic Environmental Differences 

A large variability is expected in the effects of 
thermal radiation in an arctic environment because extreme 
variations in clouds, atmospheric moisture, visibility, precipi- 
tation, and the earth's surface occur more commonly in the arctic 
region than elsewhere in the world. These variations create con- 
ditions that can as much as double significant thermal effects 
or reduce them by even larger factors. The following briefly 
Gescribed effects will Le discussed further in subsequent 
sections of the handbook. 


4.1.1 GHB :sivsiity 


= S282 we see ele lee wlll le  —_— = 


4 

be na ats : . r ; 

y! TT? Surface visibility in arctic and subarctic climates dur- 
\e ing clear seasons is often exceptionally good because of the low 
5 

s humidity coincident with cold temperatuces and the absence of 

. 

J dust in the air. The northern and coastal regions during the 

5; warmer months are subjected to extensive sea fog and low cloud- 
’ 4-1 
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iness. Falling or blowing snow can reduce the visible range to 


2 wir 


U 

s 

! less than a mile during stormy periods. These extremes of atmos- 
( 

pheric comditions require examination of their effects on trans- 
mission of thermal energy. 


| 
‘ 4.1.2 QR roo ] : 
1 


4 Ld Arctic fog and precipitation generally reduce the range 
ermal phenomena. At the very low temperatures of the arctic 


‘ig winter the atmosphere is capable of holding very little mosture. = 
Such low temperatures as a rule are accompanied by minimal , 

“ surface wind, and these conditions together are favorable to ~ 

4 the formation of fog. Ice-fog crystals consist of many spherical - : 


particles and some hexagonal plates and columns of 2u to 304 dia-~ 

meter formed at about -40°C in high concentration that reduces -- 
visibility significantly. In addition ice fog can cause extinc- . 
tion of the infrered beam of an infrared guidance system. "4 


3 
¥ 
4 
* 4.1.3 OR arveco Surfaces ‘ 


Ground surfaces covered by snow and ice have a 
much higher albedo than bare ground in temperate climates. P, 


4 ver eeee : i. 
% ; The transmission of thermal radiation is considerably i | 
a 

‘ i enhanced by the presence of these high albedo surfaces. 

a 


Layers of cloud, smoke, or haze are other common albedo : 4 
surfaces, ! 


¥ 4.1.4 Cloud Cover 
WP c21019 cover 


\ fg The low dense cloud cover characteristic of arctic 
i areas can result in significant enhancoment of thermal environ- 
ment for targets at low altitudes from a low aititude burst. 
A high albedo ground surface is very likely in the arctic. The 


combination of a high albedo ground surface and a low cloud 


i 
i 
i 
i 
’ 
H 
i 
¢ 
cover results in a definite channeling of thermal energy and ee 
a marked increase in thermal fluences. 
i] 
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4.1.5 aa Humidity 


Even though t'e relative humidity is generally high 
in the arctic espezially over the ocean areas because of the 
low temperatures, the absolute conce cation of water vapor 
is much lower than in temperate areas. This results in less 
absorption of thermal radiation in the important infrared 
water vapor absorption bands and tends to increase the thermal 


transmission. 


4.1.6 PP Temperatures 


ed The low temperatures in themseives do not result in 


changes in the thermal transmission except for their influence 
in prodycing ice fogs, ice/snow surfaces, etc. 


Materials used for clothing and supplies in arctic 
climates do not possess the same vulnerability to thermal effects 
as materials used in less severe temperatures. Furthermore, 
cold temperatures reduce somewhat the vulnerability of most 
materials to thermal effects. 


—_ Besides a reduction in flarmability with reduction 


7 74 


in temperature, combustible material is less susceptible to 

thermal] damage when protected by snow and frost covering. 
Gharacteristic low humidity of the arctic air will somewhat 
mitigate the reduction of combustibility. In some of the 

tundra, expanses of coarse vegetation growing on a thick ' 
peaty layer might be subject to surface fires started by 

nuclear detonations. 


4.2 GBB ansmission Effects 


The quantity and effectiveness of thermal radiation 
that reaches a target is dependent on a large number of param- 
eters whose variability in an arctic environment is sufficiently 
great to produce a significant change in thermal radiation trans- 
mission. “fRie parameters to be considered here may be grouped 
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into two general categories; first, those parameters which 
determine the manner in which thermal radiation is scattered 
and absorbed and in the atmosphere are discussed here, and 
secondly, those that determine the manner and extent of its 
reflectiom will be considered in Section 4.3. 


Lio, | 


The thermal irradiance, H, received at a distance R> 
from a nuclear burst is given by the expression 


— Po Sa 


PT 
B= —z COS 2d, (4.1) 
4nR 
where P is the total power radiated by the burst as a function of 7 


time, cos A= ] if the receiver area is normal to the burst, and 
T is the transmission factor. T is used here in a very general iS 
sense and includes such effects as atmospheric attenuation, - 
surface albedo effects (ground, water, or clouds), and source 


asymmetries. 


The radiant exposure, Q, during time from zero to t is 
then defined to be 


t 


Q, ae a Gt, (4.2) 


° 


where in general P, T and R may depend upon time. If we assume 
that the transmission T and the radius R to a unit area facing 


the burst are constant, then 


Q > (4.3) 


where W is the yield of the weapon in calories and f is the 
thermal partition or efficiency, and fW = E is the thermal 
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yield. Exen though both T and R are seldom truly independent 

of time, sufficient accuracy may often be obtained by using 

the values corresponding to the time of peak radiated power. 
Often it is assumed that defining T = 1 will give the worst 

case environment. This assumption is usually good unless 

the conditions include source asymmetries and albedo from clouds 


and ground surfaces. 


In order to determine the amount of thermal energy 
actually transmitted to the receiver, allowance must be made 
For the attenuation of the radiation by the atmosphere. This 


“attenuation is mainly of two forms - absorption and scattering. 


There are no strong absorption processes for the visible wave- 
lengths but strong absorption bands exist in the ultraviolet 
and infrared. Scattering occurs with radiations of all wave- 
lengths. The state of the atmosphere in the visible region 
can be ~epresenteG ty what is known as daylight visibility. 


4.2.1 QP acctic visipizity 

7 There are several highly variable climatological 
characteristics that could significantly change the absorption 
and scattering of thermal radiation in the arctic and subartic 
atmosphere. As indicated although the arctic and subarctic 
regions are generally areas of comparatively low absolute 
humidity and little industrial dust, characterized by good 
visibility; the northern and coastal regions are e-*sonally 
subjected to extensive sea fog and low clouds, and falling or 
blowing snow can reduce visibility to less than one mile for 


extended periods in some areas. 


There is no single correct value of the atcenuation 
coefficient pw for any given set of atmospheric conditions. The 
value of uw is a function of both the nature and distribution of 
the scattering and absorbing particles, and also of the wave 
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= of the radiation involved. There is no simple average 


value of w because the spectral distribution of the radiation 


will change with th? 7istance "R" involved. In spite of the 
variable nature of uw, an assumption is often made that reasonable 
average values of uw can be determi:.ed in terms of visibility. 
This is not too unreasonable an assumption since that portion 

of the spectral distribution of radiated energy which penetrates 
any conziderable distance in the atmosphere is concentrated 
mostly iin the visible and near visible wave lengths. The con- 
ventional visibility as given in weather forecasts is generally 
the distance at which the transmission is reduced to 5.5%, i.e., 


pee Mes 1055. 


ey It is convenient to calculate the visibility on the 

asis of the above approximation to illustrate the dependence of 
atmospheric transmission and attenuation on visual properties of 
the atrosphere. There is wide discrepancy among vaiuves assumed 
for the distance called visibility and in relating that param- 
eter tc the optical properties of the atmosphere one should 
be aware of its very approximate and necessarily subjective 
nature. In most technical literature on atmospheric transmis- 
sion the term meteorological range is defined as that distance 
where the transmission is 2%, i.e., T= e = .02. Unfortunately 
some authors use the terms visibility and meteorological range 
interchangeably. 


The international code for correlating the condition of 
the atmosphere with visibility is given in Tahle 4-1. In temperate 
climates one may see variations in the visibility from the highest 
to lowest visibilities depending upon the concentration of aerosol 
particles from pollution sources. In most areas of the arctic the 
background level of pollutants at the surface is low leading 
to very high visibilities. However, the occurrence of‘certain 
weather conditions such as blowing snow, ice fog, etc. resuit in 
very low visibilities of one mile or less. 
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| ee 4-1 INTERNATIONA, VISIBILITY CODE 


Code Numter Description Visibility 


From 


Dense Fog ---- {55 yds) 
Thick Fog 50 m (220 yds) 
Moderate Fog 200 m (550 yds) 
Light Fog 500 m (9.6 mi) 
Thin Fog Y; kn 2 km (1.2 mi) 
Haze 2 kr (2.5 mi) 


| Light Haze 4 km (6 mi) 
Clear 10 km xm (12 mi) 
Very Clear , 20 km » (30 mi) 
Exceptionally 50 km 280 km (170 m:) 


Clear (Glasstone, 1977) | 


When one is conSidering the transmission for a bread 
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spectrum as results from @ nuclear weapon, the relation between 
the visibility and transmission is not as straightforward as 
indicated above. Scattering and buildup effects occur which 
result in a nen-exporential falloff of thle transmission. The 
various interaction cross sections vary as a function of wave- 
length so that integration of results across the broad wave- 
length must be considered. Extensive discussion of all aspects 
of this problem are presented in the DNA Thermal Sourcebook 
(Keith, 1973) and EM-1 (DNA, 1978) which is currently under 
revision. 


Fg Pigure 4-1 (Keith and Sachs 1977) shows predicted trans- 
mission as it varies with visibility of one to:30 miles. In the fig- 
ure the variation of transmission with the ground level visibility 
is noted as a function of ground range for a large yield weapon 
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Getonated at 5000 ft above the ground. The atmospheric profiles 
for every quantity except thr aerosol concentration are unchanged 
for the various calculations. The aerosol concentration at : 
ground level is adjusted to give the desired visibility; then a 
constant exponential lapse rate is defined between ground level 
and 5 km. The aerosol concentration profile, as in most atmos- 
pheric models, is assumed to be unchanged above 5 km altitude 
regardless of the visibility at the surface. This implies that 
essentially none of the particulate contaminants are carried 
above 5 km altitude. 


= The calculated transmission increases as expected with 
increasing visibility. These curves represent a burst at 5000 ft 
with no albedo surfaces present. The representative arctic humid- 
ity is 1 a/n?, and the spectrum corresponds to a high-yield burst. 
At a range of 10 km the transmission for 30-mile visibility is 
about 50% higher than that for 6 miles and 1¢C times that for 
visibility of 1 mile. At a range of 30 km the transmission for 
30-mile visibility is about 3.5 times that for 6-mile visibility, 
ana transwission for l-mile visibility is practically negligible. 


The transmission curves given above refer to a ground 
level absolute humidity of 1 g/m. As shown in Table 1-5 the 
absolute humidities of the standard atmospheres are .46 in 
January and 5.6 in July. The effect of the humidity for a 
particular visibility is very small in the visible region of 
the spectrum but can strongly affect the portion of the weapon 
energy emitted in the infrared. The clear visibility curves 
from reported results correspond to a water vapor concentration 
that is considerably in excess of that normal for an arctic 


winter. 
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ee In Figure 4-2 the effect of changing the absolute ; 
umidity at ground level is shown. All quantities except the : 
water vapor concentration are held constant. The ground level 
concentration is set egual to the desired value and an exponen- 

tial lapse rate is defined between ground level and 5 km. The 


profile above 5 km is assumed to remain the same for the various | 
humidities. As expected the effect is not as large as noted 
a 


with changing the visibility. A surface burst where the entire 
path is along the ground should show the maximum effect. The 
overall slope of these curves is determined by the visibility, 
which here is 30 miles, and the relative placement shows the ] 
effect of the water vapor. For the higher concentrations the 
water vapor absorbs strongly within the first km of the path ] 
until the infrared energy is depleted; then the transmission 
versus range is determined by the visibility. Note that increas- 
ing the concentration beyond 5 g/m? has a relatively small effect. 


Tre large yield bursts have a relatively larger fraction 
of the yield in the infrared region where humidity effects are 
important, therefore, these results represent a reascnable upper 
limit to the effect of humidity on the transmission. The dif- 
ferences shown are not of importance considering the uncertainties 
in the other meteorological parameters. 


4.2.2 OP i<- Fos 


Some of the effects of ice fog on infrared transmission 
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in Alaska have been reported (Kumai and Russel, 1969). Besides 
reducing visibility significantly, ize fog can cause attenuation 
of the infrared beam in an infrared guidance system. The optical 
properties of fog depend on the number concentration and size 
distribution of the particles, which can vary significanily 
during different meteorological conditions. ice-fog crystals 
appear as initial stages in the formation of snow crystals at 
about -40°F. At the very low temperatures often occurring during 
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and surface wind is almost invariably calm. Since conditions 
are frequently conducive te formation of ice fog from any source 


A of water vapor, the frequency of ice foq has increased with human 3 
p activity in the arctic. 
f Table 4~2 from the report summarizes the important M 
" physical properties of the ice fog at -39°C and -41°C. The ice- - 
\ fog distribution at -39°C is shown in Figure 4-3. The number of Fs, 
' particles and the mass of fog per unit volume are shown as a “1 
r function of particle diameter. * 
P it. 
iat 
TABLE 4-2 7M 
PHYSICAL PROPERTIES OF ICE FOG AT FAIRBANKS, ALASKA , 
: Air “a 
.) Observations N ; ¥ ode adh aan Ar temp vines ° ; . 
; (no.feom")  (#) (wu) te) Ce) UEC) g/m) i 

No. 1(Fig 4-3) 140 3.0 1.5 12.0 0.5 <-39 0.08 

; No. 2 90 1.5 1.5 12.0 0.5 -41 0.02 


N = total concentration UNCLASSIFIED 


r = mode radius = radius corresponding to the maximum number 


mode 
. of ice-fo¢ crystals 


' r_.. = minimum radius 


t) r = maximum radius 


3 
o 
* 


Ar = radius interval containing n(r) crystals, where N =D alr) 
x 


. iF L.W.C. = liquid water content 
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BBD Figure 4-3. Size and mass distribution of ice-fog crystals 
formed at ~-39°C ambient temperature. 


Computer calculations of the attenuation and back~ 

a seattering of radiation by ice fog alone showed them to be 
within the same order of magnitude as those for water fcg of 
equivalent foa0 concentrations and observed wavelengths. The 
opticai constants used in the calculations were considered to be 
known much less exactly for ice than for water. Calculations 
made for the distribution in Figure 4-3 were presented and are 
reproduced as Table 4-3 (Kumai and Russell, 1969). 


omen e 


These calculations were done for narrow wavelength 
i intervals in the infrared and do show detailed differenves at 
specific wavelengths due to the different scattering character- 
istics of the ice crystals end water droplets with equivalent 
. amounts of water involved. The same size distributions were 
assumed for these calculations which may not be a realistic 


a et 


assumption since the size distribution of the ice fog is con- 
Siderably different from other types of fog. 
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TABLE 4-3 
INFRARED ATTENUATION COEFFICIENTS (m2) 


Water Wavelength, Microns 
Content 2.7 4.5 5.75 9.7 


g/t? 


201153 | .01212 { .0132) | .01304 | .00513 | .00679 
.02306 | .02424 | .02642 | .02608 | .01026 | .01359 
04611 | .04848 | .05283 | .05216 | .02052 | .02717 
-06917 | .07272 | .07925 | .07824 | .03078 | .04076 - 7 


} 


Water-Fog - 

201153 | .01217 | .01314 | .01282 | .00784 

, 202305 | .02434 | .02628 | .02564 | .01568 
; -04611 | .04868 | .05256 | .05127 | .03137 
: -06916 | .07303 | .07885 | .07691 | .04705 


Another paper on visual range in polar regions 

(Mitchell, 1958) also states that the visual range in ice 

fog is characteristically very low, frequently less than a 

quarter of a mile. The total particulate water content of 

ice fog is comparable to that of other fogs, as shown in 

y Figure 4-4, put the average ice-fog particle is smaller. Thus 
the ice-fog contains more particles and favors greater optical - 
scattering. Because ice fog is often a man-made phenomenon, 

J it is e@ particular problem at arctic military installations. 


The details of the transmission of thermal energy : 
in the atmosphere depend upor the aerosol particle distribu- 
. tion as well as the concentrations. The scattering phase 


functions determine the details of the thermal transport includ- a 
) ing the angular distributions and transmission factors. The [ 
4-14 2 


~~ ‘ y Breet. ACES ORS Panes tal Rot noe IB 
m te . , 
e fence ba Se ee eee deo cea tnd elie ate dni. Ne cal tee ce ata aS 


~| 


=f 


His PAGE 1S BEST QUALITY PAu a CABLES 
Pru COPY FERAISID TO EEA nae 


= Smee ww 


t 

! 

t 

| 

| 

| 

at. é . S 

os ° 

. be ee, ete ae Se. et re 
~ r Taste, 645 

N if a 


s 
ea 
we ee eee 


PE, ATtvE eM veut 
t-J e 
oe e 
T 


fc) 
~“ 


o 
oa 


ema Pen pet ae caren ts So seein ao ale Sie aS ee 


Set atnene aieeakel ine heeeenaeoee 


See 
od 


aes 
oe AIRED on 


6 © 2 30 «0 SC 60 70 80 90 100 
PARTICLE DIAMETE® MICRONS 


= Figure 4-4. Comparison of typical particle volume distri- 
butions for ice fogs and water fogs (Mitchell, 
1956). 


ae functions depend upon the particle properties and size 
distribution. No calculations of the Mie scattering functions 
have been made for distributions peculiar to the arctic region. 
All transmission curves in this section were prepared by using 
available phase functions derived to represent the aerosol 
distributions found in temperate polluted areas. 
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4.3 GP arreco Surface Effects 


The nature of the clouds and the condition of the 
surface in the vicinity of a nuclear detonation can considerabiy 
alter the amount and direction of thermal radiation reachiny a 
target. Although there is great geographical and seasonal 
variability throughout the world in these parameters, more 
extreme variations occur in arctic and subarctic areas. Th¢ 
number of variables involved and the range of possible varia-~ 
tions in these parameters make a comprehensive consideration 
of the reflection or albedo problems very comple». 


|g In terms of thermal radiation, surfaces or atmospheric 
anomalies which reflect radiation are known as albedo surfaces. 
Surface albedos range from 0 to 1, where the value 1 indicates 
a perfect reflector. Typical albedo surfaces are the ground 
plane, especially when covered by snow, ice or water, a cloud 
layer, and a smoke or haze layer. Dense clouds may have albedos 
as high as 0.9. 


Even if albeco surfaces are not present, ail sides of 
an object will receive radiation even though the side facing the 
firebzl1 will usually receive the dominant exposure. A portion 


of the radiation traveling upward is lost to space with relatively 


littie being scattered downward unless clouds are present above 
a berst. If a cloud layer is present a large portion of the 
incident radiation will be diffusely reflected from the surface 
with a small fraction being diffusely transmitted. A typical 
ground surface also reflects the radiation diffusely with the 
albedo varying from near zero to near unity for snow or ice 
surfaces. Some materials such as water or ice may also have a 
fairly large specular reflection. The thermal exposure for 
«.rgets bounded by clouds and a high albedo ground surface can 
be several times the vacuum exposure value. 


i+ bt omet ome GES ARS 


mA Se Nn 


we) he he ie as [7 


all” of PALS poet 


per 


= tn iia eee 
>, B 


- ft yay- 


LOS OL FL 


To 
4-4, 4s 


EEE EEE LEE AE NEE CE OEE Re ee 


Searonal aspects of the albeco of arctic surfaces north 
of 65° N latitude h-ve been examined (carssen and Orvig, 1962) 
from data in the literature. At high latitudez seasonal varia- 
tions in albeto are largely determined by the presence or absence 


Mt Pe: 


of snow cover. In the tundra zone the contrast as greatest. In 
forested tundra albedo from snow cover beneath the trees is 
approximately twice as great as from the ground vegetation 
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except in a chose-crown forest where snow is caught on the surface 
for a relatively short time. Different tynes of ice reflect 
differently, and although open water and low cloud cover usually 
are found coircidentally, little information is available cn the 
albedo of opea water containing ice. Albedo stereograms in 

Figure 1!-!6 represent the seasonal and latitudinal changes in 


arctic surface albedo. 


{ 

| 4.3.1 a Experimental Results 
] a A series of arctic transmission measurements were 
made by the United States Army Electronics Command, Fort Monmouth, 

New Jersey (Cantor and Petriw, 1964) in Greenland under various 
| weather conditions with particular emphasis given to the albedo 
effects of the ground surface and cloud layers. 


a Under many atmospheric and surface conditions, the 
indirect radiation effects can equal or exceed the direct radia-~ 
tion. A 2w detector is used to simulate a flat plate receiver 
so that the full indirect as well as the direct transmission 


can be measured, 


Sy Such a receiver, near ground level, was employed with 
a 6500°K point light source about 400 feet above the surface, 
under generally hazy atmospheres on the New Jersey shore from 
October 1960 to February 1961. These tests indicated sharp 
increases of radiation under relatively high surface albedos. 
The maximum surface albedos, however, are readily obtainable 

in the arctic or antarctic regions. This led to studying trans- 
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mission effects of an energy source located between two high 
albedo surfaces, a snow-covered surface and an extensive cloud 


cover at Camp Century, Greenland. 


Cantor and Petriw give a very complete description of 
the experiments of thermal transmission made at Camp Century, 
Greenland in March, October, and November 1962, where measure- 
ments were made on the Greenland icecap to try to maximize the 
effect of albedo. The terrain was essentially flat and had an 
albedo of alzost unity extending for about 100 miles in all 


directions. 


The light source was a xenon flash tube surrounded 
by a 10"~diameter opal-covered sphere about 100 feet above the 
surface. The blackbody temperature of the source was about 
6500°K. Photomultiplier tubes were used as 27 detectors with 
no variation in the field of view possible. Occulters were 
used to block the direct irradiance. Measurements were made at 
ranges of 0.13, 0.5, 1., 4.5, 7.6, and 10.3 miles whenever 
conditions permitted. 


er conditions, experimental configuration, and data obtained each 
night measurements were made. Plots of the total and scattered 
transmittance for each night and many of the signal variations as 
@& function of time are also given. Very large, short-term varia- 
tions of up to 160% were noted in the intensity in time intervals 
of less than 30 seconds. These largest variations occurred in 
periods of high visibility and steep temperature inversions with 
smaller variations occurring for smeller temperature inversions. 


fg Figure 4-5 is derived from the summary graph from the 
report and gives the transmission as a function of range. The 
labels refer to the visibility in miles followed by the altitude 
’ of the cloud layer in feet. Several interesting effects can be 


4. 


Their report gives a detailed explanation of the weath- 
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seen a comparing the transmissions for the different atmos- 


pheric conditions. First note that the case referred to as good 


to excellent visibility with no clouds shows essentially no trans- 


mission variation with range over the fairly short range of the 
experiment. The addition of a cloud layer at 2000 feet did 
introdace an increase due to the ducting effect but the change 
is of the order of only 25%. A much larger increase is noted for 
the case with 2000 ft cloud layers and a lesser visibility of 

16 miles. This general t,pe of behavior is noted in Monte Carlo 
calculations of this effect as will be discussed in the next 
subsection, but the calculated effects are much less extreme. 
Note that for a visibil’ ty of 1 1/2 miles and clouds at 2000 ft 
there ts an initial increase in the transmission above unity 
followed by a decrease in the transmission. 


| dl For the 1/2 mile visibility case with no clouds a 
rapidly decreasing transmission is shown. For a cloud layer 


at 300 feet anc 1/2 mile visibility the transmission increases 
to over 5 at a range of about 1 km beyond which no measurements 
were made. One would then expect a precipitious drop with 
increasing range. 


a The dashed curves refer to cases in which the visi- 
bility reduction was due to ice fog rather than water fog. 
There is an indication that for the same visibility the trans- 
mission may be higher for ice fog indicating larger scattering 
contributions and differences in the scattering phase functions. 


The experimental uncertainty is fairly sarge; so 
quantitative measures of the effects should not be derived 
from these experiments. The measurements do not extend to 
long ranges as are necessary for nuclear weapon thermal 
prediction methods. The general trend of the results does 
agree with results of Monte Carlo calculations of thermal 
transmission including the effects of albedo surfaces. 
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4.3.2 a Monte Carlo Calculations 


d oat This section will summarize work at RRA (Wells, 
| ollins, Marshall, 1969; Wells, Collins, Cunningham, 1966) 
begun im the mid 1960's and the further work at KSC (Keith, 1973) 
in developing Monte Carlo codes describing the thermal radiation 
j transport model atmospheres. Given a complete specification of 
the atmospheric parameters a calculation of the transmission 
] with these codes will probably be the most accurate that can 
be obtained theoretically. 


Monte Carlo calculations (Wells, Collins, Marshall, 
1969) have been made of the transmission for a 6000°K blackbody 
source at 1 km altitude for the model atmospheres representing 
meteorological conditions for an arctic case and three mid- 
latitude cases - summer, winter, and a winter inversion. 
Calculations were run with and without a ground albedo factor. 
: Table 4-4 lists some parameters used in the four atmospheres 
. compared on the graph. The only difference in the summer and 
winter midlatitude case is in the absolute humidity. The winter 
case with the inversion added the very low visibility region of 
verosols below 2 km. The arctic case was chosen to provide 
exceptionally high visibility and zero humidity, which as nected 
previously will increase the transmission by a relatively small 
: factor. A snow-covered surface with cloud distribution was also 


. assumed. 


\ 
] 
] 
‘| 
] 
I 


In Figure 4-6 the results of their calculations for 
a target on the ground surface are summarized and replotted. 
The differences between the summer and winter midlatitude clear 
visibility cases are relatively small as expected. The excev- 
tionally clear visibility assumed for the arctic case gives a 
much larger transmission factor. The winter inversion case for 
a haze visibility of 2.2 km or 1.4 miles results in significant 
attenuation even at fairly small ranges. The effect of changing 
‘the ground albedo from 0 to .9 is seen to result in a significant 


4-21 


Ne eg eee ee OE 


Re eee EE Oe TS EES ERY OT SEE TES NTE ET RET AROS OT J TIRE rn ge mem ramen sve sams 


er 


= ~ 
~~. Frew 
Winter, 
Summer 
Mid-Latitude 
5 \ Winter 
Gg Inversion 
“10 
e 
wv 
c 
oO 
E — Cl 
N 
\ 
\ 
\ 
N\ 
‘X 
» 
\ 
=e 
-3 *S « 
10 Sau 
Ground Albedo = 0 
— —~ Ground Albedo = +? 
-4 
10 
30 60 90 


Horizonto! Range (km) 


Gey icve 4-6. TRANSMISSION IN MODEL ATMOSPHERES 
4-22 (from Wells, 1969) 
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TABLE 4-4 


METEOROLOGICAL CONDITIONS 
USED FOR MODEL ATMOSPHERES 


« 
4 
: 
: 
§ 
; 
| 
i 
bs 
> 


Ground Level Ground Level 

Atmospheric Absolute Visibility Aeroso) Size 
Model Humidit m3 kn Distribution 
Summer Midlati tude 12 18 N(r)er* 
Winter Midlatitude 3 18 N(rjer”- 
Winter Midlatitude 

(Inversion Profile 3 2.2 N(rjer”> to 2 km 

to 2 km altitude) ., altitude 

N(r)er ~ above 2 kn 
altitude 
Arctic 0 148 N(r)er"" 
(Exceptional ly 
clear) 


* ry - radius of aerosol particles 


N(r) = aerosol concentration as a function of r 


buildup of the transmission for the midlatitude cases considered. 
The effect is only about 5% for the arctic case so that the 
effect of one high albedo surface for exceptionally high visi- 
bilities is not large. For targets above the surface a larger 
effect is noted, 


cmd In a recent study (Kaman Sciences, 1978) calculations 

with Monte Carlo computer codes were done including the effect of 
introducing albedo surfaces. In Figure 4-7 the effects of various 
combinations of albedo surfaces are compared. The burst conditions 
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| Figure 4-7. VARIATION OF TRANSMISSION WITH ALBEDO SURFACES 
AT 12-MILE VISIBILITY. 
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~ oe Same as in the transmission figures in Section 4.2 for 


: varying visibilities and humidities. The ground level visi- 

. bility is i2 miles and the absolute humidity is 10 g/m’. The 

addition of a single albedo surface, either snow covered ground 

or a cloud layer, causes a modest increase in the transmission 

of about 20% to 30%. The transmission with only snow cover is 

seen to be somewhat higher than with only cloud cover. The 

effect of a clouc layer at 8 kft becomes less with increasing 

range, and for ranges greater than 30 km a reduction in trans- 

mission will occur. The effect of a 16 kft cloud layer is 

; Similar except the increase and decrease occur over a longer 

{ range. Tne combined effect of snow cover and a cloud layer can 
be quite large as noted by the two upper curves on the figure. 
An increase of up to a factor of about 2.5 is possible. The 
effect is larger at the small ranges for the lower cloud layer. 
At much larger ranges than shown on the figure a reduction in 
transmission will result. In Figure 4-8 similar curves are 
given for a visibility of 30 miles. The same general trends 
relative to the curve with no aloedo surfaces are obtained with 
the albedo effect being somewhat larger with the higher visi- 
bility. The presence of the albedo surfaces causes a large 
effect in the transmission factors but since the effect depends 
upon the cloud height, it must be quantified for specific cases 


of interest. 


ee The effect of two albedo surfaces on the transmissior 
is seen to be very large. The calculations of RRA were done 
! with too large a visibility to represent realistic Arctic 
: conditions and the KSC calculations were done with a visibility 


lewer than can be expected in the Arctic. Neither set of calcu- 
{ lations was done with aerosol scattering functions representing 
aerosol concentrations appropriate for arctic conditions. 


wee 


7 5 5 . CH 
: : RRR gat gad 5 4p", da®. O38, KU a 
PRT eu wa Lh Re ee) ew Ce ov ta* at A ee ee he 
ver cece re Se re emmmmume orem 
exer OTE, ERIE SE: wencTer -arernne-vegs commen reer nes me rm rewwmnruenn or om ern men — 
’ 


se 


ww FIGURE 4-8. 


* 


¢ Lg 
MITTAL Ne Kn Mai Man Yat 


Transmission 


Cloud of 8 kft 


5000-ft burst 
(Keit', 1978) 


0 10 
Ground Ronge (km) 


30-MILE VISIBILITY 


4-26 


20 


No Snow or 


Cloud ot 8 kfi with Snow 


Y Cloud at 16k Ft with Srow 


Cloud at 16 kft 


Cloud 


ee ence te aa ee 


30 


VARIATION OF TRANSMISSION WITH ALBEDO SURFACES AT 


oS 


eH fet eel eee 


-_—- 


SOR ORL TAS OT OK WY 


Oca TTT COLE LS SP ETI 


a 
Ly 


4.4 | | Example Thermal Exposures _ 

Figure 4-9 gives curves of radiant exposure from 1 to 
1000 cal/em* versus ground range from the point beneath a high-~ 
yield burst at 5000 feet altitude. For a very clear atmosphere 
(30 mi visibility) the exposure is 10 calories per square centi- 
meter at a range of about 15 miles. This same damaging exposure 
is seen to occur approximately half as far from the burst ina 
thin fog with visibility of a mile. The upper curve shows the 
exposure for a visibility of 30 miles with snow cover and a 
Cloud layer cottom at 8 kft altitude. The 10 cal/em? exposure 
ececurs at a ground range of about 22 km which is about 50% 
larger than for the 30 mile visibility case with no albedo sur- 
faces and a factor of three larger than for the 1 mile visi- 
bility case. 


A recent study (Keith, 1979) considered the thermal 
environment for a selection of Soviet cities considering the 
wide variation of meteorological conditions that may result in 
this area. Representative and extreme days were defined which 
happen to be of some interest for considering the thermal 
environment in the arciic. The results are shown in Figure 4-10. 
The extreme low day refers to a .3 mile visibility with no cloud 
cover and medium ground albedo which corresponds roughly to a 
heavy ice fog with no cloud. The extreme high day refers to a 
50 mile visibility with a 16 kft complete cloud cover and a 
medium ground albedo which corresponds roughly to a clear winter 
arctic day with no fog or haze. For these extreme cases the 
range corresponding to 10 cal/em@ varies from about 25 km to 

! 4 km which is a somewhat larger spread than noted in Figure 4-9 
and is much larger than typical for temperate climates. 
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—_ Al) of the examples given previously refer t> the 

thermal exposure received by a target on the ground surface. 
An important consideration is the thermal exposure on air- 
eraft in the vicinity of the nuclear burst both for consider- 
ation of a safe escape range for a delivery aircraft as well 
as determining the thermal damage potential for attacks on 
Soviet airbases in the Arctic. No calculations are available 
for the specific Arctic cases of interest but the transmission 
results presented in Figure 4-11 shows the type of effect that 
will be experienced. The transmission is given as a function 
of horizontal range and altitude from a burst at an altitude 
of 1 km. The solid lines give the exposure contours for a 
case with a cloud bottom at an altitude of 6 km and with a 
medium ground albedo corresponding to desert sand. The trans~- 
mission for the Arctic albedo case would be somewhat larger. 
The dashed lines refer to the same burst conditions without 
albedo surfaces. 


= The ratio of the solid to dashed contour values at 
the same point in space give an indication of the buildup 
introduced by the albedo surfaces. Note that the dashed con- 
tours are all less than unity and tend to decrease with 
increasing range and decreasing altitude in a regular pattern. 
The introduction of the albedo surfaces leads to a much more 
complicated spatial Cependence for the transmission because of 
the complicated interaction occurring between the attenuation 
and scattering properties of the atmosphere and the diffuse 
scattering of the two albedo surfaces. 


rf Vulnerability levels of about 100 cal/em@ are realistic 
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for aircraft which in previous figures would occur at a range 
of about 7 km near the ground for the albedo case. xatios of. 1.8 
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are noted near the ground for these ranges. At higher altitudes 
ratios of I.5 are obtained. These represent large differences 
in exposure levels and could represent the difference between 
sure~Safe end sure-kill environments. 


a Birectly above the burst ratios of 2.5 are obtained. 
, t larger ranges representing lower exposures ratios near 3 
are experienced at the lower altitudes. In general thermal 
effects became more important as the yield is increased and as 


shown in the last three figures, the effects of the albedo 
surfaces are largest at the longer ranges. 


4.5 we fhermal Effects of Underwater Bursts 


Thermal effects of underwater nuclear detonations 
are generaliy ignored. In Chapter 3 of DNA EM-1, only land 
surface and subsurface bursts are treated in any detail. It is 
stated that in tne case of underwater bursts, thermal effects in 
the atmosphere are usually insignificant, and the tact that a 
20 kt burst in 90 feet of water produced negligible thermal 
radiation is cited. (20 kt at 90 ft is a very shallow detona- 
tion -- ~3340243, ) The presence of ice in the Arctic would tend 


to reduce thermal effects in the atmosphere even more. 


Since the thermal energy of an underwater detonation 
is largely ehsorbed by the water, the question arises as to 
whether there will be left a body of heated water sufficient to 


era ine 


ereate and maintain an ice~free pool in a region of otherwise | 

total ice cover. Neither DNA EM-1 nor the Underwater Handbook 

addresses this question. A limited amount of experimental data 

has been collected on the temperature changes produced in water 

by underwater explosions. While most of the Cata have been a, 

acguired on experiments conducted with a steam-generating explo- 
a 
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sive (Lithanol), developed for the purpose of simulating the i 
bubbie behavior of underwater nuclear detonations, a few parallel 
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—_ eaducted with Pentolite, a conventional high explo~ 
sive, and water temperature data were collected during Operation 
Wigwam. The igh explosive tests were conducted in Chesapeake 
Bay, Puerto Rico, and Panama City during the period 1965 to 1969. 

Lithanol charges up to 13,000 pounds were used. 


7 None of these data ha: demonstrated any significant 
eating of the water. The results of the non-nuclear tests have 
been publishet in the open literature (Young, 1971 and 1973). 
Young and Scott, 1970, summarized the existing experimental and 

theoretical mowledge of the heating of water by underwater 


explosicns and examined phenomena that hac not been treated in 
earlier studies (e.g., Young, 1968 and Young and Scott, 1968). 
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fg A simple calculation will show that it is not sur- 
prising that significant heating effects of underwater explco- 
sions have not been observed. Assume an explosion deep enough 
that the first bubble at its maximum radius does not penetrate 
the surface, Sav d = 24owl/4 , Assume that 100% of the available 
energy remains in the bubbie and is used to heat the water in 
the cylinder less half sphere that is between the bubble and the 
surface. It can be shown that the temperature rise in that 
volume of water is about 2.4wi/4 °C, where W is in kt. A 10 kt 
detonation would heat this water less than 5°C and a 100 kt 


explosion less than 8°C. 


The conclusions of Young and Scott, 1970 provide the 
best summary of what has been found in the investigations of i 
the heating effects of underwater explosions: 
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4.6 GP sem Damage Effects 


= Thermal damage effects result from absorption of the 
ermal energy on the target accompanied by a surface temperature 
increase during the delivery time of the nuclear weapon pulse. 
The temperature reached in the target depends upon the thermal 
characteristics of the material, the thermal pulse amplitude 

and duration, the thickness of the material and the absorptivity 
of the surface. The range of magnitude of thermal effects ranges 
from personnel burns that can occur at levels as low as 2 cal/em? 
to massive melting and ablation of metals in blast hardened 
structures requiring a thousand cal/em? or more. 


= Thermal burns on personnel in the Arctic will be reduced 

cause of the amount of exposed skin will be much less than in 
temperate climates. The degree of incapacitation depends upon 
the fraction of the body burned as well as on the severity of 
the burn. First degree burns can result from exposures as low 
as 2 cal/om? for low yield weapons, but a first degree burn must 
occur over most of the body to produce a casualty. Thus, one 
would not expect casualties from such burns in the arctic. A 
less extensive second degree burn may cause a casualty but this 
takes about twice the exposure for a first degree burn. A miti- 
gating factor is that the parts of the body most likely to be 
exposed are the face and hands, and a burn on these portions of 
the body affects performance more than other parts. 
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Flash blindness casualties may be affected in the artic. 
During the day, the casualties will probably be less than in 
temperate climates because the high level of light will generally 
require some eye protection which will reduce the thermal radia- 
tion received by the eye, even under the possibly higher arctic 


* transmission. Even if no eye protection is worn, the pupil of 


the eye will be very small and will reduce the problem. Dur- 
ing the long winter nights, however, the eye will be more 
sensitive to flash blindness. 


The type of clothing worn in the arctic will tend 
to reduce the effects of burns. In temperate climates only 
thin layers of material are usually worn, and if these ignite 
or melt as can happen with synthetic fabrics, an extensive burn 
can result. Dark materials will be affected at levels as low 
as 10 cal/cm* 
require twice as much exposure because of their smaller absorp- 


for low yield weapons while white materials may 


tion. In the Arctic since one would probably be wearing several 
thickness of material, a burn would be much less Likely to 
reach the skin. 


ae The effect of the frost covering surfaces in the arctic 
w 


111 be to reduce the absorption of the thermal energy because of 
its high reflectance. However, the thickness of the frost would 
be crucial because the fraction of the energy absorbed during 
the early part of the pulse could melt the frost and expose the 
underlying surface to the later portion of the pulse. No dis- 
cussion of this effect was noted in the literature, but the 
magnitude of the effect can be estimated as follows. 


od Consider the surface loading expressed as the equivalent 
g/em” of water on the surface. Then with the assumed standard 


temperature of -24°C it will take about 24 cal/em? to bring the 
frost to melt and another 861 cal/em? to completely melt it for 
each g/cm” of water Inading, for a total of 105 cal/om” of absorbed 
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energy. The albedo can be as large as .9 for fresh snow and 
frost. Therefore, about 1000 cal/om* of incident energy is 
required to completely remove a frost equivalent to 1 g/em* of 
water. ‘The hoarfrost and snow buildup can be easily this large 
under arctic conditions; so that this would result in a very 
effective thermal shield even at high exposure levels. Note 
that a frost equivalent to a .01 cm rain would require about a 
10 cal/on* exposure to disperse and expose the underlying 
material. This could be of importance in considering the 
damage ta combustible inaterials and other materials whose 
damage level is low such as canvas tents and truck tops. 


Fg The ignition threshold for materials such as leaves, 


and newspapers is typically defined for conditions represent- 
ing a nominal humidity of 30% ~ 40% in temperate climates. The 
relative humidity in the Arctic is generally much higher than 
this, bot the absolute humidity or the moisture available in 
the air is low in the Arctic temperatures. This might result 
in an effective lower humidity for these mate ials and a lower- 
ing of their ignition threshold. 


> During the warm weather months, one would expect the 
ignition of fires and their subsequent spread in inhabited areas 
or forest and dry tundra areas to progress much as in temperate 


areas, During the cold winters, however, the lower temperatures 
imply a larger heat input to raise materials to the ignition 
temperature and suStain burning. Thus, one would expect exten- 
sive fires to be a less significant damage mechanism than in 
temperate Climates. Because of the extre’2 conditions existing 
in tt@ arctic, however, loss of shelter becomes a-very signif- 
icant factor in survivability and in retaining the capability 
of performing a mission, : 
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The thernal damage threshold for hardened targets 
such as radars which have been designed to survive large blast 
overpressures would probably not change under arctic conditions. 
The amount of energy necessary to raise the temperature froin the 
low arctic value would be negligible compared to that required 
to cause the relevent material damage. It is possible that for 


specific desians the high surface temperatures in conjunction 
with the lower overall structure temperature might result in 
larger thermo-mechanical loading and an increased warping force. 
If this occurred in materials that became more brittle at the 
lower temperatures, then there might be some chance for the 
structure to suffer damage at a lower exposure due to the lower 
temperatures. No Studies of this applied to weapons effects 
were found. 


4.7 a Conclusions and Recommendations 


ay The arctic is characterized in general as a region of 

elatively large surface visibilities with the high probability 
of. high albedo surfaces in the form of snow or ice covered terrain 
and low cloud iayers. This combination leads to a very high 
transmissio1 of thermal radiation as compared to average condi- 
tions in tenperate locations. At the same time erctic meteor- 
@logica? conditions result in the large probability of ocecur- 
rence of water and ice fogs and blowing snow which tend to reduce 
the visibility to less than 1 mile when these conditions are 
existing. This is a much smaller visibility than will be found 
in heavily polluted temperate climates. Thus, an extreme varia- 
bility in possible thermal damage ranges must be expected in the 
arctic depending upon the specific meteorological conditions at 
the burst point. 


The low absolute humidity characteristic of the arctic 
does result in an increase in the transmission of infrared energy 
as compared with temperate cases with the same visibility. Be- 
gause of the relatively small amount of infrared energy in the { 
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output spectrum of nuclear weapons the increase in overall 
transmission is small. The uncertainties in the other meteor- 
ological parameters are larger than this effect. 
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The experimental data showing the transmission of 
thermal radiation under various arctic conditions are very 
limited amd can not be used independently as a prediction 

method for general transmission calculations in the arctic. 


The data do indicate a very pronounced ducting effect for cases 
in which a high albedo snow layer exists in conjunction with a 

cloud layer. Enhancements as large as a factor of 100 over the 
direct exposure were noted for low visibility cases. However, 

because the direct exposure may be low for these cases, the 


position for a high visibility. Enhancements of a factor of 
two were noted in the high visibility cases. 


Two different Monte Carlo calculations have heen made 
of thermal transmission including the presence of two reflecting 
surfaces as well as attempts at handling the problem analytically. 
The simple techniques result in significant over-estimates of the 
enhancement due to poor handling of the attenuation in the atmos- 
phere even under high visibility conductions. The Monte Carlo 
calculations indicate possible enhancements of as much as 2.5 
due to the ducting effect over the general region of interest 
for thermal damage. 


Reliable Monte Carlo calculations require careful 
specification of atmosoheric parameters includino detailed infor- 
mation on the suspenaed particulate matter. The current standard 
atmosphere profile tables include a 75° latitude model which is 
suggested for use at ail higher latitudes. For the purpose of 
low altitude nuclear effects this is not a serious problem since 
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] total exposure may be less than would be noted at the same 
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the available meteorological information gives a good descrip- 


tion of the arctic atmospheric conditions near the surface 
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which is of mest importance. 


& few studies have been made of the particle distri- 


Yaa 


butions and concentrations for various conditions in the arctic. 


These have usually been oriented towards providing information 


on the visibility for aircraft operation. Recent measurements 


my 


Ss 


have indicated that layers..of aerosols at altitudes of 1 km to "i 

3 km exist often during the winter portion o. the year leading "e 
to lower visibilities in these layers than observed on the i 
surface. This could lead to predictions of a higher thermal Pra! 
fluence than actuzily exists for these conditions if the current “A 
techniques of using the surface visibility to specify the aerosol » 


distributions are used. One calculation has compared the rela- 


tive transmission of ice fog and super-cooled water fog. There , 
has been no general study made of the transmission for various i 


wavelenaths for the observed particle size distributions and 
the effects expected in the transmission of rvciaar weapon 4 


thermal energy. 


No studies were found relating to the change of thermal 
damage thresholds due to arctic conditions. The discussion in 


Section 4.6 reveals that one would expect the thermal damage ae 


ww 
threshold to be larger under arctic conditions except in par- I 
ticular cases where brittle materials conceivably would be ™ 


‘ 
13 
subjected to a more damaging thermo-mechanically induced force me 


than under temperate conditons. , 
: ie 

4.7.2 JB Reconnendations : 
Se 
Analyses of the meteorological information available ‘te 
arctic sites should be made to define the relative occur- ot) 


rence of high and low visibility conditions on a seasonal basis. 
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Cloud layer definition, aerosol concentrations and distributions 
and ice and water fog concentrations and distributions. 


The parameters of interest include visibility, ground albedo, 
From these analyses reasonable model atmospheres 

could be defined for purposes of thermal transmission studies. 

Calculations should be made of the thermal transmission for 

these model atmospheres emphasizing those cases which result 

in an enhancement over that expected in temperate climates 

such as those involving the ducting effect in high visibility 


conditions. 


Fg From these calculations, figures showing the predicted 
exposures under arctic conditions should be prepared. Compar-~ 
isons should be made with blast HOB charts since in general 
blast and thermal are competing nuclear damage mechanisms. It 
may be that considering the reduction that occurs in blast effects 
over snow that thermal will be of more importance in the arctic 


| even considering the low visibility conditions that can occur. 
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SFCTTON § 
NOCLFAR RADIATION 


> The nuclear radiations considered include the prompt 
gamma rays, the prompt neutrcns, the secondary gamma rays from 
neutron interactions with air and ground, ground activation 
pruducts, and the radiation from the fission products from the 
weapon. The last two components are typically treated together 
and considered in two time regimes. The initial radiation occurs 
within a minute or so following detonation while the residual 
radiation zis that which is contained in the rising debris and 
subsequentiy is distributed over a wide area e3 fallout. For the 
underwater bursts the initial radiation is associated with the base 
surge, and there will also be some radioactivity remaining in the 
water which should be considered when considering possible ship or 
sumarine contamination. 


§.1 Arctic Environmental Differences 


es fhe primary atmospheric parameter affecting the prompt 
radiation dose is the density. As shown in Section 1.2 the 


January 75° standard atmosphere has a density 16% greater than 
the midlatitude standard typically used for weapons effects 
studies. Por the temperature extremes noted in Section 1.2 the 
density will be even larger. For these cases the radiation will 
tend to be decreased relative to temperate areas. 


| af The ground composition can have an effect on the neutron 
and gamma ray transport in the atmosphere primarily involving the 


secondary gamma rays. The amount of water in the ground is 
important. 


Under arctic conditions involving a snow or ice cover 
changes might be noted in the neutron and gamma ray dose. 
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| earl Fallout depends upon many parameters which are signif- 
antily different in tne arctic. The particle size distribution 


and compcsition of the swept up debris cloud will be significantly 
different for the snow/ice Situations. The atmosphere profiles of 
pressure, density and temperature may change cloud rise character- 
istics. The meteorological conditions in the arctic including wind 
and precipitation patterns could affect the fallout distributions. 


The freezing conditions that occur in the arctic area 
may be important in terms of retaining radioactivity from the base 
surge on ships near the area. 


5.2 WP From Radiation Effects 


oe The characteristics of prompt nuclear radiation under 
ac ¢c conditions will be discussed with regard to effects of the 
air density, the ground composition, and the depth of burst. 


5.2.1 | St: Density Effects 


P| The vast majority of the predictions that are made of the 
effects of prompt radiation use scaling relationships applied to 
infinite uniform air transport results. The techniques can involve 
codes such as ATR (Harris, 1972) or graphicai techniques as con- 
tained in EM-l. (DNA, 1978). 


ay The basic transport resultsS are typically presented as 
a 4WR° dose as a function F(PR) of the amount of air between 


source and receiver (PR). F(pR) will depend upon the particular 
source spectrum of interest and the particular dose response func- 
tion desired. The dose at a particular range R is then given by 


F( PR) 
4eR 


D= 


where for a uniform density (p) case the amount of air is just 
the product of pand R. If the density varies over the path, an 
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integral of tne density over the path is used in the above expres- 


¥ 
: 


sion. At the same range in atmospheres of different densities the 


a dose will be related by the expression 

z Dy a F(P,R) 

. Dy) F(P,R) 

* The function F(PR) typically shows a buildup from zero as con- 


tributions from multiple scattering increase the dose, then a 
decrease with increasing pR as various extinction processes begin 
to dominate the transport. ; 


: At ranges for typical environment levels of interest 


ene ee ENE A ee en ET A Saar SS A - 


near sea level the F(pR) is a decreasing function of PR. Thus, 

if P, > Py then the dose Do is less than Dj. Thus, for the winter 
arctic conditions where the density is greater than in temperate 
climates a decrease in the radiation will be expected. 


Lg Example curves are given for several different prompt 
radiation doses of interest in military systems. In all caces a 
1 MT nominal thermonuclear weapo.. is assumed for the source function 
for the radia’ion. The neutron fission heating is used to assess 
the vulnerability of warheaus. In Figure 5-1 tne neutron fission 
heating is shown as a function of range for several different 
ratios of the density to the density of sea level standard mid- 
latitude atmosphere. Ignoring surface effects, the top curve 
represents the neutron fluence versus range for the weapon chosen 
neac sea level for temperate climates. The curve marked i.157 
represents the neutron fluence expected for the arctic winter 
Standard case. Note that for a 1025 
about 10% is noted in the range. The higher ratios refer to 


heating level a reduction of 


: densities corresponding to more severe coid temperatures that 

might occur in the arctic For the extreme case corresponding to 

‘ a temperature of about -80° which occurs very rarely a reduction 
of about 21% is noted. 
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es. The neutron fluence expressed as the 1 MeV Silicon 


bility of systems. Figure 5-2 shows the 1 MeV Si Eq fluence 
For a level 
; of 1022 the nominal arctic reduction is about 12%. The 

: extreme difference at 10}2 is about 25%. At higher levels 


versus range for tne density ratios considered. 


corresponding to vulnerability criteria for harder systems, the 


: difference is slightly less. at 1024 


s are 10% and 22% respectively. 


The prompt gamma ray dose rate is a common damage mech- 


range of densities considered previously and for yields of 10 KT 
rad (Si)/sec 
and fo’ the extreme 


. and 1 MT. For the 1 MT case and for a level of 
the reduction for the arctic winter case is 11% 
case the reduction is 23%. For the l0KT burst, 
reductions are 10% and 20%, 


wi In Figure 5-4 the total dose from the 
rays is shown for a 1 MT source for the density 


For the 10° level a reduction in range of about 
the arctic winter case and a reduction of about 
the extreme case. 


Note that all reductions found for the arctic winter 
case are about 10% and the reductions for the extreme caSe are 
between 20% and 25%. These correspond to area coverage reductions 
of 20% and 38% to 44% respectively, which might not be insignificant 
for specific system considerations. Note also that if one enters 
; the curves at a certain range, tor instance 1.2 km on Figure 5-2. 
* the 1 MeV Si Equivalent fluence decreases from about 3.5x1 


* temperate climates to 1.5x19013 


changes in the fluence levels and could easily span the range 
i from sure kill to sure safe for a system. 
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192° 


the correspondina 


prempt gamma 


ratios chosen. 
9% is noted for 
19% is noted for 


for normal arctic winter to 


Equivalent fluence is more commonly used to assess the vulnera- 


the reductions in range 


anism for TREE. In Figure 5-3 gamma dose rates are shown for the 


913 


: 3.5x1024 for arctic extreme weather. These are very significant 


OL Ws 


for 


Ne ee. ce cet eal eel ERE wae an SY | 


aa 


+ 


° 
_ 


oe en eT 
aie ne te te te 


a see 


-_ 
t—) 
z 


—_ 


NEUTRONS/CM 


% 


0 2 4 6 8 1.0 1.2 1.4 
Range (km) 


Qpicu 5-2 SILICON EQUIVALENT NEUTRON FLUENCE .- 
AS A FUNCTION OF RANGE - 


ee coded OO ee PRIAFLLLA i. 


x é = 7 : ‘i ' + . + * 
Te ee eT FO FU RUA WR WPL BRUM WAU WPS We Br RO RAN BP AP BR da” tat. ba" or 


: ' 


: PORE. 
RP SS IE FETT TIT TOE, 

- a o_—_ 

> 

¥ oe 
# A 

4 pe 
if : 
é . 

‘ A 

. \ 

; fo 

‘y ke 

i ve ‘ 

¢ ws! 

¢) e w' 
5] 102 At t 
’ . 4 
4 ; WHA s 

oa wt 
Y “ i 


: i 
Roe : 10 }— 10 KT : 
2 7 
Ny) Cg 
2 ~~ 
Ns a t 
vi Pas 
: = 
Ny — 
2 
A 2 
0, 
{ oe 
* ww 
! 6 10 
oe 10 . 
? 
Y + 
‘ 
i : 
§ ; ; 
r) 
\ 
U 
- 
Y 10° 


ai 
“TT 


y 
’ \ 
\ 
f { y 
| 
e 
‘ t- 8 
4 : 10 
m \ 0 2 4 6 T.0 1.2 2.0 : 
% l Ground Range (km) i 
‘ 1 GP cure 5-3 PROMPT Y DOSE RATE AS A FUNCTION OF RANGE 
4 - ie : 
’ 5-7 . 
t) : t 
‘ ; | 
‘ | as e f 
d 7 t veep ente mee 0 cee ramets Le attain dee need y 
‘ Pe i ay i 2, faikeig: Gina Go ESS, Ne Sa ey OE tg oe ‘ 
0 ee sree. : 
.) q 
a $ 


UCU TAT eR CHER Ce eke tit hata testes Si itil Can ek On aL ail eS SAE SAN NTs 


ad OD A 


RRR at a CXF RTI TERUG wt UU CP Ue PU ON eM PL LN TKY 


- - \ 


106 


: % j05 
oe rn 
(a) 
& 
& 
10° 
: » 
103 . 
225 0 2 4 ie.) 8 1.0 1.2 1.4 #1.6 #$%J1.8 2.0 2.2 
Ground Range (Km) 
yaa 5-4 PROMPT GAMMA DOSE AS A FUNCTION OF RANGE 
4 ry i 


; 
« * , ie > ft iy inn ee 
x é ee ee rer Tey cee sab ie tN Si att enn lai aha Re Mead / a, 
a | : : ~~ aoe 
Ae j Y : : \, oes be cae 
~ ? af . as ae 


! : igs pou, . 5 ae CT hw a) Pi ne OP Ry MF wa ev ia ia” ate LY bgt ia’ ae oat in” ae” a ne “th 
ial rene aaa ote Set as te MeN TY Ha ohn <, aN AA on Soar Mn eee es as “st,” a hy 


a a a en, 


| a The initial fission product radiation depends upen 
e density in a more complicated fashion than doves the prompt 
radiation. The time scale is such that the atmosphere is signif~ 
icantly perturbed by the weapon. The fireball within which the 
debris is contained is growing and rising. Tne shock wave is 
moving outwards through the air altering the integrated density. 
Several techniques have been derived to handle this ccmponent of 
the radiation dose as described in EM-1 and computed in ATR. Es- 
sentially, infinite air results are scaled to account for the fact 
that there is a hydrodynamic enhancement *ecause of the low density 
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fireball region. If the initial dose is ¢.sired, scaling is a 
difficult tack since the time dependence of the various effects 
must be considered and integrated. The above two techniques 
use a scaling method due to Mooney and French (1965) to incluée 
the ground effects on the dose. 


Another method developed (KSC, 1974) uses the density 
profile defined by the HULL code as a function of time and 
transports the gamma rays through the highly perturbed air by 
Monte Carlo methods. The HEAT code has been exercised several 
times to compute specific cases of interest for Ballistic Missile 
Defense cystems. A sample calculation is shown in Figure 5-6 
showing the tissue dose as a function of ground range for a 5 iiT 
burst at an altitude of 3.4 ft. The HEAT results are seen to be 
considerably less than the RRA results which were thought to be 
compatible with EM-l. However, note that there is significant 
difference berween the RRA model and EM-1. The 500-rad tissue 
dose level is a very severe personnel dose. At this level and 
below the difference between the HEAT and EM-1 results are less 
than 10% while the RRA results are much larger. 


el In comparing HEAT cases run, it was found that the 

ity scaling results using the integrated density froin the 
firevall to field point agreed with HEAT results especially for 
scaling between comparable cases. For the higher density arctic 
cases, the firebal] will be slightly smaller. Cube root scaling 
is used to define the ficeball radii with increasing altitude 
(decreasing density). If we use the scaling for higher densities 
also then reductions of the fireball size of about 4% and 10% 
would be seen for the arctic and severe density cases. These 
are within the uncertainties of the fireball modeling itself. 
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P| FIGURE 5-6 TISSUE DOSE VS GROUND RANGE, 5 MT, 3.4 KFT HOB 
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Therefore, the HEAT results have been scaled by the 


game scheme as before to show possible dose reductions in Fig- 

. ure 5-6, At the 500 rad tissue level, the reduction is about 13% 

7 for the arctic winter case and about 25% for the severe case. 
Again, these represent significant reductions in area coverage for 

1 the arctic cases. 


: §.2.2 Wpecrects of Ground Composition 


: Loe For sources and/or detectors near the ground surface, 

e radiation fluences are depressed below the infinite air 
results because there is absorption in the ground and loss of 
radiation from the atmosphere. The ground composition can affect 
the production of secondary gamma rays. A set of such ecalcu- 
lations (Campbell and Sandmeier, 1973) has keen made for several 
sources and several detector and receiver altitudes above the 
surface. Surface compositions of dry ground, wet ground and 
sea water were used to determine the effects of composition. 
Figure 5-7 shows the results as a function of slant range for 
a source at 5 m and a detector at 2.5 m above the surface. The 
results are presented as tissue dose, and the neutron and second- 
ary gamma contributions are shown separately. The neutron dose 
is seen to be much less than in the free ai- case and to be 
essentially independent of surface composition. The secondary 
gamma doses for the surface cases are larger than for the free 
air case at the snaller ranges and drop lower at longer ranges. 
Note that there is some variation in the secondary gamma ray dose 
for the three compositions and that the difference becomes less 
with increasing range. The dry ground gives the highest secondary 

i gamma dose and the larger the fraction cf water in the ground 
the smaller the secondary gamma dose. 
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Arctic soils are said to represent the same general 
classes that are observed in temperate climates; so no marked 
departure in the radiation doses would be expected. In the 
tundra and muskeg areas, the water content is higher than in 

- the temperate climates; so the wet ground curves would be 

more representative. In those areas with heavy snow or ice 
cover or over the open sea the sea water curve should be used. 
It is probable that the curves over fresh water would show some 
slight decrease below that shown for sea water because of the 
absence of the salt contribution to the secondary gamma rays. 
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In conelusion, the effect of the ice or snow cover on 


E 
es ne 


G 
the prompt radiation would not be expected to be large and would a 
generally be represented by the values that were computed for " oy 
Sea water. The wet ground results should be used for prediction i a 
purposes rather than the dry ground as being more representative 2 
of arctic econditons. : *¥] 
rat 
$.2.3 as Depth of Burst Effects : ft 
I£ the burst occurs below the surface of ground or <4 ne 
water, the prompt radiation from the device is strongly affected on rr 
by the surface material. Only a few feet of material is neces~ aa we 
Sary to markedly reduce the amount of radiation reaching the i 
atmosphere and being transported in the manner described above. - | f 
No particular differences in this effect would be expected in | us 
the arctic as compared with other underground and underwater i ra 
bursts. In the case of a burst beneath snow, the depth should es 
om be measured as the equivalent water depth since the density can t 
be much less than water. at; 


The initial radiation from the early time fission 
products can be an important contributor to the radiation dose 
even for bursts under the surface. The fission products and 
activated materials are ejected above the surface and form a at 
radiation source which may be highly anisotropic because cf the 
Surrounding surface material ejected into the atmosphere. Calcu- 
lations of this effect have been made for shallow-buried munitions “ 
in the ground, and no major changes would be expected for arctic os 

f ground condi:.ons. For those cases involving bursts in ice or 
snow, no comparable calculations have been made. An estimate of | 


cs 


zy +58 


the effect could be made by measuring the snow/ice depth as 
equivalent soil mass. There is some evidence, however, that for 


oe ae ee 


equivalent conditions a larger amount of snow or ice could be 
ejected into the air resulting in a reduction in the radiation. * 
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5.3 rg Hesidual Radiation Effects 


Residual radiation is that radiatic1a that is emitted 
later than one minute after the cxplosion. The sources and char- 
acteristics of this radiation vary depending on the extent to 
which fission and fusion reactions contribute to the energy of the 
weapon. Residual radiation from a fission weapon arises mainly 
from fissiom products and, to a lesser extent, from radioactive 
isotopes formed by neutron reactions in weapon materials and from 
uranium and plutonium that have escaped fission. Other sources 
of residual radiation hazard are the activity induced by neutrons 
that interact with ' arinus elements present in the earth, sea, air, 
or other substances :n the explosion environment. The most impor- 
tant of these sourc:.s is the neutron-induced activity in soils. 
The radioactivity f-:m a thermonuclear weapon will not contain the 
Same quantity of fission products that are associated with a pure 
fission weapon of c-e same yield; however, the large numper of 
high energy neutcons w.ll produce larger quantities of neutron- 
induced activity in weapon components and the surroundings. The 
total radicactivity from such a weapon will, however, generally 
be less than f£: 7m a pure fissior weapon of the same yield. 


ay The residual radioactive contamination (fallout) that 
results from fission ¢ sducts that are distributed subsequently to 
a contact Surface or subsurface burst is much greater than the 
radioactive contamination that results from the induced neutron 
activity. Thus, the neutron-induced activity may be neglected 
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for contact, surface, and subsurface bursts. 


ea If a weapon is burst in the transition zone (burst height 
<100w 35 feet) as far as faliout is concerned, the neutron-induced 


ee 


activity generally can also be neglected if the burst height is in 
‘ the lower three-quarters of the fallout transition zone, i.e., if 

the burst is below about 75W°*3° geet. If the height of burst is 
. in the upper eee of the transition zone (between about 754°° 35 
° feet and 1oow? 235 feet), the Heurronrangeees eee may net he 


mee x oi + . 


vj negligible compared to fallout. 
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5.3.1 oa Induced Activity 


| nd The type, intensity, and energy distribution of the 
nduceé activity produced by the neutrons will depend on which 


isotopes are produced and in what quantity. These factors depend 


on the number and energy distribution of the incident neutrons and 


the chenical composition of the soil. Induced contamination con- 


tours are independent of wind, except for some wind redistribution 


of the surface contaminant. The contours can be expected to be 
roughly circular. 


Examination cf several thousand analyses of the chemical 
composition of soils and the relative probabilities of neutron 
capture by the various elements present in the various samples 
has indicated that sodium, manganese, and aluminum generally wiil 
contribute most of the induced radioactivity. Small changes in 
the quancities of these materials can change the activity signi- 
ficanily. Other elements can also influence the radioactivity. 
Some elements have a relatively high probability for capturing 
neutrons (cross section), but the isotope that is formed after 
the capture either is not radioactive, @oes not emit gamma rays, 
or has such a long half life that the low activity does not pro- 
duce a hazardous dose rate. The presence of such elements in 
the soil will tend to lower the hazard from neutron-induced ac~ 
tivity. 


eal Calculations (Pugh and Galiano, 1959) have shown that 
se induced activity in sea water .s about a factor of 1600 less 


than in Nevada Test Site soil for times after burst of 1 hour or 
greater. At early times the contribution of the very short half 
life 28a) in the soil makes the ratio even larger. From the fact 
that sea ice has slightly less salt content of sea water and ice 
over land has a low salt or mineral content, one can assume that 


induced activity in an snow/ice layer will be less than in sea 
water. 
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Less induced activation should occur in ice/snow/soil 

configurations compared to bare soil due to the shielding effect 

} of the ice or snow layer. The magnitude of the effect depends 

s upon the depth of ice or Snow. In areas such as the Greenland 

7 ice cap, where ice is over 1000-feet thick, little or no neutron 

i" induced radiation should occur. The same result shouid prevail 
where snow packs of 30~50 feet may occur. Data to support 

7} quantitative conclusions with respect to ice and snow have not 

been discovered. 


5.3.2 Noe Radioactive Fallout 


A yet unpublished report (Spencer, Chilton and Eisenhauer) 
contains an excellent discussion of fallout gamma rays from nuclear 
detonations with exhaust‘ve literature references in all phases of 
the fallout problem. In these lists, there are no references to 
work involving the effects of the arctic environment on fallout. 


The source of fallout is a combination of the fissicn 
Products, weapon activation products and activation products from 
surrounding materials such as soil. For low altitude or surface 
bursts much of the activated materials will be vaporized or 
fragmented by the strong shock interactions and swept up into the 
rising cloud and will contribute to the total fallout dose. 


Fg At very early times, the fireball is nearly spherical or 
hemispherical for a surface burst and is beginning to rise as the 
. blast wave begins to move outwards but there has been no signif- 
icant movement of material. The mixing and entraining of the swept 
rr up soil materials is occurring. Calculational models of this 
development have been deScribed (Huebsch and Olken, 1976) using 
the HULL hydrodynamics code to describe the flow field from the 
low altitude burst. Routines are added which describe the growth 
of water, ice, soluble salts or insoluble particles in the rising 
- expanding cloud. It is possible to include the effect of different 
meteorological conditions in this model. No ¢ases have been run 
a for conditions approximating the Arctic. 


eee 
Tot te AM me ee ee ee a ne ee at ltl ate ln Bate tnt eanienlaal Be . 


oe tee 


rr 


ee 


Fa ee a tae Oe Wah Het Os? On" Fa? te fe" Ba eM Be be ta Bat 8 s7 Oat be! Bev oa” 6% Wad bev Geb Gav be wet gou gow ota baa te” Oa? 


A toroid.l circulation builds up with very large upward 
velocities and the gamma source is strongly radiating and moving 
upwards rapidly. While the movement is occurring the vaporized 
materials will be cooling and will condense into very small 
particles, If solid or liquid particles are present, the material 
will partly diffuse into the surfaces of the particles. 


For altitudes above which significant solid materials 
are drawn into the cloud only very small particles will be produced 
and there will be little localized fallout. In this case, the 
residual dose near the burst will be primarily due to the induced 
activity, if any. For lower altitudes of burst, large quantities 
of material will be entrained and there will be a wide range of 
particle sizes in the cloud. The larger sizes will precipitcte 
out to form the local intense fallout field and the smaller sizes 
will remain in the cloud for a long time ana may be dispersed 


over a wide region. 


F The spectrum of debris particles tends to be representa- 


tive of the soil composition and type. Thus, for bursts in the 

arctic over land no significant differences would be expected in 
the cloud loading except possibly a slightly smaller fraction 

of soil material as compared with water since the water content 

of the soils in the arctic is typically larger. 


A burst over snow or ice in the arctic would not contain 
any of these large solid particles and the average size of the 
cloud particles would be very small. This would lead to less 
intense local fallout patterns and a large amount of the radio- 


active material would be swept to high altitudes and widely dispersed. 


This case may be nearly the same as occurs with a burst over sea 
water where the cloud material consists of weapon debris, salt and 
water. The particles are extremely small but highly hygroscopic. 
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us, their size buildup can be very sensitive to local meteor- Ny 
Olegical conditions. The fallout from water bursts is described 
in EM-I. A much less intense local fallout is expected unless “4° 
rainost occurs. a 


The salt swept into the atmosphere may have a seeding 
effect and result in a weapon induced rainout of local material. 
Because of the lower temperatures and the very high humidity in 
the arctic, these seeding effects may be enhanced. The moisture 
Capacity of the atmospnere is much less and precipitation may be 


much more likely than in temperate climates. The particle growth Y 
in the clond may also be affected by the atmospheric parameters. 
No work in this area was found in the literature. y 
| 4g For low yield weapons, most of the moisture comes from i 
air entrained by the developing cloud. In the arctic, the absolute ° 
humidity is very low because of the low temperatures so that the 
water available for preducing the larger sized particles may be 4 
less. Thus, under these conditions there may be fewer large : 
particles produced and a less intense local fallout. x4 
Modeis have been developed which describe the rising 
debris cloud and its dispersal by winds. The diffusion of the a 
radioactive material and the influences of precipitation, * 
converging and settling of the particies by gravitation and 4 
diffusion are considered. The meteorologic parameters including es 
precipitation, wind patterns, and humidity have a strong effect } 
on the late time fallout. : 
a 
iad Much of the Arctic region is quite arid. Annual precipi-~ HW 
tation over the Polar basin, i.e., the Arctic Sea, ‘s less than six , 
inches of water equivalent per year. Most of the precipitation ie 
ce 
iM 
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alls as snow. There are local areas such as north of Hudson x) 


Bay and eastern Canada with moderately high snowfall (50-100 in.). ” i 
Lowest smowfall is in the northern Canadian islands and in north we 
i , Greenland, where total annual precipitation is frequently less 
than fomr inches of water. 
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{ a Ci Contrary to some popular opinion, surface winds in the . nd 
‘ retic are on the average very light. Observations from Soviet ‘ 


drifting stations in the central Polar Gcean indicate that monthly 


mean speed at the surface is about 8-10 knots. However, in well a“ 
\ developed storms wind observations show speeds on the order of 4 
50 knots. The mean wind speed increases rapidly with altitude ] 


and just below the tropopause (7-8xkm) the highest monthly mean 
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wind speed may reach 40-50 knots. Maximum wind speeds may be a 
om | much higher. Wind variability is larger in the Arctic. 
! ‘ae ' 
! The fallout prediction models range in complexity from x 
ae empirical fits to fallout measurements made on the U.S. nuclear ss 
oo test series to very sophisticated numerical models which attempt NN 
to describe the development and dispersal of the cloud from first * 
~ principles, WSEG is an example of the first type of code which it. 
has seen wide usage in system codes because of its fast running 
time and realistic results. The yieid of the burst and the wind 
‘ 
speed description are the major input parameters. There is no Hl 
t 
provision for adjusting to other meteorological parameters. . 
< 
F oe 
i) : > one NUCROM (Baum et al, 1974) is an example of a code with it, 
intermediate complexity. It is a simplified rainout model which 
= 2 ° 2 < r 
1 tt allows some freedom in introducing meteorological parameters. A ‘ 
, stabilized debris cloud model is used which is then separated into 4 
‘ : segments as a functic. of altitude. Diffusion and migration under 4 
wre the influence of wind are considered and scavenging by precipita- 5! 
. tion events is allowed. The scavenging efficiencies are nandled ~ 
‘ xs in a gross manner and detailed particle distributions are not i, 
; ' considered. a 
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DELFIC (Maloney and Klemm, 1975) is an example of < 
code which attempts to calculate nuclear fallout from basic 
physical principles without recourse to empirical modeling of 
the test results. Detailed calcuiations are made of the debris 
cloud rise and loading by particulate matter. Particle size 
distributions and their evolution with time are computed. A 
detailed atmospheric definition is used. The distribution of 
activity with particle size is considered. Precipitation clouds 
are defimed and detailed descriptions of the in-cloud and below- 
cloud scavenging as a function of debris particle size are given. 
This code is a long running, expensive code to use but is capa- 
ble of including the full range of Arctic meteorological param- 


eters. 


Representative calculations were made (Normant, 1974) 
of the scavenging by rain and snow clouds for tactical nuclear 
conditions in Germany. In these cases the burst altitude is 
high enough to minimize the local strong fallout and the activity 
is contained in micron-sized particles which would disperse over 
a very wide area with a low intensity. Scavenging by precipita- 
tion can, however, under the right conditions produce a very 
intense local fallout field. The effect is most important for 
low yields because the stabilization altitude for large yields 
is higher than normal cloud altitudes. 


In Figure 5-8 the rain and snow washout coefficients 
are compared. The precipitation events were picked to fit 
European conditions with the rain precipitation rate being 
20 mm/hr and the snow being 10 mm/hr. For the larger particle 
sizes the rain coefficient is seen to be almost an order of 
Magnitude larqer than the snow but may be comparable in the 
micron region. 


5-23 


in a, tne ee see ne a a 


NRURRUU TRI STVAERO LST REE Dab ole" oFe ole Miat ica g' a8 6 Rp Wel Nal Ol Ol 


“7 


oye", 


ee 


ee 


et 


ee 


ee 


to 


ee ee 


Washout Coefficient, A (sec ') 


r) 5 10 15 20 
Particle Diameter (pm) 


pice: 5-8 COMPARISON OF RAIN AND SNOW WASHOUT COEFFICIENTS 


PELL Le 


Pa ie 


a ee 


5 Ka 


» 
- . Sa Oa woe ied tec ae ae Co Go Cal Cel Sal =. ~p =a "pm R® aa” he “" ." au" iad eal ." oan, 
fag ON, ie eee eee te Ned +} al Ay Pe ere wi ne ~ aa? ide ag? Fagtoe, Patout 


bad Pies bad i. nf Oe 


5 y Sess ey e : e 
cs, We we ee ee RP a°e, eo tea * Bet atl 


i a Se, on 


ee 


bie"2 "2%" 


In Figure 5-9 the fraction of particies scavenged 
is noted for interactions of 10 minutes and 1 hour with rain 
and snow clouds. Note that for particle sizes above about .3 
micron all of the particles are scavenged even in a 10 minute 
interactiam with a rain cloud. The critical diameter is about 


jo FOS ay 


-03 micron for snow, and a large proportion of the particles 

q = with sizes above this diameter are scavenged for a 10 minute 

h : interaction. Normant’s calculations of the snow scavenging 
_ efficiencies Go not agree with other calculations and experiments 
: as will be described in Figure 5-11. 


fhe main point of these results is that a fairly short 
‘ interaction of the debris cloud with a precipitation event can 
result in essentially complete removal of activity from the 
cloud and deposition on the ground with the precipitation. 
Calculations made with typical Arctic precipitation rates would 
{ be of interest. Considerations of induced precipitation events 
ib and the interaction with the debris cloud would be of great 
if interest and may represent the major difference expected in the 
X Arctic. 


In the previous figures the submicron particle size 

is seen to be a region where the scavenging efficiency shows a 

i large dependence upon particle size. This is due to the detailed 
physical interactions that are occurring and their relative 
importance. For very small particles Brownian diffusion is 

i ‘ very important, and for large particles inertial accretion of 

particles dominates. The interactions are very complex for micron 

size particles, and electrical forces may play a major role. 


An extensive study of scavenging was performed at 
t Illinois Institute of Technology Research Institute (Knutson, 
{. 1974) with theoretical and experimental studies of the relative 


rain and snow scavenging efficiencies being considered. Repre- 
| H 1 sentative results are shown in Figure 5-10 where the snow 
¥ ' 
he Se. 
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is FIGURE 5-10. EMPIRICAL MODELS FOR PARTICLE SCAVENGING BY SNOW AND 
RAIN (KNUTSON, 1974). 
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= is predicted to be much larger than rain for 


’ 
particle sizes above .5 micron. This contradicts the results | 
given in Figures 5-8 and 5~9 and represents the current uncer- : | 
tainties in these types of data. 


A recent extensive review (Pitter, 1976) considers 
snow and ice scavenging and the detailed physical interactions 
that occur. In Figure 5-11 results from various investigations 
are compared for ice crystal scavenging. Very wide variations 
are noted for the submicron to micron size range. The left 
hand portion of Figure 5~l1] summarizes theoretical and laboratory 
results for the snow scavenging coefficient. The Pitter (labeled 
present results), Starr and Mason, and Sood and Jackson results 
all show general agreement with the minimum in the scavenging 
efficiency occuring at about .5 micron. The Knutson results 
are considerably larger at the .5 to 2 micron size range but 
are decreasing rapidly at .2 micron with no indication cf an 
increase at smaller sizes. The DELFIC~PSM results equivalent 


wee 


to the results given in Figure 5-9 are also shown and are seen 
to be drastically different above .05. An abrupt increase in 
the efficiency at .05 micron is noted then no variation with 

increase in particle size. The differences in the results are 


ee eres ae 


due to different ways of handling the physical processes between 
the particles and the snow flakes. Also noted on the right hand 


i} 


portion of the figure are the results of various field measure- 
ments of aerosol scavenging coefficients. 


ae The main point emphasized by this figure is the current 
extreme variation in values of the scavenging efficiency for snow. 
The micron particle size is interesting for airbursts and surface 
bursts over water, snow and ice. This is precisely where a large 
uncertainty exists, and the differences noted would cause a large 
difference in the local fallout from a burst when natural or 
induced precipitation was occurring. 
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§.3.3 Underwater Bursts ny 

§.3.3.1 t 


Whe sources of radioactivity from an underwater ; 
nuclear esplosion are (1) the fission product activity in the Oi 
column an@ crown, or plumes, (2) the base surge, and (3) the 
residual zadioactivity deposited in the ocean, or radioactive 
pool. These phenomena are described for a burst in temperate Ma, 
regions im Chapter 5 of DNA EM-l and in Chapters 7, 9, and 10 wy 


* 
of the Underwater Handbook, and methods for predicting the isd 
Magnitude and duration of the effects are presented. Existing be 
Manuals, towever, do not address the modifications to these MS 
effects that might result from an underwater detonation in | < 
Arctic regions. i 
fhe Arctic environment can affect the sources of 
radioactiwity from underwater nuclear explosions in several n, 
ways: yi 
@ Ice cover may modify the characteristics of the Y 
surface phenomena and attendant radiation fields. These 
phenomena are variable at best, depending as they do on the N 
state of the bubble as it reaches the surface. Part of the | nN 
energy remaining in the bubble is expended in breaking through tn 
solid ice or imparting upward motion to blocks of ice. s] 
= ° Depending upon the depth of burst, an underwater | aT 
explosion that vents may form a radioactive column, plumes, or M4 
base surge. The cold Arctic temperatures will cause freezing ' wt 
of some of these products, producing a local fallout field on | 
the ice and the formation of radioactive ice on the weather | ) 
decks of ships close to the detcnation. Ne 
@ The presence of ice and the typical Arctic water 
density wrofile may affect the formation and migration of the MY 
radioactive pool, which normally rises to the surface and hg 
diffuses fairly rapidly. Heavy ice cover may contain the pool Sa 
Wa 
‘y 
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from a veating explosion to the area cleared by the explosion. 
In the case of a very deep, non-venting explosion, a solid ice 
cover may contain the pool below the ice so there is no above- 


surface evidence of its existence, 
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Wnere the extent of the ice permits operations by 
Surface ships, those operating in the vicinity o: an underwater 
burst may encounter a buildup of radioactive material on the 
superstructure and weather decks due to the freezing of the 
base surge and spray from the radioactive pool. ‘thirasawa and 
Bjerke, 1968 studied this problem in some detail, using the 
same computer model (DAEDALUS) used to compute dose rates and 
total dose for the various conditions presented in the examples 
of Figures 5-46 through 5-75 of DNA EM-1. The report of 
Shirasawa and Bjerke is the source of the material in the 
remainder of this section, including the figures. 


The factors considered in the study include the 
transit radiation exposure due to radiation emitted directly 
from the base surge or the contaminated pool, and the deposit 
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‘radiation exposure due to emission from fission products after 


they have settled on ship surfaces. The transit exposure esti- 
mates were based on the CAEDALUS model output (exposure rate) 
with the assumption that the Arctic environment does not alter 
the radiation characteristics of the base surge and radioactive 
pool. The deposit exposure estimates were based on the rate 
and extent of shipboard icing and the concentration of fission 
products. A number of previous studies were consulted to 
determine the rate of formation and distribution of ice and its 
associated fission products from the base surge and pool. A 
typical destroyer was selected as the representative ship for 
the calculations, a nominal 10 kt ASW weapon was chosen, four 
explosion depths were selected (65, 150, 500, 1000 ft in 

5000 ft of water), three post-detonation entry times (10, 20, 
30 min), and three ship transit speeds (10, 20, 30 kt). 


The rate of ice increase is a function of temper- 
ature and wind speed. Figure 5-12 is a semi-quantitative 
presentation of the relation of these parameters. The various 
regions surrounding the conditions for icing shown in Ficure 
5~1 indicate the following: 


Region I Wind force is not sufficient to biow 
spray over ship. 

Region II Temperature is not ow enough to cause 
spray to freeze on surfaces. 

Region III Wind force is so high that green water 
covers decks and keeps melting ice (ice 
is possible on higher superstructure). 

Region IV Temperature is so low that spray 
freezes before striking ship. 


The central area of the diagram where a “heavy” icing 
rate is indicated corresponds to 2 tons per hour or greater. 
The “iight’ area denotes rates of i ton per hour or less. i 
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Figure 5-12. Icing Rate as Related to Wind Force 
and Temperature 
’ 
; _- (Shirasawa and Bjerke, 1968) 
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; a ; this presentation, the situation described for Region I is the by 
io. least clear. However, in the studies cited no icing was af 
a ; obserwed when the wind force was less than Beaufort Force 3 in ry 
\ * spite of the low temperatures. f. 
eI , Se Serious icing occurs where temperatures below 29°F “ 
NY ‘ -1.5°€) are combined with flying spray, which forms only i, 
wf . : with winds of 17 kt or more (U.S. Navy Hydrographic Office ie 
i F Pub. Bo. 705). The areas within which these conditions occur Ky 
ee are mere restricted in latitude than is generally realized, My 
eae owing to the modifying influence of water temperatures on i 
ce surface air temperatures; however, a significant proportion of ! 
at ocean operating areas may be subject to serious icing condi- a 
tions. Figure 5-13 shows areas of probable superstructure fi 
*% icing for January-March, May, and November, based on a southern st 
limit of 10% frequency of temperatures below 32°F (0°C). " 
, | Shireésawa and Bjerke present the results of the 
*S computer calculations for a number of combinations of param- in 
‘ eters, but in view of the conclusions of the study these will MY 
not be given in detail here. To provide for maximum ionizing ¢ 
. radiation exposure, they assumed an early (10 min) entry time 
PAs: into the base surge followea by a traversal of the pool, under a 
a no-wind no-drift condition and with concentric base surge and ng 
a pool still underyoing dynamic expansion during the traversal. Ne 
Figures 5-14 and 5-15 summarize the calculations and present a om 
- comparison between the deposited and transit contributions to fe 
. the total exposure, for a 65 ft and a 500 ft depth of burst > 
respectively. It is immediately apparent from the figures that 
a ' the activity entrapped by ice accretions, regardless of source, 
snk a is not a major contributor to the total radiological hazard. 
; In each case, the exposure contribution made by base-surge 
wee fe deposit is only 8% or less, while the pool spray deposit is 
: z Kao negligible. 
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The discussion and conclusions of the Shirasawa and 
Bjerke report are quoted: 


aa "In our preceding consideration of contamination 

by ice-entrapped fission products, the possibility of 
countermeasures has been wholly ignored. This was 
purposely done to permit a maximum hazard evaluation. 
However, it is obvious that several immediate possibil- 
ities exist for reducing degradation of personnel and 
ship capability following contamination of weather 
Surfaces. Three immediate countermeasures which might 
be considered are: 


1. Reduction of number of personnel in high 
exposure areas and rotation cf personnel. 

2. Use of ship washdown system, 

3. Initiation of ice removal procedures. 


The importance of countermeasures is borne out by 
consideration of the exposure rates existing on ship- 
board after traversal maneuvers, because of continued 
exposure from deposited activity, however small. 
Implementation of countermeasures will also prevent or 
Minimize contamination ingress as well. 

|g "Operation of the washdown system aboard a 
destroyer was shown to be feasible in freezing 
weather.* Initiation of this countermeasure upon 
leaving the radioactive pool would contribute to a 
significant reduction in the ice deposited by pool 
Spray and/or base surge contact. Though icing may 
continue during the use of the washdown system, the 
relatively warmer water from the sea would serve to 
melt and rinse away the contaminated ice accumulated 
@uring the pool and/or base Surge traversal. It was 
estimated * that use of the system for &0 minutes or 
more under conditions of an air temperature of 10°F 
and a wind velocity of 21 knots would produce a max- 

_imum of l-ineh of ice. This of course wouid be "clean" 
ice. It has been estimated that 6 inches of ice, or 
an ice accumulation of 200 tons, on hori: iital and 
vertical surfaces would interfere with the operation 
- of a destroyer in an 80-knot beam wind. 

ie are “Removal of slush ice after washdown cessation 
can, if ship mission and stability permits, be accom- 
plished quite successfully by personnel with shovels, 
brooms, boards and buckets. These procedures would 
effect the most direct and efficient removal of 


Editor's Note: The reference is to a report by 
Perkins, W. W. and Railey, R. M., Operation of 


Shipboard Washdown in Freezing Weather, U.S. Naval 
Radiological Defense Laboratory USNRDL-TR-972, 31 


December 1965, Unclassified 
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contaminated ice. However, these procedures may 
create the attendant problem of personnel exposure and 
the potential hazard of tracking activity inside the 
ship with subsequent danger of ingestion. 

"The following conclusions would appear to be 
justified within the general limits of this study. 


The radiological consequences to naval ships 
of coming in contact with the post-detonation 
formations typical of underwater nuclear 
explosions are not siqnificantly changed by 
an arctic environment. 

Radiation exposure resulting from freezing 
spray of radioactive-pool derivation does 
not present itself as a problem insofar as 
interference with the tactical missions of 
ships is concerned. 

Radiation exposure from the freezing spray 
ef base surge aerosol exceeds that of pool 
spray deposit, but it is well pelow levels 
which would threaten degradation of ship's 
effort. 

Initiation of ship washdown operation and/or 
Manual removal of slush ice can reduce the 
amounts of deposited fission products to 
levels comparable to those existing in more 
moderate environments. 

The limiting radiation hazard for involve- 
ment in poSt-detonation ship maneuvers wili 
be the transit exposure as a result of 
encountering the base surge and the pool. 
There is no reason to believe that this 
exposure will be significantly changed by an 
arctic environment,” 


It may be concluded, on the basis of the foregoing, 


that the estimates of Figures 5-46 through 5-75 of DNA EM-1, 
prepared for use in temperate climates, may also be used for 


Arctic environments. 


5.3.3.4 


GPpissicactive Pool 


Chapter 10 of the Underwater Handbook contains a 


' -. technical review of the literature as of December 1966 


on the distribution of the radioactive debris and associated 
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nuclear radiation from underwater nuclear explosions. It was 
concluded at tha: time that no adequate comprehensive radio- 

- logical prediction system existed ir the literature. With 

. respect to the radioactive pool, a review of the literature 

: Since that time reveals little reason to alter that conclusion 
Significantly. Rinnert, 1967 and 1968 has developed FORTRAN IV 
computer programs to estimate the exposure rate history and 
total exposure for surface and subsurface traversals of a 
radioactive pool, but these are based on the pool model of 

; Ksanda, 1963 and the work of Pritchett, 1966, both of which 
were available when Chapter 10 was written and were referenced. 


= The Rinnert reports are concerned mostly with the 
ocumentation of the programs and do not present results of 
calculations for ranges of input parameters. They each have an 
example, however, of information that can be derived from the 
programs. These examples are presented here, since neither DNA 
EM-1 nor the Underwater Handbook contains estimates of exposure 
for a submarine traversal of a radioactive pool. Figure 5-16 
shows the calculated exposure rate history for a single set of 
parameters, and Table 5-1 Shows the total exposure for traverses ~ 
as calculated Ly the modified Ksanda model (Rinnert, 1967) and 
by the modified Fritchett model (Rinnert, 1968), for several 
sets of parameters. Both examples are for unshielded detectors. 
The submarine's hull and internal piping systems would reduce 


these exposures by varying amounts, which may be calculated 
from standard references (e.g., DASA 1892). 
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Examples of Exposure Rate History of Un- 
shielded Detector Traversing Radioactive Pool 


Bands indicate range of estimates for circular 
pool and for elliptical pool whose minor axis 
is half the major axis. 
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Gen 35-1. TOTAL EXPOSUKE FOR SUBMARINE TRAVERSE 
OF A RADIOACTIVE FOOL 


(Rinnert, 1968) 


¥ Total Exposure for Traverse (Rcentgens) 


Speed of Parameters* 
Traverse dnsdte tesa Scherr Modified Ksanda 

(Knots) Model Model 

5 4.19 5.05 - 7.30 = 
10 5.30 6.29 - 9.03 DOB = 300 
15 6.22 7.13 - 10.2 SO = 4000 

20 6.89 7.79 - 11.2 2zD = 50 

30 8.28 8.82 - 12.7 
40 9.27 9.54 - 12.8 

5 0.48 0.89 - 1.26 Ww =l 
10 0.52 1.10 - 1.55 pop = 300~=CsS! 
15 0.59 1.20 - 1.70 so = 4000 | 
20 0.677 1.27 - 1.79 2D = 300 ' 
30 0.12 1.37 - 1.92 
40 1.7 x 10727 1.43 - 2.01 

5 14.4 15.7 - 22.9 w = 100 
10 17.7 19.3 - 28.3 DOB = 949 
15 20.3 21.9 - 32.0 SO = 10,000 
20 22.2 23.9 - 34.9 2D = 50 

30 26.2 27.0 - 39.5 
40 30.8 29.5 - 43.1 

5 1.68 2.81 - 4.07 Ww = 100 
10 1.89 3.58 - 5.19 DOB = 949 
15 1.97 3.89 - 5.63 SO = 10,000 
20 2.11 4.09 - 5.92 zD = 300 

30 2.58 4.35 - 6.30 
40 4.12 4.53 - 6.56 


WwW Yield, kilotons 
DOB Depth of burst, feet 
so Stand-off distance, yards (traverse beygins at stand~off 
distance at time of burst and proceeds across pool, 
passing through surface zero) 
2D Depth of traverse, feet 
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It should also be noted that the DAEDALUS calculations ‘ 
referred to in the previous section, presented by Shirasawa and 
Bjerke, 1968 and in Chapter 5 of DNA EM-1, provide predictions 
of the free field radiation 15 ft above the ocean surface, and 
therefore consider radiation coming from only a thin surface 
layer of the pool. There is little experimental evidence to 
confirm model predictions of the rate of growth and migration of 
a radioactive pool, of the fraction of its radioactive content 
left at warious @>_*4; as the bubble oscillates and migrates 
toward the surface ~ Subsequent history of deep pools ieft 
behind, or the corn -ci8S under which a substantial amount of 
radioactivity may ne--: appear at the surface. 


| od In view of the uncertainties surrounding radioactive 

pool foraation and behavior under non-Arctic conditions, and the 
oat lack of experimental data of this nature under Arctic con- 
ditions, an assessment of the possible effects ot the Arctic 
environment must be regarded as conjectural. Two character- 
istics of A.ctic regions that might alter the behavior of the 
radioactive pool from what might be expected elsewhere are the 
presence of ice cover and the strong density gradient in the 
water coiunn. 


= An underwater explosion at a depth that vents will 
normally cause the development of a radioactive surface pool 
initially centered at surface zero. In the case of scattered or 
broken ice, which is usually only a few meters thick, it is 
expected that the pool would undergo its normal expansion and 
diffusion, and there would be little effect of the ice except to 
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provide a certain amount of shielding against that portion of 


the radioactivity that is below it. In. the case of consolidated a 
: pack ice, with extensive pressure ridges and ice keels, the w) 
. venting explosion will blast a hole in the ice and the radio- Mt 
e active pool will initially be centered in the ice-free water of a. 
“ the hole. Since a 10 m pressure ridge will be accompanied by 
Se about a 50 m ice keel, and in extreme cases ice keels may extend be 
to 150 m (Bowditch, 1977, Chapter 36), the horizontal migration we 
; and diffusion of the pool may be impeded. If such should occur, AM 
° the radioactivity in the exposed pool would be smaller in extent " 
and more concentrated than might be expected from the DNA EM-1 A 
: examples in Chapter 5. However, since it appears that solid ‘ 
pack ice is necessary for this condition to arise, the effect ttt 
would be apparent only to aerial observation, and of signif- bt 
icance only to the determinaticn of the location of an under- 

water burst sometime after the fact. a 
As is discussed in Section 7, it is not known how much '*) 
ot the explosion energy is required to break through solid ice o 

c.7et and vent However, for very small yields or very deep 4 

ex s\losions, »t may be that the bubble will have insufficient x 
erergy when 4. ceaches the surface to fracture whatever thick- wy 
ness of .ce is gresent. In this case a radioactive surface pool ay 

wizl not be formed and the pool will be trapoed below the ice 5 

layer. Tiis would prevenc the detection of the explosion by 5 
aerial survey methods, although the radioactive pool would 4 
remain a hazard to submarine operations. m4 

In Section 1.2.7 it was noted that, in general, strong a 

- ° positive de: sity gradiznts exist in the upper few hundred meters it 
: of the Arctic water column. This region (pycnccline) severely Ra 
« impedes the voward migration of heat and salt and effectively m4 
‘« insulates the surface from the water masses below. This charac- mm 
teristic of the region leads to the speculation that much radio- x 

] active debris may often be trapped below the surface, whether or , 
not an ice layer is present. AN 
! 4 
4 
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ss Qucting from Chapter 10 of the Underwater Handbcok: 


“For deep and very deep explosions, where the 
bubble experiences several oscillations as it 
migrates toward the surface, radioactivity may be 
ejected from the bubble at minima .... Measurements 
at Operation Wigwam ... indicate that there are bcth 
a radioactive surface pool and random lens-like poois 
of debris in the thermocline layer. These deep 
pools were measured some days after detonation and 
were found to be small and quite stable. Whether 
these deep deposits represent radioactivity that was 
left behind by the migrating bubble or material 
Carried to the surface by hydrodynamic flow and 
returned to its original stability itevel, is not 
known." 


Except for the shallowest evplosions that vent most of 
their fission products to the atmosphere with a water column or 
Plumes, ard then fall back to form a radioactive pool with sur- 
face waters, underwater nuclear detonations in the Arctic will 
Cause a substantial amount of the highly saline, deeper water 
to mix with the radioact‘ve material. It is conceivable that 
this water mass with its crapped fission products will sink 
to and remain in the pircnucline. In addition, any radioactive 
pools left at depth by the pulsating bubble would have their 
ascent stopped at the pycnocline. Thus the radioactive effects 
of the surface pool could be substantially less than those 
predicted by existing models. The pools would remain a sub- 


Marine hazard, however. 
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On the basis of the work of Kaulum and Bennett, 1971, 
it may be concluded that there are combinations ot yield and 
‘ depth that give a high probability that the radioactive debris 


from an mCerwater nuciear explosion may be contained beneath 
the surfece for periods of time long enough for the radio- 
activity to decay to undetectable levels. Figure 5-17, besed 
On calculations for a wide range of density gradients, includ- 
ing typical Arctic gradients, provides a reasonable basis for 
estimatiags the conditions under which a radioactive surface 
pool would not be formed for yields of 100 kt or less, whether 
or not ice cover were present. Ice cover would prevent such 
formatica fur any detoration deep enough not to rupture the 
ice. The subsurface pools would, however, be a hazard to 


Submarines. 


Le There is no reason to expect any changes in the radia- 
ion darage vulnerability levels in Arctic conditions. The 
possible exception is some slight enhancement of effects on 
personnel. The severe winter Arctic environment imposes a 
heavy strain on personnel at best so that radiation effects 
might have a more deieterious reaction at lower levels. 


Bunkers or personnel shelters buried under snow or 
ice would provide slightly better protection than concrete on 
an overburden weight basis. Information on the protection 
factors is widely available such as given in EM-1l. 
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Consinusions and Recommendations 
Conclusions 


No studies considering the effect of arctic environment 
on prumct raciacion envircnmencs were found except the Ft. Bliss 
study fOSWD, 1969). However, the techniques that have been 
develoged to compute radiation environments in temperate environ- 
ments can be used with no basic changes except using the proper 
mocel atmosphere. The density is the only important parameter 
of interest. 


<a Scaling the available infinite air transport results 
to the arctic winter ground level density indicates that the 
environment levels corresponding to typical damage criteria for 
hardened electronics occur at about a 15% smaller radius under 
arctic conditions than under temperate conditions. At a par- 
ticular range the fluence or dose level can be 1/2 to 1/3 as 
large for arctic winter conditions as in temperate conditions. 
Thus, prompt radiation effects tend to be depressed in the 
arctic which is an advantage for considering damage to U.S. 
instailations from Soviet bursts. However, the reduction in 
prompt radiation effects should be considered when considering 
the effectiveness of U.S. bursts against Soviet systems. 


ey The presence of the surface layer under the atmosphere 
tends to reduce the radiation environment in the air as compared 


with the free air values, No calculations of this effect have 
been made for arctic surfaces and conditions. Inspection of the 
available calculations indicates that there is essentially no 
difference in the neutron dose as measured close to the surface 
for wet or dry ground or sea water. The dose from the secondary 
gamma rays for wet ground is about 20.6 lower and for sea water 
is about 30% lower than for dry ground. The dose over fresh 
water or ice might be somewhat lower still. Since most of the 
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oa 4s 
° retic is covered by wet tundra, fresh water ice or sea and 


sea ice, one would expect a small reduction in the dose result- Nt 
. . ing from neutrons from this effect. 7 
: i The dose from the early time fission products can be 
. important in contributing to the total dose received by reentry it 
vehicles, airplanes and ground installations as well as an 2 
‘ important contributor to personne! casualties. In addition to fi 
- the reduction due to the increased density there would be some re, 
effect due to the smaller fireball. No calculations of this “ 
effect have been made. Scaling estimates indicate that a 1 
reduction of about 15% in the range corresponding to a tissue ind 
! , dose of 500 rad can be expected in arctic conditions. This is 9 
about the same uncertainty that exists in current inodeling of "* 


s 
a 


this radiation component, 


x, aes Differences in the fallout in arctic conditions could 
arise in several ways: 


ek, 


differences in the induced activity, 


differences in the size of the particles the active particles 4 
are attached to, differences in the debris cloud development NK 

and dispersal due to the meteorological conditions, and for ~ 

, : ri : P \) 
< underwater bursts differences in the radioactive material ¥ 

ejection into the air due to the ice cover. : 
‘Sy 
The induced activity in bursts over arctic soils will me. 

eo probably be essentially the same as in temperate climates since * 
there is in general the same range of soil types there. For . 

bursts over sea water or sea ice the induced activity is much 
an 

less than over ground and for bursts over fresh water, snow, ‘is 
or ice the induced activity is zero. Thus, in many situations 4 
in the arctic the residual radiation source is due only to the Nl 
fission yield of the weapon and no induced activity from the AY 

; j ) 
ao thermonuclear component will exist, a 
3 ” 
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For bursts over arctic soil the particle sizes which 
result in the rising debris cloud would not be expected to be 
different than existing in more temperate climates since the 
basic soil types are comparable. For bursts on snow, ice or 
sea water, however, one would expect considerably different 
debris cloud characteristics. One probably would not expect 
the arctic case to be much different from a sea burst in 
temperate areas. 


Debris cloud develupment and dispersal has not been 
considered for arctic conditions. One might expect somewhat 
different development because of the different density and 
temperature profiles. Winds are not significantly different 
in the Arctic except perhaps being more variable; so no signif- 
icant differences in fallout predictions would be expected 
except an increase in the uncertainties of such already uncertain 
predictions. 


Relative scavenging efficiencies of snow and water have 


= measured and analyzed with conflicting results. Some studies 


indicate a much larger scavenging efficiency for snow than water 
while other studies indicate no difference. The Arctic has a 
much lower precipitation rate than most temperate areas so that 
one might expect less of the activity to be scavenged and might 
expect therefore a more wide ranging and less intense fallout 
pattern that might occur in temporate areas if precipitation 
occurs. No studies have been made of induced precipitation by 
nuclear bursts in the Arctic. 


aa The major uncertainties in predicting the effects of 


nuclear radiation from an underwater burst result from a lack 
of knowledge of the amount of explosive energy that is required 
to break through an ice layer and that is therefore lost as far 
as the development of surface and above-surface phenomena are 
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concerned. This in turn leads to uncertainty concerning the 
amount of radioactive material ejected to the atmosphere, the 
extent of its initial dispersion, and, in the cases of very 

small yields or great depths of burst, whether the ice will 
contain the effects of the detonation so that there will be no 
atmospheric phenomena. It is considered that available studies, 
though unverified experimentally, are adequate to the understand- 
ing of the effects of radioactive products freezing on exposed 
surfaces and the probable effects of the arctic environment on 
evolutior of the radioactive pool. 


fF Studies have been made of the accumulation of activity 


on ships in icing conditicns. This does not seen to be a very 
important mechanism of damage. 


5.5.2  ] Recommendations 


The effects noted in the prompt radiation environments 
were not very large but covld be of significance for specific 
systems. It is recommended that currently available air transport 
results be scaled to provide isofluence and isodcse profiles 
for neutrons, gamma rays and x-rays from selected weapon classes 
as a function of burst altitude. This could be done for the _ 
Standar’ arctic conditions as well as for other extreme condi- 
tions which can exist as indicated in Section 1.2. 


F | These predictions should be incorporated into the 


appropriate chapters of EM-1 and perhaps could be a part of a 
more general section relating to the effects of atmospheric 
departures from standard on radiation transport. 


The effects on the early fission product dose should 
be determined for a few selected cases including possible 
fireball and cloud development changes. These calculations 
should be used to indicate scaling procedures so that inexpensive 
predictions can be prepared for a range of practical cases. 
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i Fallout prediction is at best very uncertain. The 


additional complications introduced by arctic conditions for 

which the U.S. will never have empirical data make the predictions 
even more untrustworthy. Since fallout is usually treated as a 
collateral damage mechanism ir military situations, there may 

be no need to have accurate prediction techniques in the Arctic. 


Basic studies are required to specify the size dis- 
tribution of the debris particles in the Arctic. Resolution of 
the discrepancies that exist in the analyses of the relative 


scavenging efficiencies of Snow and water is necessary before 
predictions of the Fallout under arctic conditions can be 


nade. 


Computer models exist that could be used to compare 
the fallout from arctic and temperate climates. These models 
require as inputs such information-as the debris cloud develop- 
ment, loading and particie size distributions, wind patterns 
as a function of cltitude, precipitation patterns and rates, 
scavenging efficiencies, and particle diffusion character- 
istics. It is recommended that preliminary studies in these 
various areas be performed to identify the maximum differences 
that might exist in these parameters between the arctic and 
temperate climates. Predictions of the fallout using the minimum 
parameter differences should be made. If militarily Significant 
differences occur between the arctic and temperate conditions, 
then additional research may be required in specific areas. 


Sar Fallout predictions from underwater bursts are very 
uncertain. The ultimate destiny of the radioactive materials 


Ae | - 


for variows DOB is very uncertain and specifically the fraction ) 
that appears above the surface to contribute to fallout is S 
unknown. It has been conjectured that the forces associated | Xs 
with the range of yields and depths of burst that are likely a 
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to 2 of interest for underwater bursts are so great that the 


energy loss in breaking through ice would have minimal effect 

on the development of surface and above-surface phenomena. 

Since the fallout is in general less than that expected over 
land, there may be some question about its importance except 

in very specific cases involving nearby surface ships. If hydro- 
@ynamic calculations are made of the bubble development and shock 
interactions with the water and ice layers, tracer particles 
should be introduced in an attempt to determine the distribution 
of the radioactive particles for various DOB conditions. 
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SECTION 6 
COMMUNICATIONS AND EMP 


ad This study was nominally limited to low altitude 

ursts and effects, and no large effort was to be expended on 
high altiteade effects. Ouring the literature seacches one study 
(Jordano, et al, 1976) was found which specifically addressed the 
latitude dependence on HF absorption and the effects of model 


atmosohere differences on debris cloud development. This report 
was reviewed and a summary of the results is included. 


Some chanyes are expected in the EMP on the surface 
due to the differences in the magnetic field intensities and 
direction. The variation in EMP SMILE diagrams for various 
latitudes and various longitudes in the polar region are given. 


6.1 OP xcric Environmentai Differences 


me The different profiles for the atmospheric parameters 
(density, pressure and temperature) for the arctic region “can 


effect the debris cloud development and stabilization altitude. 
The delayed gamma-ray source function may then be ditferent, which 
can cause differences in tne fonization levels and attenuation 
peoperties of the atmosphere for electromaynetic wave propagation. 
The reaction rate constant. that determine the sustained ioniza~ 
tion levels are a function of tne temperature and particle concen- 
trations. The concentration of minor species can be important 

tn determining deionization rates and may he altered at high 
latitudes due to the differences in energetic particle effects 
noted in this region at high altitudes. 
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The greater intensity of the magnetic field in the 


J 

» 

! 

) 

, 

) 

| 

; 
polar areas will increase the magnitude of EMP in the region. 
The direction is more vertical than at lower latitudes which 
will affect the relative magnitudes of the horizontal and 
vertical components of the EMP and may affect the coupling of 


the EMP into targets. 


The coupling of energy into structures and cables 
buried under the ice and snow may be affected because of the 
differences between the conductivity, capacitivity and 
permeability of ice and snow and of more typical soil materials. 


6.2 GF sccruation of HF Communication 


The changes in HF absorption from a near surface burst 
due to the different profiles of density, pressure ard temperatures 
has been considered by Jordano et al (1978). The calculations 

were done by defining high latitude atmospheric models, incorpor- 
ating the models iz existing communication codes, and comparing 

the effect on HF communication links passing through the D region. 
The eftect of the atmospheric differences on the debris dynamics 
was also considered. The influence of the atmospheric parameters 
on the deionization kinetics wat considered, but the effect of 
differing concentrations of the minor species was not included. 
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. for the July and January 60° N models are shown in Figure 6-1. 
The July and January extreme profiles were defined by adding 

| . to or subtracting from the 60°N profiles a component represent-~ 
| ing the diurnal variation plus a two-sigma random variation in 


: 5 #4 
such a manner as to increase the variation from the mean standard Pt. 
. ; ; ist 
profile. The circles plotted for altitudes below 30 km represent 7, 
Me 
the 75°N January temperature profile described in Section 1.2 and if 
: are seen to agree with the defined January extreme model. The b 
WEPH VI/ROSCOE system defines the pressure and density profiles ai 
from the temperature profile using hydrostatic equilibrium and A 
x 
match neasureé mean molecular weights. ee 
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It iS recommended thet additional studies be made on 
os detcis cynamics for Arctic atmospheres. The results of these 
studies have inplications in the determination of fallout under 
Arctic conditions as well as in the effect on communication black- 


out. 


The concentration of minor species in the high aititudes 
should be considered and their effect on the deionization kinetics 
should be determined. Prediction of the communication blackout 
expected by bursts in the Arctic including a11 of these effects 
should be made for a wide range of frequencies. 
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| i The horizontal and vertical components of EMP 
shoul’ be determined for the polar region. Changes in coupling 
to surface based systems should be considered. Predictions 


should be made of the coupling of the EMP to cables and facilities 
buried in frozen ground or covered with snow. 
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